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NOTATION 

A list of the symbols used throughout this document and their definitions is provided 
below for convenience. 

Roman Symbols 

a . . . speed of sound 

c p . . . gas specific heat at constant pressure 
c. v ... gas specific heat at constant volume 
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e . . . total internal energy 
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p . . . pressure 

r . . . radius or radial coordinate 
t . . . time 

v x . . . velocity in the Cartesian coordinate system x direction 

v y . . . velocity in the Cartesian coordinate system y direction 

v z . . . velocity in the Cartesian coordinate system z direction 

v r . . . velocity in the cylindrical coordinate system radial direction 

v 0 ... velocity in the cylindrical coordinate system circumferential direction 

w re t . . . relative velocity in the circumferential direction (= Vo - ru>) 

x . . . Cartesian coordinate system coordinate 

y . .. Cartesian coordinate system coordinate 

2 . . . Cartesian coordinate system coordinate 

ADPAC07 . . . Advanced Ducted Propfan Analysis Code Version 07 
ADS PIN . . . ADPAC post processing program 
ASCII . . . American Standard Code for Information Interchange 
CFL . . . Courant-Freidrichs-Lewv number (At/Nt max, stable) 

D . . . diameter 

DOC . . . Direct Operating Costs 
EEE . . . Energy Efficienct Engine 

F . .. i coordinate direction flux vector 
FPS . . . Flight Propulsion System 

G . .. j coordinate direction flux vector 

GRIDGEN . . . Multiple block general purpose mesh generation system 

H . .. k coordinate direction flux vector 
Piloted ■ ■ ■ total enthalpy 

HP . . . high pressure 

ICLS . . . Integrated Core/Low Spool 

ISA . . . Industry Standard Atmosphere 

K ... cylindrical coordinate system source vector 

L. .. reference length 
LP . . . low pressure 

LPS . . . low pressure subsystem 

M . . . Mach number 

MFLOP . . . million floating point operations per second (CPU speed factor) 

N . . . Number of blades 

Q . . . vector of conserved variables 

R . . . gas constant or residual or maximum radius 

IZ. .. turbulent Reynolds number 
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Re. . . Reynolds Number 

P ... turbulence kinetic energy production term 
Pr . . . gas Prandtl Number 

S . . . surface area normal vector 

SDBLIB . . . Scientific DataBase Library (binary file I/O routines) 

T . . . Temperature 

U . .. Freestream velocity (units of length/time) 

V . . . volume 


Greek Symbols 


7 . . . specific heat ratio 
A . . . calculation increment 
e . . . turbulence dissipation parameter 
V . . . gradient vector operator 
u . . . vorticity 
p . . . density 

/i . . . coefficient of viscosity 
t . . . fictitous time or shear stress 
Iljj . . . fluid stress tensor 


Subscripts 


[ ]i . . . inlet value 
[ ] 2 . . . exit value 

[ ]ax • • • pertaining to the axial (;c) cylindrical coordinate 

[ ] coarse • • ■ coarse mesh value 

[ ]e / fective ■ ■ ■ effective value 

[ ]/ine ■ • • fine mesh value 

[ ]freeslream ■ ■ • freestream value 

[ • • • grid point index of variable 

[ ]iaminar ■ ■ ■ laminar flow value 

[ ]max ■ • • maximum value 

[ } m in • • • minimum value 

[ ]nearwait ■ • • near wall value 

[ ] non— dimensional • • • non-dimensional value 

[ ]r • • • pertaining to the radial (r) cylindrical coordinate 

[ }ref ■ ■ ■ reference value 

[ ] stable * • • value implied by linear stability 

[](... turbulent flow value 

[ ]to£a/ * • • total (stagnation) value 

[ ^turbulent ■ • • turbulent flow value 

[ ]waii ■ ■ ■ value at the wall 
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[ jj, . . . pertaining to the y Cartesian coordinate 
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[ ] n . . . Time step index 
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[ ] . . . Time-averaged variable 

[ ] . . . Density- weighted time-averaged variable [ ] . . . Vector variable 
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Chapter 1 
SUMMARY 


The objective of this project is to provide the capability to analyze the aerodynamic 
performance of the complete low pressure subsystem (LPS) of the Energy Efficient 
Engine (EEE). The analyses were performed using three-dimensional Xavier-Stokes 
numerical models employing advanced clustered processor computing platforms. The 
analysis evaluates the impact of steady aerodynamic interaction effects between the 
components of the LPS at design and off-design operating conditions. Mechanical cou- 
pling is provided by adjusting the rotational speed of common shaft-mounted compo- 
nents until a power balance is achieved. The Xavier-Stokes modeling of the complete 
low pressure subsystem provides critical knowledge of component aero/mechanical in- 
teractions that previously were unknown to the designer until after hardware testing. 
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Chapter 2 


INTRODUCTION 


2.1 Description of Engine Balance Design Problem 

Competitive market conditions in the gas turbine industry have placed stringent 
demands on engine manufacturers to respond to customer requirements with efficient, 
cost effective products with significant reductions in development time. During the 
engine development period, component efficiencies often fall short of desired goals by 
significant margins. The engine cycle rebalance which results causes other components 
to operate at non-optimal (off-design) flow conditions, further reducing efficiency 
and complicating the identification of the original source of inefficiency. Expensive, 
multiple build rig testing, representing a major portion of the overall development 
cost, has, in the past, been required to balance component performance and optimize 
the engine system design. 


2.2 NASA Programs Addressing Advanced Engine 
Design 


Effoits to attack the problems associated with aircraft gas turbine engine develop- 
ment have been addressed through several NASA Programs. The Advanced Subsonic 
Technology (AST) program specifically supports technology development to improve 
the performance of subsonic aircraft, both in flight characteristics and propulsion. 
The High Performance Computing and Communication (HPCC) Program and more 
specifically, the Computational Aerosciences (CAS) Project are directed to acceler- 
ate the availability of high performance computing technology for use by the U.S. 
aerospace community. Finally, the Numerical Propulsion System Simulation (XPSS) 
project represents an interdisciplinary program to unite the various disciplines used 
in gas turbine engine design. A primary goal of the XPSS program is to numerically 
solve the entire flow through a realistic gas turbine engine using high fidelity compu- 
tational tools. Additional details of these programs and their relation to the present 
work is given in the sections below. 
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2.2.1 NASA High Performance Computing and Communi- 
cations (HCCP) Program [1] 

NASA is a key participant in the Federal High Performance Computing and Com- 
munications (HPCC) Program. As a key participant of the Federal Program, the 
primary purpose of NASA’s HPCC Program is to extend U.S. technological leader- 
ship in high-performance computing and communications for the benefit of NASA 
stakeholders: the U.S. aeronautics, earth and space sciences, and spaceboine ieseaich 
communities. As international competition intensifies and as scientists push back the 
frontiers of knowledge, leading-edge computational science is more important than 
ever. 

The NASA Program is structured to contribute to broad Federal efforts while 
addressing agency-specific computational problems called Grand Challenges. NASA 
provides resources to develop tools to solve Grand Challenges in four HPCC project 
areas: Computational Aerosciences (CAS), Earth and Space Sciences, Remote Ex- 
ploration and Experimentation, and Information Infrstructure Technology and Ap- 
plications. The NASA Research and Education Network (NREN) also supports the 
four projects. 

2.2.2 The Computational Aerosciences (CAS) Project [2] 

CAS is a computer science-related program oriented around the needs of the aero- 
science community. The CAS project [2] directly supports other NASA aeronautics 
programs and is driven by the needs of the aeronautics industry. The CAS goal is to. 

“Accelerate development and availability of high-performance computing technol- 
ogy of use to the U.S. aerospace community, to facilitate adoption and use of this 
technology by the U.S. aerospace industry, and to hasten emergence of a viable com- 
mercial market for hardware and software vendors to exploit this lead. 

The science and engineering requirements inherent in the NASA Grand Chal- 
lenge applications like aeronautics require orders of magnitude improvement in high- 
performance computing and networking capabilities over the capabilities that existed 
at the beginning of the NASA HPCC Program in FY1992. Without an accelerated 
development program, this level of improvement may not be available for many years. 

CAS has traditionally been oriented around the longer term thrust of the ex- 
ploration of future high-end supercomputing for aerospace needs - extreme high- 
performance computing (TeraFLOPS). As a result of increased interactions with in- 
dustry, CAS has added research efforts in a new direction - the use of networked 
workstations in the design environment. Networked workstations is a shoitei teim 
thrust oriented around the effective use of current generation computing hardware to 
reduce costs. 

The goals of the Grand Challenges in Computational Aerosciences are to: 

1 Provide focus for the entire HPCC Program by providing requirements for Testbeds 

and Networks and Systems Software 


NASA/CR— 1998-206597 


4 



2 Assuie relevance of HPCC Program to the U.S. aerospace community by provid- 
ing base computational technology for multi-disciplinary analysis and design of 
aerospace vehicles and propulsion systems on HPCCP platforms, demonstrating 
superiority of HPCCP systems for solution of relevant aerospace problems, assist- 
ing U.S. aerospace industry in implementing HPCCP technologies within their 
organizations and realizing real improvement in their design cycle process and 
final products. 

The goal of the Computing Testbeds of the CAS Project is to support the ac- 
celeiated development of cost-effective, high performance computing machinery from 
domestic vendors in order to benefit the aerospace industry through the: 

1 . Creation of “beta-test” prototype computing facilities scalable to TeraFLOPS, and 
through the evaluation of the functionality and robustness of associated system 
software. 

2. Creation of prototype networked workstation clusters, that are representative of 
existing clusters in aerospace companies, in order to provide the environment to 
develop and test the software necessary to make clusters an alternative to the 
traditional supercomputer. 

The goal of CAS System Software research is to identify, define, and provide the 
critical software and tools not available from others sources that will enable the ef- 
fective utilization of networked, heterogeneous, high performance scalable computing 
environments. 

2.2.3 The Numerical Propulsion System Simulator (NPSS) [3] 

The Numerical Propulsion System Simulator (XPSS) [3] is an interdisciplinary project 
to unite the various disciplines used in gas turbine engine design. The project is 
coordinated by NASA Lewis and is designed to bring together different groups and 
codes in order to create a system to engineer advanced jet engine designs. XPSS can 
be lepiesented by thiee main concepts: “zooming,'’ “ coupling, ’ and “integration.” 
Zooming enables the simulation of complete engine systems at a level of analysis 
required by the physics. Coupling refers to the joining of the various disciplines in a 
single analysis. Integration refers to the integration of the various engine component 
simulations. 


2.3 Integration of NASA Program Objectives and 
the Engine Balance Design Problem 

A jet engine can be characterized by a number of different components working to- 
getliei very efficiently within a range of demanding operating conditions. Several of 
these engine components are sensitive to interactions with neighboring components. 
F°i example, the efficiency of the compressor is very sensitive to steady inlet and out- 
let conditions, outlet pressure fluctuations can unstart a supersonic inlet and expel the 
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shock, substantially increasing drag and reducing engine performance. Consequently, 
during the design process it is important to consider not only isolated components 
but the engine as a system of components which influence each other. 

Historically, the design process has started with a study of the complete proposed 
engine using performance maps and one-dimensional analysis. Then individual engine 
components are simulated and designed in detail by component design teams. Some 
engine components are designed by the airframe manufacturer. These results improve 
the performance maps and one-dimensional analysis, which helps address component 
interactions. These components are experimentally tested in isolation, progressively 
integrated, and adjusted to finalize the engine design. 

Component design teams depend on numerical analysis techniques to achieve the 
best performance. Streamline curvature methods continue to be extensively used to 
analyze multistage turbomachinery. More recently, the trend has been to apply ad- 
vanced 2-D and 3-D numerical techniques [4] to engine components to understand 
the details of their operation in isolation. These applications range from quasi-tin ee- 
dimensional blade calculations which predict the behavior of a transonic blade to 
multistage compressor calculations which simulate the behavior of transonic com- 
pressors to simulation of nacelles and combustor chemistry. Multistage analyses for 
turbomachinery are also becoming increasingly more valuable [5], [6]. 

These advanced component analysis techniques do not systematically account for 
inter-component interactions. Multistage analyses may someday provide adequate 
representation of interaction effects between blade rows in an axial compressor, for 
example, but do not presently provide information related to inter-component inter- 
actions (HP/LP turbine systems, e.g.). 

One goal of XPSS is to create a system which will allow these individual component 
codes to be coupled to create a full engine simulation system. This system will 
allow analysis at different levels of accuracy by coupling codes of all levels from 1- 
D models to full 3-D computational fluid dynamics codes. This system would then 
allow the design engineer to “zoom'* between levels of detail, while still pioviding 
some indication of the overall system interaction effects. 

2.4 NPSS and the Energy Efficient Engine (EEE) 

Several examples of multistage turbomachinery aerodynamic performance prediction 
techniques [7], [5], [6] exist which demonstrate the viability of large scale simulation 
in the gas turbine engine design environment. Unfortunately, many of these models 
have only explored aerodynamic interaction effects for a specific subcomponent of 
an engine (HP compressor, or fan section in isolation, e.g.). Improvements in the 
power and availability of high speed processors, and a streamlining of the problems 
associated with large scale simulation data management has afforded the opportunity 
to perform large scale simulations of coupled subsystem components, and peihaps 
even an entire engine. 

In 1976 NASA initiated the Aircraft Energy Efficiency (ACEE) Program to assist 
in the development of technology for more fuel-efficient aircraft for commeicial aiiline 
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use. The Energy Efficient Engine (EEE) Project of the ACEE program was intended 
to lay the advanced technology foundation for a new generation of turbofan engines. 
This project, planned as a seven-year cooperative government-industry effort, was 
aimed at developing and demonstrating advanced component and systems technolo- 
gies for engines that could be introduced into airline service by the late 1980’s or 
early 1990’s. In addition to fuel savings, these new engines offered the potential for 
being economically atti active to the airline users and environment all v acceptable. 

The goals of the EEE program were: 


• 12% reduction in installed specific fuel consumption compared to a CF6-50C at 
maximum cruise thrust, M=0.8 at 35,000 ft ISA 

• Comply with I'AR 36 (1978) with provisions for growth 

• Comply with EPA Proposed (1981) Standards for new engines 

• 50% reduction in the rate of performance deterioration in-service as compared to 
the CF6-50C 

The EEE Program consisted of four major technical tasks structured as follows: 

1 Piopulsion System Analysis, Design, and Integration (Establish the component 
design and performance requirements for future tasks). 

2 Component Analysis, Design, and Integration (Design, fabrication, test, and post- 
test analysis of the components and supporting technologies). 

3 Core Test (Design, fabrication, test, and post-test, analysis of the core test vehicle 
(HP compressor, combustor, and HP turbine). 

4 Integrated Core/Low Spool (ICLS) Test (Design, fabrication, test, and post-test 
analysis of the ICLS turbofan ground test vehicle). 

The data obtained during these tasks yielded insight into the evaluation of core 
components operating in isolation and in the engine environment, and also permitted 
accuiate measuiements of important internal conditions which would be impractical 
in a complete turbofan engine. 

The EEE provides a natural vehicle for the type of large scale simulation planned 
for this study due to the availability of both subcomponent test rig data, as well as 
fully coupled, assembled engine test data. 


2.5 Objectives of the Present Study 

This project represents a consolidation of industry goals, NASA vision, and the grow- 
ing maturity of computational tools for predicting gas turbine engine flow physics. 
The overall objective of this project is to provide the capability to analyze the aerody- 
namics in the complete low pressure subsystem (LPS) of the Energy Efficient Engine 
(EEE) using three-dimensional Xavier-Stokes numerical models. The analysis eval- 
uates the impact of steady aerodynamic interaction effects between the components 
of the LPS at the design and at off-design operating conditions. The LPS modeling 
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capability will be integrated into the XPSS and be available as an option to the de- 
signer. The long range goal of the LPS modeling project, and XPSS, is to provide a 
tool that can significantly reduce the design, development, and ceitification time of 
gas turbine engines. 

The approach for creating the LPS model is to select a validated Xavier-Stokes 
(X-S) analysis code for accuracy and to minimize turnaround time. The EEE LPS 
model was developed from the geometric components in the LPS including: external 
flow, nacelle, inlet, fan blades, bifurcated bypass and core inlet, bypass vanes, core 
inlet guide vanes, quarter height booster stage, low pressure turbine blades, mixer, and 
exhaust nozzle. Initially, the engine core components were modeled using appropriate 
boundary conditions. At a later stage in the program, an engine cycle peifoimance 
deck was used to set the core operating conditions for the analysis. The complete LPS 
analysis was constructed following a verification of the performance of the individual 
subcomponents in the LPS. The X-S analysis of the fully coupled LPS enabled a 
torque balance on the low pressure spool at quasi-steadv state opeiating conditions. 
This study was divided into five major milestone areas: 

1 Geometry Definition: Detailed geometry definitions of the components of the 
Energy Efficient Engine primary gas flowpaths were assembled. 

2 Mesh Generation: Geometry definitions described above were employed to de- 
velop discrete mesh systems suitable for CFD analysis. 

3 Component validation study: Block components of the LP and HP subsystems 
were analvzed using CFD tools to verify the accuracy of the geometiv definitions, 
and to validate the CFD analysis with available rig test data. 

4 LP Subsystem Analysis: Various components were assembled to form the dis- 
crete representation of the LP Subsystem, and a quasi-steady CFD analysis was 
applied to predict both the aerodynamic and mechanical coupling of the LP Sub- 
system. 

5 Core Cycle Specification: An engine cycle performance model was coupled 
with the 3-D CFD analysis to represent the operating parameters for the engine 
core in the LP Subsystem Analysis. 

Each of the five milestone topics are described in more detail in the chapters which 
follow. The ultimate objective of this study was to develop a simulation capability for 
the LP Subsystem of modern high bypass ratio turbofan engines which would address 
the goals of the XASA XPSS program. 
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Chapter 3 

DESCRIPTION OF THE 
ENERGY EFFICIENT ENGINE 

The Energy Efficient Engine (EEE) program [8]-[103] was developed to create fuel 
saving technologies foi transport aircraft engines which would be introduced into 
service in the late 1980’s and 1990’s. The EEE development cycle included candi- 
date engines from two manufacturers: Pratt & Whitney and General Electric. Both 
manufacturers designed and tested various components as part of the technology 
demonstrations necessary to validate the final engine designs. The General Electric 
design was selected for engine testing, and included separate tests of the core [92] and 
integrated core/low spool (ICLS) [99] configurations. In the course of this discussion, 
iefeience will also be made to the flight propulsion system (FPS), which is essentially 
the integrated core/low spool with a flight-ready nacelle and inlet, rather than the 
bellmouth arrangement used in the static propulsion tests. 

An illustration of the General Electric EEE flight propulsion system layout and 
some of the design feat ures is given in Figure 3.1. A table of cycle characteristics for 
the EEE FPS are given in Table 3.1. Based on corrections to test data, the flight 
propulsion system was projected to have a thrust specific fuel consumption of 0. 551 
lbm/hr/lbf at the maximum cruise design point (35,000 ft. ISA). The ICLS achieved 
a static corrected take-off thrust of 37,415 lbf. 

An illustration of the major subcomponent arrangement for the EEE Low Pres- 
suie spool is given in Figure 3.2. The elimination of the high pressure spool from 
the proposed analysis is illustrated in the sample numerical mesh system depicted in 
this figure. The analysis of the LP spool entails considerable detail in managing both 
aerodynamic and mechanical performance of the fan section, LP turbine, exhaust 
mixei , and inlet/nozzle/ external flowfields. Detailed presentations of design expecta- 
tion, test measurement , and CFD prediction of the individual component performance 
data are presented in the following chapters. Coupled analysis of the HP/LP spool 
systems are presented in the final chapters dealing with the LP subsystem analysis. 
Numerous sources of information related to the EEE program are provided in the 
reference section of this report. 
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Figure 3.1: Energy Efficient Engine layout and design features. 


Cycle Pressure Ratio at Max Climb* 

Bypass Ratio at Max Climb* 

Fan Pressure Ratio at Max Climb* 

Turbine Rotor Inlet Temperature at 
Static Warm Day2 Take-off Power 

Specific Fuel Consumption at Max Cruise 3 , 
Bare Engine 

Specific Fuel Consumption at Max Cruise 3 , 
Installed Engine 


38 
6.8 
1.65 
2450 F 

0.542 lbm/(lbf-hr) 
0.564 lbm/(lbf-hr) 


1 Max Climb is the aerodynamic design point, M=0.8, 

35,000 ft, standard day +18 F. 

2 Sea level static warm day refers to a standard 59 F. 

3 Max cruise is the performance evaluation design point, 

M=0.8, 35,000 ft., standard day. 

Table 3.1: Energy Efficient Engine Flight Propulsion System cycle characteristics. 


N AS A/CR— 1 998- 206597 


10 





High Pressure Core Engine from Cycle Analysis 


3-D 

Navier-Stokes 

Analysis 


Figure 3.2: Energy Efficient Engine component description and CFD mesh represen- 
tation. 
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Chapter 4 

GEOMETRY MODELING 


Detailed geometry for the EEE model was extracted from the NASA Energy Efficient 
Engine Master Geometry Database. This database was developed specifically for 
XPSS-related applications which employ the EEE design for demonstration. The 
database consists of NASA IGES curve-based and surface-based entities describing 
the major components of the engine core and bypass gas flowpaths. Exact geometric 
definitions of the EEE LPS are employed, with the exception of the outer nacelle and 
inlet, which have been designed consistent with the Energy Efficient Engine design 
philosophy in order to take the place of the test rig bellmouth. A picture of the 
Energy Efficient Engine test rig hardware is given in Figure 4.1. 

The geometry database consists of individual elements (separate blade rows, for 
example) as well as “assembled” systems, which consist of more complete coupled 
collections of components. Familiarization with the geometry database package was 
facilitated by using the PATRAX [104] geometry modeling software package. A PA- 
TRAX representation of the EEE hardware geometry is given in Figure 4.2. 

Certain enhancements to the database will almost certainly be required for this 
type of geometry definition to be useful during the engine design process. The ability 
to reset blade stagger angle, for example, is a common operation in gas turbine engine 
design and test, but is still an overly complex operation with the current database. 
The EEE HP compressor employs variable geometry on several stators, for example, 
and rig test results were obtained with various stator settings which are difficult to 
reproduce in the current database arrangement. The blade restagger capability will 
require definition of the rotation axis in the database, and specific built-in stator reset 
schedules could be imported as “off-design” geometry definitions. 

The current database does not contain any indication of rotor tip clearances, 
and this was essentially approximated from experience during most of the course 
of this study. The database should include at least a reasonable approximation of 
rotor tip clearances (and other important clearance dimensions as additional geometry 
components are incorporated into the database). 

The flow in the primary gas flowpath of a modern turbine engines is complicated 
by the various networks of secondary flow systems for cooling, bleed, etc. Compres- 
sor flowpaths are affected by leakage flows through inner-banded stator seals, while 
turbine flows are complicated by the stepped, overlapped hub flowpath and inner 


NASA/CR— 1998-206597 


13 






Figure 4.1: Energy Efficient Engine test rig hardware. 



Figure 4.2: PATRAX representation of the Energy Efficient Engine Master Engine 
Geometry Database. 
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wheel purge flow. The contributions of these secondary flows in the engine play an 
important role in determining the overall performance of the machine. During the 
course of this study, the database flowpath representations were smooth walls with no 
representations of secondary system leakage flows. As experience is gained with large 
scale system aerodynamic analyses, and the analyses become more sophisticated, the 
influences of these systems must be included to accurately model the overall engine 
flow performance. 

Finally, consideration must be given to the overall operational life of the engine 
and the potential for performance degradation through component erosion and wear. 
Over time, clearances become larger, blade leading edges can become warped due 
to foreign object damage, and erosion, in general, alters the blade surface quality 
and even the blade profile. If the database is to be useful for the overall engine life 
performance analysis, then these effects must necessarily be incorporated in some 
manner. 

A UNIX tar archive listing of the geometric components included in the Master 
Engine Geometry Database for this study is included below for reference: 


Engine Assembly Components: 


rwxr-xr-x 2788/100 

May 22 

08:05 

1996 

rw-r — r — 2788/100 

May 

10 

13:30 

1996 

rw-r— r— 2788/100 

May 

10 

13:24 

1996 

rw-r— r— 2788/100 

May 

10 

13:40 

1996 

rw-r— r— 2788/100 May 10 
HP Compressor Components: 

13:26 

1996 

rwxr-xr-x 2788/100 

May 

22 

08:06 

1996 

rw-r — r — 2788/100 

May 

20 

13:04 

1996 

rw-r — r — 2788/100 

May 

20 

13:04 

1996 

rw-r— r— 2788/100 

May 

20 

13:04 

1996 

rw-r — r — 2788/100 

May 

20 

13:04 

1996 

rw-r— r— 2788/100 

May 

20 

13:04 

1996 

rw-r— r— 2788/100 

May 

20 

13:04 

1996 

rw-r— r— 2788/100 

May 

20 

13:04 

1996 

rw-r— r— 2788/100 

May 

20 

13:04 

1996 

rw-r — r — 2788/100 

May 

20 

13:05 

1996 

rw-r — r — 2788/100 

May 

20 

13:05 

1996 

rw-r— r— 2788/100 

May 

20 

13:05 

1996 

rw-r— r— 2788/100 

May 

20 

13:05 

1996 

rw-r— r— 2788/100 

May 

20 

13:05 

1996 

rw-r — r — 2788/100 

May 

20 

13:05 

1996 

rw-r— r— 2788/100 

May 

20 

13:05 

1996 


assembly/ 

assembly/ eee_eng_3d_symall .nigs 
assembly/ eee_eng_3d_symf low . nigs 
assembly/ eee_eng_3d_unsymall .nigs 
assembly/eee_eng_3d_unsymf low . nigs 


comp_ schmidt/ 

comp_schmidt/cmpr_igv_srf .nigs 
comp_schmidt/cmpr_rotor 10_srf . nigs 
comp_schmidt/cmpr_rotorl_srf .nigs 
comp_schmidt/ cmpr_rotor2_srf . nigs 
comp_schmidt/cmpr_rotor3_srf .nigs 
comp_schmidt/cmpr_rotor4_srf .nigs 
comp_schmidt/cmpr_rotor5_srf .nigs 
comp_schmidt/ cmpr_rotor6_srf . nigs 
comp_schmidt/ cmpr_rotor7_srf . nigs 
comp_schmidt/cmpr_rotor8_srf . nigs 
comp_schmidt/cmpr_rotor9_srf . nigs 
comp_schmidt/ cmpr_statorlO_srf . nigs 
comp_schmidt/cmpr_statorl_srf .nigs 
comp_schmidt/cmpr_stator2_srf . nigs 
comp_schmidt/cmpr_stator3_srf . nigs 
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rw-r — r — 2788/100 May 20 13:05 1996 comp_schmidt/cmpr_stator4_srf .nigs 

rw-r — r — 2788/100 May 20 13:05 1996 comp_schmidt/cmpr_stator5_srf .nigs 

rw-r — r — 2788/100 May 20 13:05 1996 comp_schmidt/cmpr_stator6_srf .nigs 

rw-r — r — 2788/100 May 20 13:05 1996 comp_schmidt/cmpr_stator7_srf .nigs 

rw-r — r — 2788/100 May 20 13:05 1996 comp_schmidt/cmpr_stator8_srf .nigs 

rw-r — r — 2788/100 May 20 13:05 1996 comp_schmidt/cmpr_stator9_srf .nigs 

HP Turbine Components: 

rwxr-xr-x 2788/100 May 22 08:06 1996 hpt/ 

rw-r — r — 2788/100 May 10 08:48 1996 hpt/hpt_rotorl_srf .nigs 

rw-r — r — 2788/100 May 10 08:50 1996 hpt/hpt_rotor2_srf .nigs 

rw-r — r — 2788/100 May 10 08:50 1996 hpt/hpt_statorl_srf .nigs 

rw-r — r — 2788/100 May 10 08:51 1996 hpt/hpt_stator2_srf .nigs 

LP Turbine Components: 

rwxr-xr-x 2788/100 May 22 08:06 1996 lpt/ 

rw-r — r — 2788/100 May 10 09:04 1996 lpt/lpt_rotorl_srf .nigs 

rw-r— r— 2788/100 May 10 09:04 1996 lpt/lpt_rotor2_srf .nigs 

rw-r— r— 2788/100 May 10 09:04 1996 lpt/lpt_rotor3_srf .nigs 

rw-r — r — 2788/100 May 10 09:04 1996 lpt/lpt_rotor4_srf .nigs 

rw-r — r — 2788/100 May 10 09:04 1996 lpt/lpt_rotor5_srf .nigs 

rw-r — r — 2788/100 May 10 09:04 1996 lpt/lpt_statorl_srf .nigs 

rw-r — r — 2788/100 May 10 09:05 1996 lpt/lpt_stator2_srf .nigs 

rw-r — r — 2788/100 May 10 09:05 1996 lpt/lpt_stator3_srf .nigs 

rw-r — r — 2788/100 May 10 09:05 1996 lpt/lpt_stator4_srf .nigs 

rw-r — r — 2788/100 May 10 09:05 1996 lpt/lpt_stator5_srf .nigs 

Fan + Quarter Height Booster Components: 

rwxr-xr-x 2788/100 May 22 08:07 1996 qtr.stage/ 

rw-r — r — 2788/100 May 14 10:01 1996 qtr_stage/booster_rotor_srf .nigs 

rw-r — r — 2788/100 May 14 10:01 1996 qtr_stage/booster_stator_srf .nigs 

rw-r — r — 2788/100 May 14 10:02 1996 qtr_stage/bypass_stator_srf .nigs 

rw-r — r — 2788/100 May 14 10:02 1996 qtr_stage/core_guide_vane_srf .nigs 

rw-r — r — 2788/100 May 14 10:02 1996 qtr_stage/fan_srf .nigs 

Lobed Exhaust Mixer Components : 

mixernew . igs 
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Chapter 5 


COMPONENT MESH 
GENERATION 


The development of the EEE/LPS CFD analysis requires numerical discretization of 
the Master Engine Geometry Database geometry definitions described in the previ- 
ous chapter. The nature of this discretization is defined by the requirements of the 
ADPAC CFD flow solver, which is described in more detail in the following chap- 
ter. Numerous meshing strategies are possible with the ADPAC code, the simplest of 
which is simply to use a single sheared H-type mesh for each blade row (see e.g. [5]). 
This meshing strategy also has the direct benefit that, the resulting mesh could also 
be used for other XPSS-related multistage turbomachinerv flow analyses such as AP- 
NASA. 

A key element of the meshing strategy in this project was to employ the Master 
Engine Geometry Database IGES entities directly in the grid generation process. 
Many mesh generation codes require discretized point data as input to define the 
geometry of interest. This discretized definition, and the subsequent interpolations 
which occur during the mesh generation process can lead to errors in the coordinates 
of the final mesh. One focus of the XPSS geometry definition has been to employ 
analytical definitions of geometric components in the form of IGES or XURBS-based 
entities. These analytical definitions would then form a consistent geometric database 
for all applications (aerodynamic, stress, heat transfer, etc.) and significantly reduce 
errors due to interpolations and interpretations of discrete point data. In order to 
address the mesh objectives described above, a procedure to generate meshes for the 
EEE LPS analysis directly from the NASA Energy Efficient Engine Master Engine 
Geometry Database was developed and is described in the paragraphs below. 

Detailed geometry for the EEE engine was extracted from the NASA Energy 
Efficient Engine Master Geometry Database. This database was developed specifically 
for XPSS-related applications which employ the EEE design for demonstration. The 
database consists of NASA IGES curve-based and surface- based entities describing 
the major components of the primary gas flowpath. Exact geometric definitions of 
the EEE LPS are employed, with the exception of the outer nacelle and inlet, which 
have been designed with the Energy Efficient Engine design philosophy in order to 
take the place of the test rig bellmouth. 
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Mesh Generation 


X(s,t) ={x, r, 0 } 



Minimum x NOT the Same as 
the Leading Edge 



Identify radial profiles of xmin and xj^ax from 
NURBS surface. 


© 


Use xmin and x MAX profiles to create 
meridional (x,r) point distribution (GRIDGEN) 




Sweep constant (x,r) points through NURBS to 
determine 9’s. Every point satisfies NURBS 
equation. 



Distribute points across pitch of blade 
holding constant near-wall spacing. 



Figure 5.1: Component mesh generation procedure for EEE LPS analysis. 


The construction of the numerical mesh system for each individual component is 
performed in a manner which permits a simple coupling of the component meshes 
for the complete LPS analysis. H-type computational meshes are employed for this 
purpose, although the analysis need not be limited in this fashion. A primary focus 
of the XPSS research is to employ a consistent geometry definition during all phases 
of the engine analysis. As such, a mesh generation strategy was developed whose 
only direct geometric input is the XASA-IGES based geometry of the Master Engine 
Geometry Database. A graphical illustration of the mechanics of the mesh generation 
procedure is given in Figure 5.1. 

The procedure is initiated by defining the exact geometric axial extents of the 
blade elements in the axisymmetric projection of the flowpath. This procedure was 
accomplished by interrogating the geometric elements for each individual blade row, 
and extracting the geometric leading and trailing edge outlines (in this sense, the 
geometric leading and trailing edges are represented by the minimum and maximum 
axial coordinate locations, respectively). In essence, the radial profiles of the blade 
minimum and maximum axial coordinates were extracted from the blade IGES surface 
definition. These new entities are themselves represented in GRIDGEN database 
segment format and are added to the geometry database. The SEARCH program 
was developed for this purpose. Source code for the SEARCH program is listed in 
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Figure 5.2: Illustration of the GRIDGEN user interface display for the meridional 
projection of the EEE fan section mesh system. 


the Appendix for this report. 

Once the blade row extents are defined, standard XASA-IGES capable mesh gen- 
eration schemes (GRIDGEX [1] was used for this exercise) can be employed to define 
the meridional projection of the H-type meshes. A snapshot of the gridgen user in- 
terface screen for the EEE fan section is illustrated in Figure 5.2. The blade leading 
and trailing edge elements define the positions of the blade rows in the axisymmetric 
projection, while the Master Engine Geometry Database flowpath definitions define 
the endwalls. The GRIDGEN program (which can read in the IGES entities as a 
geometry database) is then used to define the axial (a:) and radial (r) point distri- 
butions in the meridional projection. Typical mesh dimensions for the axisymmetric 
components of the meshes employed 49 points radially along the blade span, and 65 
points axially along the chord of the blade. 

X’ext, the (x,r) coordinate pairs from the meridional mesh projection are swept 
through the airfoil IGES surface definition to determine the blade surface circum- 
ferential (0) point distributions. The remaining points in the circumferential direc- 
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Figure 5.3: Axisymetric mesh projection for the EEE LP turbine. 


tion (between airfoils) are defined using a simple hyperbolic distribution routine (see 
e.g. [105]). The circumferential distributions were constructed to maintain a fixed, 
specified near wall spacing in the circumferential direction (see Appendix for coding 
details). 

A sample meridional mesh projection for the EEE LP turbine is given in Fig- 
ure 5.3. The blade outlines are visible due to the mesh cluster near the leading and 
trailing edges. An illustration of the 3-D mesh system for the same turbine is given 
in Figure 5.4. 
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Chapter 6 


ADPAC CODE DESCRIPTION 


The aerodynamic predictions for the cases described in this study were obtained 
using the ADPAC analysis code. The ADPAC code is a general purpose aerospace 
propulsion aerodynamic analysis tool which has undergone extensive development, 
testing, and verification [106]. Detailed code documentation is also available for the 
ADPAC program [107]. 

The ADPAC analysis solves a time-dependent form of the three-dimensional Reynolds- 
averaged Navier-Stokes equations using a proven time-marching numerical formula- 
tion. The numerical algorithm employs robust numerics based on a finite-volume, 
explicit Runge-Kutta time-marching solution algorithm derived from the develop- 
mental efforts of Jameson et al. [108], Adamczyk et al. [7], and Arnone et al. [109]. 
Steady-state flows are obtained as the time-independent limit of the time-marching 
procedure. Several steady-state convergence acceleration techniques (local time step- 
ping, implicit residual smoothing, and multigrid) are available to improve the overall 
computational efficiency of the analysis. A pseudo-time iterative implicit algorithm is 
available to permit large time steps for time-accurate flow predictions (see e.g. Mel- 
son et al., [110] ). A relatively standard implementation of the Baldwin and Lomax 
[111] turbulence model with wall functions was employed to compute the turbulent 
shear stresses and turbulent heat flux. 

An attractive feature of the ADPAC code is the versatility and generality of mesh 
systems upon which the analysis may be performed. The ADPAC code permits the 
use of a multiple-blocked mesh discretization which provides extreme flexibility for 
analyzing complex geometries. The block gridding technique enables the coupling 
of complex, multiple-region domains with common (non-overlapping) grid interface 
boundaries through specialized user-specified boundary condition procedures. An il- 
lustration of the wide variety of problems which have been analyzed using the ADPAC 
code is given on Figure 6.1. 

ADPAC supports coarse-grained computational parallelism via block boundary- 
specified message passing. Interprocessor communication is controlled by the Ap- 
plication Portable Parallel Library (APPL) [112] with optional programming layers 
using the Parallel Virtual Machine (PVM) [113] and Message Passing Interface (MPI) 
[114] communication protocols. Both serial and parallel computations were employed 
during this study utilizing a wide range of high speed processors, workstation clusters, 
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Propfans 



Ducted-Propfans 



UHB -Turbofans 



Sample ADPAC Applications 



Labyrinth Seal 



Linear Cascade 




Inlet Distortion Angle of Attack Rotor/Stator Interaction Boriv Forres 



Figure 6.1: Summary of variety of problems which can be analyzed using the ADPAC 
code. 
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and massively parallel computing platforms, depending on availability. 

Steady-state aerodynamic predictions for multistage turbomachinery are performed 
using a specialized boundary procedure known as a “mixing plane”. The mixing plane 
strategy was developed to permit numerical simulations based on only a single blade 
passage representation for each blade row, regardless of the differences in circumfer- 
ential spacing for each blade row. This simplification is afforded by circumferentially 
averaging data on either side of the interface between blade rows (the mixing plane), 
and then passing that information as a boundary condition to the neighboring blade 
row. This strategy is illustrated in Figure 6.2. 

Several formulations for the mixing plane have been tested in the development of 
the ADPAC code, including varying the choice of variables to integrate and the use 
of “non-reflecting” boundary procedures (see e.g. Saxer [115] ). A robust scheme 
results by simply averaging the conserved variables from the numerical scheme (p, pu, 
pv, pw , and pe) where p is density, u, v , w are the axial, radial, and circumferential 
velocity components, and e is the total internal energy. This scheme has the advantage 
of being numerically robust, conserves mass and momentum, and tends to preserve 
velocity triangle information across the interface plane more accurately than other 
approaches. A disadvantage of this scheme is that neither total pressure or total 
temperature are numerically conserved across the mixing surface. In practice, it was 
found that these conservation errors were detectable, but very small, and this scheme 
was therefore employed for the present set of calculations. 

A graphic illustration of a mixing plane analysis for a multistage compressor is 
shown in Figure 6.3. The mixing planes are represented by the circumferential lines 
approximately midway between blade rows. 

The solution procedure for the ADPAC analysis requires the definition of the nu- 
merical mesh, boundary conditions, and solution control input files. The meshing 
strategy for the EEE/LPS analysis was described in the previous chapter. The AD- 
PAC boundary data file were created through a combination of hand construction, and 
data provide by the PATCHFINDER ADPAC tool program. The PATCHFINDER 
program interrogates the mesh system, and through a rigorous coordinate search 
routine determines where neighboring mesh blocks share coordinates and outputs the 
specific ADPAC boundary specifications to couple the neighboring block aerodynamic 
solutions. The ADPAC input file is essentially constructed by hand, and determines 
solution specific parameters such as reference pressure and temperature, number of 
iterations, etc. 

The numerical solution proceeds with an initial flow specification from which the 
solution is advanced forward in time until the desired convergence criteria has been 
reached. The initial data is normally specified as a uniform flow, or may be read in 
as a “restart” of a previous existing solution. During the EEE/LPS simulations, the 
solution initialization procedure was complicated by the large range of pressures and 
temperatures encountered when doing large scale simulations of gas turbine engines. 
To ease the numerical problems with these wide variations in flow properties, the 
solution was initially started with very low pressure and temperature specifications 
in the boundary conditions, and was then iteratively restarted with subsequently 
larger values until the desired final conditions were achieved. The “full” multigrid 
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Geometry 



Computational Domain 



passage representation 
(or each blade row 


Boundary Treatment 



Phantom cells 
loaded with 
circumferentially- 
averaged data 
from neighboring 
blade row 


Boundary conditions 
imposed through 
finite volume “phantom 
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Axial Velocity Contours at Rotor 
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Axial Velocity Contours at Stator 
Inlet Plane 




Phantom cell approach 
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. ... circumferential flow 
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Figure 6.2: Illustration of ADPAC mixing plane boundary formulation for steady 
prediction of multistage turbomachinery flows. 



Figure 6.3: Illustration of mixing plane analysis (predicted Mach contours and mesh 
system) for a 3-1/2 stage compressor. 
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startup procedure was employed during this process such that this iterative startup 
procedure employs only the coarsest mesh system available. 

Steady-state solutions were normally deemed converged when the average residual 
R was reduced by a factor of 10 -3 , or when the residual has ceased to be reduced. 
Experience has shown that pressure-driven flow quantities generally converge first 
(e.g. mass flow, lift, etc.) while viscous driven flow quantities (loss) converge after 
a larger number of iterations. It is also therefiore necessary to monitor integrated 
performance parameters such as efficiency to determine when the solution is truly 
converged. 
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Chapter 7 

PARALLEL COMPUTING 


In order to address the goals of the HPCCP program described earlier, the EEE/LPS 
simulations were developed with the application of advanced parallel processing tech- 
niques as the computational foundation. Parallel processing computations using the 
ADPAC code are performed via a coarse-grained domain decomposition, and inter- 
processor message passing via either the Application Portable Parallel Processing 
(APPL) [112], Parallel Virtual Machine (PVM) [113], or Message Passing Interface 
(MPI) [114] message passing libraries. The LPS analysis was optimized for and per- 
formed on workstation cluster computing platforms using these parallel processing 
techniques. The NASA Ames Research Center davinci and babbage workstation clus- 
ters, and the NASA Lewis Research Center LACE workstation cluster were utilized 
for the analysis. The Allison Engine Company Silicon Graphics 16-processor Power 
Challenge XL server was also employed during this program. Details of each of these 
computing platforms are given in the sections which follow. 


7.1 davinci Workstation Cluster 

The davinci cluster consists of one front-end system and eight compute nodes. The 
front-end system (named davinci ) is the host that users log into. The front-end is 
a Silicon Graphics Power Challenge L with four 75 MHZ R8000 CPUs and 384 MB 
memory, and serves as the system console, compile server, file server, user home 
server, PBS server, etc. There are eight compute nodes (four two-cpu nodes, and four 
eight-cpu shared-memory nodes) with the following configuration: 


Machine Cpu Memory Swap /tmp Use 


davinci 

4 

384MB 

1.2GB 

davinci-01 

2 

128MB 

1.2GB 

davinci-02 

2 

128MB 

1.2GB 

davinci-03 

2 

128MB 

1.2GB 

davinci-04 

2 

128MB 

1.2GB 

davinci-05 

8 

2GB 

6GB 


1.2 

GB 

user home, fileserver 

0.9 

GB 

compute node, console 

0.9 

GB 

compute node 

0.9 

GB 

compute node 

0.9 

GB 

compute node 

4 

GB 

compute node 
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Figure 7.1: NASA Ames Research Center Silicon Graphics (davinci) workstation 
cluster schematic diagram (configuration circa 1995). 


davinci-06 

8 

2GB 

6GB 

4 

GB 

compute node 

davinci-07 

8 

4GB 

6GB 

4 

GB 

compute node 

davinci-08 

8 

4GB 

6GB 

4 

GB 

compute node, console 


All the machines were connected via Ethernet, FDDI, and HiPPI. ATM network 
adapters from both SGI and Fore Technology were also tested on this cluster. The 
eight compute nodes and the front-end were running the IRIX 6.2 operating system. 
PBS 1.1.8 was the job queuing system, and MPI 2.0 from SGI was the primary inter- 
processor communication library. A schematic illustration of the NASA-Ames Silicon 
Graphics ( davinci ) workstation cluster is given in Figure 7.1. 


7.2 babbage Workstation Cluster 

The NAS SP2 babbage workstation cluster is a 160-node MIMD parallel computer com- 
posed of IBM RS6000/590 workstations. On paper (and according to some bench- 
marks), the SP2 is capable of outperforming a 16-processor Cray C90. The NAS 
SP2 resulted from the HPCCPT-1 Cooperative Research Agreement (CRA) between 
NASA and a consortium led by IBM. 

Each node has at least 128 Mbytes of main memory and 2 Gbytes of disk space. 
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Figure 7.2: NASA Ames Research Center IBM RS-6000 ( babbage ) workstastion cluster 
schematic diagram (configuration circa 1995). 


Some nodes have additional memory and disk space, as well as HiPPI or FDDI. The 
SP2 also has an external filesystem accessible by all nodes. The full 160-node SP2 
has: 


23.9 Gbytes of main memory 
485 Gbytes of disk space 
342 Gbytes/second main memory bandwidth 
42.8 Gflops peak performance 


An illustration of the NASA- Ames IBM SP2 ( babbage ) workstation cluster is given 
in Figure 7.2. 


7.3 LACE Workstation Cluster 

The NASA Lewis Research Center LACE cluster is a group of thirty-two networked 
IBM RS/6000 machines ( Iace01-lace32 ) plus one “control” node called lace. Job 
submission and queuing was moderated by the LSF (Load Sharing Facility) software. 
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Figure i .3. NASA Lewis Research Center IBM RS-6000 ( LACE) workstastion cluster. 

An image of the NASA-Lewis Research Center LACE IBM RS-6000 workstation 
cluster is given in Figure 7.3. 

7.4 Allison Silicon Graphics Power Challenge XL 

The Allison Engine Company Power Challenge XL workstation consists of a 16- 
processor shared-memory parallel computing platform. The machine consists of 16 
R10000 CPU’s with 2 gigabytes of main memory completely shared across all pro- 
cessors. The operating system during this study was IRIX 6.2, with job submission 
managed by the PBS software package Version 1.1.9b. A typical machine of this type 
is illustrated in Figure 7.4. 
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Figure 7.4: Allison Engine Company Silicon Graphics Power Challenge XL parallel 
computer. 
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Chapter 8 

NEPP CODE DESCRIPTION 


One facet of the analyses performed during this study was the desire to investigate 
aspects of the “zooming” feature of the planned NPSS engine performance analysis 
architecture. In this regard, the intention was to numerically couple detailed CFD 
simulations of the EEE LP spool with an engine cycle analysis of the EEE HP core. 
This coupled simulation, would, in fact, be a complete simulation of the two-spool 
EEE engine with varying levels of fidelity for the LP and HP subsystems. The ADPAC 
analysis described in the previous chapters was directed at the 3-D CFD portion of 
this simulation strategy, while the NEPP 1-D cycle analysis was directed at the HP 
spool simulation strategy. A brief description of the NEPP analysis code is given 
below. 

The NEPP computer program [116]-[141] performs one-dimensional, steady-state 
thermodynamic performance analysis of aircraft gas turbine or jet engine configura- 
tions. Data inputs specify a standard set of components and their interconnections, 
allowing simulation of virtually any engine configuration. As many as six modes of 
engine operation may be configured to analyze multimode or variable cycle engines 
whose flowpaths and operating components vary over portions of the aircraft flight 
regimes. Physical components which may be used include propellers, inlets, ducts, 
combustors, fans, compressors, turbines, shafts, heat exchangers, flow splitters, sub- 
sonic mixers and/or supersonic ejectors, nozzles and water injectors or gas generators. 
Two options are available for gas thermodynamic properties. The default option uses 
built-in curve fits for mixtures of air and JP4, the standard hydrocarbon jet fuel. 
Alternatively, the Chemical Equilibrium Compositions (CEC) auxiliary program can 
model nearly any propellant combination or evaluate the effects of chemical dissoci- 
ation of gases. The CEC option also permits simulation of rocket components and 
fuels in an engine configuration. 

Although NEPP was originally intended to perform only thermodynamic anal- 
ysis, additional capabilities have been implemented. Simplified aircraft installation 
effects give preliminary estimates of inlet and nozzle drag forces. A turbine cooling 
algorithm estimates the gas bleed flow required for high temperature engine oper- 
ation. An approximation algorithm computes emissions of nitrogen oxides. Engine 
operation line data is accumulated for subsequent plotting on compressor and turbine 
performance maps (this feature presently requires software for plotting which varies 
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Figure 8.1: Screen illustration of NEPP analysis with NPAS user interface 


from installation to installation). 

Two additional auxiliary programs further extend the capabilities of the NEPP 
system. \fyATE estimates engine component weights. INSTAL gives more accurate 
estimates of inlet and nozzle drag forces and inlet weights, provided design details are 
available. 

There are several steps for putting together a NEPP input file to analyze an engine 
system. 

• Select the engine cycle. 

• Convert the cycle into a block diagram for NEPP. 

• Define the compressor and turbine performance maps. Exact maps for the 
application are not required, the program can scale maps as required. 

There exists a graphical user interface front-end for the NEPP code referred to as 
NPAS which simplifies the use of the code. An illustration of the NPAS user interface 
scheme for the analysis of the complete EEE engine is given on Figure 8.1. This 
complete engine model formed the basis for the reduced models described later in 
this report. 
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Chapter 9 

Component Performance 
Validation 


Component performance validation was considered a necessary milestone both in 
validating the accuracy of the analysis as well as verifying the accuracy of the ge- 
ometry specifications in the EEE Master Engine Geometry Database. During this 
phase of the program, specific subcomponent geometries were selected and analyzed 
in isolation from the other major subcomponents of the overall EEE LPS analysis. 
The engine subcomponents analyzed during the component performance validation 
phase of the program were the fan/bypass/booster compressor, the LP turbine, and 
the lobed exhaust mixer. Design point validations were also performed for the core 
compressor and the HP turbine for completeness. 

The component validation phase also served two additional purposes: the resulting 
simulations could be used as the initial conditions for the coupled EEE LPS analysis, 
and the results could also be used to evaluate, at least to some extent, the steady 
aerodynamic interaction effects resulting from subcomponent coupling in the EEE 
LPS analysis. Results from the component validation studies are summarized in the 
following sections. 


9.1 EEE Fan Section Analysis 

9.1.1 Description of Design 

The EEE fan section design is based on a unique split flow configuration selected to 
minimize mission fuel burn and direct operating cost. An illustration of the EEE 
fan section flowpath and blade arrangement is given in Figure 9.1. The EEE fan 
section design employs a full span fan rotor with a design corrected tip speed of 
1350 ft/s. and an inlet radius ratio of 0.342. The fan employs a part span shroud 
to improve structural rigidity. The fan rotor exit flow is split radially by an island 
splitter. The inner annulus of this island splitter is designed to capture 22% of the fan 
flow and employs a 1/4-height booster stage. The 1 /4-height booster stage further 
supercharges the flow entering the core and enhances core protection from foreign 
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Figure 9.1: Axisymmetric projection of EEE fan+l/4-height booster stage configu- 
ration ilustrating test data instrumentation plane locations. 


object damage. The use of the booster stage also permits a lower fan rotational 
speed, increased fan efficiency, lower fan hub aerodynamic loading, and provides for 
an easier engine growth path. The flow through the booster stage is subsequently split 
by the core inlet, with 68% of the booster flow entering the core and the remaining 
42% of the booster flow reentering the bypass flowpath through the bypass vane. 
The outer annulus flow carries the remaining 78% of the fan rotor flow through the 
bypass vanes. A detailed listing of the EEE fan section design parameters is given in 
Table 9.1. 


9.1.2 Mesh System 

Three-dimensional Navier-Stokes computations were performed for the EEE fan sec- 
tion configuration using the A D PAC analysis code. The mesh generation procedure 
described in Chapter 4 was employed for this task. The analysis included the full 
height fan with part span shroud, 1/4-height booster stage, core inlet guide vane, and 
bypass vane as shown in Figure 9.1. 

The mesh generation procedure previously described was employed to define a 
1,605,000 cell mesh distributed among 8 mesh blocks. The mesh system for the fan 
section is somewhat more complicated than the other subcomponents (such as the 
LP turbine, for example) in that the mesh block structure is not a simple end-to- 
end stack of blade rows. The various radial divisions of the flow (by the part span 
shroud, island splitter and core splitter) all require block modeling, and must still 
be compatible with the mixing plane formulation and the H-type mesh structure 
developed for the individual blade rows. This complex mesh block structure is all 
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EEE Fan Section Aerodynamic Design Parameters 

Engine Design Point: Max Climb (M=0.8, 35,000 ft ISA) 


Inlet Radius Ratio 

0.342 

Fan Specific Flow 

42.8 lbm/ft**2/s 

Corrected Tip Speed 

1350 ft./s. 

Bypass Pressure Ratio 

1.65 

Core Pressure Ratio 

1.67 

Booster Mass Flow/ 
Fan Mass Flow 

0.22 

Core Mass Flow/ 
Booster Mass Flow 

0.58 

Bypass Ratio 

6.8 


Table 9.1: EEE fan section aerodynamic design parameters. 
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EEE Fan Section Component Performance 
Validation Mesh Block Size Tabulation 


Block 

1 Index 

J Index 

K Index 

# Pts 

1 

113 

81 

49 

448497 

2 

113 

21 

49 

116277 

3 

97 

65 

33 

208065 

4 

97 

65 

33 

208065 

5 

105 

33 

33 

114345 

6 

29 

33 

33 

31581 

7 

89 

33 

49 

143913 

8 

137 

37 

49 

248381 


Table 9.2: Tabulated mesh block sizes for EEE fan section component performance 
validation analysis. 

rather easily managed by the GRIDGEN mesh generation program, but does require 
some additional thought and planning by the user. The resulting mesh block sizes 
and general relationship to the fan section components is given in Table 9.2. An 
illustration of the axisymmetric projection of the mesh system is given in Figure 9.2. 


9.1.3 Design Point Analysis 

A design point analysis was performed for the EEE fan section using the mesh system 
described in the previous subsection. The EEE fan section design bypass ratio is 6.8, 
and the fan design point represents the engine maximum climb operating point. The 
analysis was performed on a 4-processor Silicon Graphics Power Challenge L multi- 
processor computer with 1 GB of main memory. A converged solution was obtained in 
a total of 6 hours (wall clock time) using all four processors. Figure 9.3 illustrates the 
predicted fan surface static pressure contours from the analysis. Numerical predic- 
tions for the EEE fan section were compared with experimental data derived from full 
scale rig tests of the fan section [74]. Figure 9.4 illustrates a comparison of predicted 
and experimental bypass vane exit and 1 /4-stage vane leading edge spanwise total 
pressure distributions. The total pressure distributions are plotted and correlated 
with the colors of the data survey stations indicated on Figure 9.1. The character of 
the spanwise pressure distribution was very accurately captured, and was well within 
the range of test data. 
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Figure 9.3: Predicted surface static pressure contours for FEE fan plus 1/4-height 
booster stage configuration. 
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A DP AC EEE Fan + Quarter Stage Flow Prediction 

Quarter Stage Stator and Bypass Vane Exit Spanwise Total Pressure Comparison 
1.0 


0.8 


* 

o 

bm 


0.6 


! 0.4 f- 


0.2 


0.0 


O Bypass Vane Exit Arc Rakes 
A Bypass Vane Exit Radial Rakes 
O Bypass Vane - Vane Mounted 
□ 1/4 Stage Stator Vane Mounted 

Prediction - 1/4 Stage Stator L.E. 

Prediction - Bypass Vane Exit 

Prediction - Bypass Vane Exit 


Part-span shroud 


Island Splitter 



0.8 1.0 1.2 1.4 1.6 

Total Pressure Ratio (Pt/Pt* inlet) 


1.8 


Figure 9.4: Comparison of predicted and experimental spanwise total pressure dis- 
tributions at bypass vane exit and 1/4-stage vane leading edge for EEE fan plus 
1/4-height booster stage configuration. 

9.1.4 Off-Design Analysis 

In order to investigate the off-design analysis capabilities of the EEE fan section 
model, a number of predictions were performed at 100% corrected speed with varia- 
tions in both fan exit static pressure and fan section bypass ratio. These off-design 
results were obtained by prescribing the flow entering the core, and adjusting the 
bypass exit static pressure until the desired fan inlet flow was achieved. Excursions 
in predicted bypass ratio ranged from 6.0 to 10.8. 

Predictions of overall performance were compared with measured data derived 
from full-scale rig tests of the fan section [74]. A comparison of predicted and ex- 
perimental overall pressure ratio and adiabatic efficiency versus corrected mass flow 
rate for the core stream flow (downstream of the core inlet guide vane) of the EEE 
fan section is given in Figure 9.5. The corresponding maps for the fan bypass stream 
flow (downstream of the fan bypass vane) is given in Figure 9.6. The data on these 
figures illustrates the overall capabilities of the EEE fan design. Bold symbols on each 
figure illustrate the test performance at extreme high and low values of bypass ratio. 
It is interesting to note that in both the test and the prediction, bypass ratio did not 
significantly alter the characteristics of the bypass stream, but does have a significant 
effect on the core stream flow. The overall character of the off-design performance 
predictions displayed good agreement with the test data. 
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Core Stream Total Pressure Ratio Adiabatic Efficiency 




Figure 9.5: Comparison of predicted and experimental total pressure ratio and adia- 
batic efficiency versus corrected flow rate for the core inlet of the EEE fan section. 
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Figure 9.6: Comparison of predicted and experimental 
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section. 
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9.2 EEE Low Pressure (LP) Turbine Analysis 

9.2.1 Description of Design 


The EEE LP turbine consists of a 5-stage design employing moderately loaded airfoils 
and a rather high (25 degrees) endwall slope. The 5-stage design was based in part 
on results obtained from studies of highly loaded fan turbine technology development 
at General Electric, and from system studies aimed at minimizing direct operating 
cost (DOC). The EEE engine LP turbine design is coupled to the HP turbine via 
a short (3 in.) transition duct. The relatively high bypass ratio (6.8) of the EEE 
fan section, and subsequent reduced core flow requires high specific energy from the 
fan-drive (LP) turbine. The design efficiency goals for the LP turbine were 91.1% 
for the integrated core/low spool (ICLS) test and 91.7% for the flight propulsion 
system (FPS) at the engine design point (M=0.8, 35,000 ft. altitude ISA). The LPT 
maximum tip diameter was set by mechanical and configuration control requirements 
at 46.5 in. The outer wall slope was also limited to 25 degrees (established as a 
maximum to maintain good aerodynamic performance) through stage 3, transitioning 
to a cylindrical outer wall at the stage 5 exit. A table of pertinent LP turbine design 
and operating parameters is given in Table 9.3. 


9.2.2 Mesh System 

Mesh generation was based on the 4-step procedure described in Chapter 5. A mesh 
system consisting of 10 mesh blocks (1 per blade row for 5 stages) containing 1,660,000 
computational cells was assembled. A table of the mesh block sizes for the blade 
passage meshes is given in Table 9.4. An illustration of the axisymmetric projection 
of the LP turbine mesh system is given in Figure 9.7. 


EEE LP Turbine Aerodynamic Design Parameters 

Engine Design Point: Max Climb (M=0.8, 35,000 ft ISA) 



Stage 1 

Stage 2 

Stage 3 

Stage 4 

Stage 5 

Energy Extraction 
Ah J/g (BTU/lbm) 

73.04 (31.4) 

79.09 (34.0) 

82.11 (35.3) 

70.25 (30.2) 

49.31 (21.2) 

Pressure Ratio (Pt/Pt) 

1.30 

1.35 

1.40 

1.36 

1.26 

Aero Loading 
Ah/2u 2 

1.71 

1.58 

1.43 

1.13 

0.80 

Flow Coefficiect 

1.25 

1.08 

1.04 

0.98 

1.07 


Vz/u 


Table 9.3: EEE LP turbine aerodynamic design parameters. 
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EEE LP Turbine Component Performance 
Validation Mesh Block Size Tabulation 


Block 

1 Index 

J Index 

K Index 

# Pts 

1 

97 

49 

33 

156849 

2 

97 

49 

33 

156849 

3 

97 

49 

33 

156849 

4 

97 

49 

33 

156849 

5 

97 

49 

33 

156849 

6 

97 

49 

33 

156849 

7 

97 

49 

33 

156849 

8 

97 

49 

33 

156849 

9 

97 

49 

33 

156849 

10 

97 

49 

33 

156849 


Table 9.4: Tabulated mesh block sizes for EEE LP turbine component performance 
validation analysis. 



Figure 9.7: Axisymmetric projection of EEE LP turbine component validation mesh 
system. 
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Figure 9.8: Predicted surface static pressure contours for EEE LP turbine. 


9.2.3 Design Point Analysis 

Design point numerical simulations of the EEE Low Pressure (LP) turbine were per- 
formed to permit comparison with 2/3 scale rig test data [101]. The analysis was 
performed on a Silicon Graphics Power Challenge L multiprocessor computer with 1 
GB of main memory. Converged solutions were obtained in a total of 3 hours (wall 
clock time) using four processors. Xote that the turbine simulation was nearly twice 
as fast as the fan section simulation in spite of the fact that approximately 20% more 
mesh points were involved. This feature results from the generally favorable pressure 
gradients involved in the turbine flow, leading to a rapid definition of the boundary 
layer flow. Conversely, the fan section flow involves predominantly adverse pressure 
gradients requiring significantly more computation time to resolve. The rapid compu- 
tation time for the turbine clearly indicates the suitability of the analysis for design 
cycle studies. In fact, more time was involved in generating suitable meshes than was 
involved in the aerodynamic analysis itself. Predicted turbine surface static pressure 
contours are illustrated in Figure 9.8. This figure illustrates the three-dimensional 
nature of the blading and the general arrangement of the LP turbine. 

A comparison of predicted and experimental spanwise variation of fifth stage exit 
total pressure and total temperature profiles is given in Figure 9.9. This preliminary 
analysis was based on a simple flat inlet profile of total pressure and total tempera- 
ture and employed the exact blade and endwall definitions provided in the original 
Master Engine Geometry Database. The correlation between rig test and calculation 
is excellent in the 20%-80% radial span region. Noticeable discrepancies exist in the 
near endwall regions. These discrepancies were assumed to be due to the fact that no 
clearance flows, turbine hub overlap geometry, or shrouded rotor cavity geometries 
were modeled in this initial prediction. 

In order to resolve differences between prediction and experiment near the end- 
walls, several additional calculations were performed to assess the effects of variations 
in geometry, flow parameters, etc. The variations tested included modifications to the 
first stage vane setting angle, modifications to the inlet flow profile, and the addition 
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ADPAC EEE LP Turbine Rig Test Analysis 


5th Stage Exit Span wise Total Pressure Variation 



ADPAC EEE LP Turbine Rig Test Analysis 

5th Stage Exit Spanwlsc Total Temperature Variation 



Figure 9.9: Comparison of predicted and experimental spanwise total pressure and 
total temperature distributions for fifth stage exit of the EEE LP turbine 2/3-scale 
test rig. 
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of a shrouded rotor endwall cavity model. Each of these variations are described in 
detail in the sections below. For each case, the effects are compared based on the pre- 
dictions from the original, smooth endwall, flat inlet profile, unmodified LP turbine 
geometry as it existed in the EEE Master Engine Geometry Database at the begin- 
ning of this study. Each calculation was performed using a common static pressure 
ratio specification at the turbine exit hub surface. 

9.2.4 Effect of Variations in First Vane Setting Angle 

GE engineers familiar with the actual test rig and EEE engine geometry recommended 
a 1 degree (open) reset of the LP turbine first stage vane. The effect of the reset on 
the LP turbine exit spanwise flow profiles is illustrated on Figure 9.10. A distinct 
improvement in the predicted total temperature distribution was observed at the 
turbine exit, particularly near the tip, for the calculation involving the modified 
geometry. Given this observation, all further calculations were based on the modified 
first stage vane orientation. 

9.2.5 Effect of Variations in Inlet Profile 

Several multistage calculations were performed with variations in the first vane input 
spanwise total pressure, total temperature, and flow angle profile distributions. The 
profiles are categorized as flat (baseline, essentially no variation across the span except 
at the tip), boundary layer (BL - 10% thick total pressure deficit at the endwalls), 
and engine (derived from a simulation of the HP turbine exit flow). An illustration of 
the spanwise variation of inflow total pressure and total temperature from the three 
profiles is given in Figure 9.11. 

Figure 9.12 illustrates the comparison of predicted and experimental LP turbine 
exit spanwise total pressure and total temperature profiles for each of the inlet pro- 
file variations described above. Note that there is not a significant change in the 
exit profile total pressure characteristics with variations in inlet profile specification. 
This is partially due to the fact that each calculation is run to the same exit static 
pressure ratio. There is some variation in the exit total temperature distributions, 
although this behavior essentially correlates with the inlet total temperature profile 
characteristics. 

9.2.6 Effect of Variations in Endwall Geometry 

The final comparison of results involved discrete modeling of the turbine shrouded 
rotor seal cavities. The calculations described above were all performed using a 
geometry model based on a smooth, continuous endwall definition. In reality, the 
endwalls are quite discontinuous and irregular due to the use of shrouded rotors 
and overlapping geometry, (see e.g. Figure 9.13) and these irregularities can have a 
significant impact on the primary gas path flow. Previous experience in predicting 
flows through compressor seal cavities suggests that the seal cavities themselves can 
often be modeled using two-dimensional techniques, and then subsequently coupled 
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ADPAC EEE LP Turbine Analysis 

5th Stage Exit Spanwise Total Pressure Variation 



ADPAC EEE LP Turbine Analysis 


5th Stage Exit Spanwise Total Temperature Variation 



Figure 9.10: Comparison of predicted and experimental spanwise variation in fifth 
stage exit total temperature distributions for EEE LP turbine analyses with variations 
in first vane reset and endwall modeling. 
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Total Pressure Ratio 


Total Temperature Ratio 


Figure 9.11: Comparison of spanwise variation of inflow total pressure and total 
temperature ratio profiles for profiles for the EEE LP turbine analyses with variations 
in inlet profile. 
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ADPAC EEE LP Turbine Analysis 


5th Stage Exit Spanwise Total Pressure Variation 



Total Pressure Ratio (Pt,exit/Pt, inlet) 


ADPAC EEE LP Turbine Analysis 


5th Stage Exit Spanwise Total Temperature Variation 



Figure 9.12: Comparison of predicted and experimental spanwise variation in fifth 
stage exit total temperature distributions for EEE LP turbine analyses with variations 
in inlet profile. 
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Axisymmetric Projection Geometry and Mesh Block Outline 



Axisymmetric Projection of H-Type Mesh System 



Figure 9.13: Illustration of realistic LP turbine endwall irregularities and CFD mod- 
eling techniques for the EEE LP turbine component validation study. 
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Inlet 

First 

Endwall 

Mass 

Pt, 

Tt„ 

Adiabatic 

Profile 

Vane 

Type 

Flow 

Exit 

Exit 

Efficiency 

Type 

Reset 

(Ib/s) 

(psia) 

(deg. R) 

(Mass- Averaged) 

Flat 

0 

Smooth 

67.652 

10.299 

514.01 

91.72% 

BL 

0 

Smooth 

67.366 

10.292 

514.49 

91.60% 

Engine 

0 

Smooth 

68.228 

10.314 

512.70 

91.99% 

BL 

1 open 

Smooth 

67.896 

10.308 

514.57 

91.64% 

Flat 

0 

Cavity 

67.146 

10.304 

526.86 

86.74% 

BL 

0 

Cavity 

66.705 

10.284 

522.92 

88.27% 

BL 

Notes: 

1 open 

Cavity 

67.784 

10.316 

522.60 

88.53% 


1 . Nominal inlet total pressure « 45.0 psia 

2. Nominal inlet total temperature « 750 deg. R 

3 Approximate variation in computed mass flow from blade row to blade row: 
Smooth Endwall: 0.3% 

Cavity Endwall: 2.0% 


Table 9.5: Comparison of predicted overall performance parameters due to variations 
in inlet profile, endwll model, and first vane reset for the EEE LP tuibine 2/3 scale 
test rig. 

with the 3-D blade passage flow through averaging techniques similar to a mixing 
plane. This was the approach adopted in this study to minimize the computational 
effort involved with modeling this more complicated flow case. 

Figure 9.14 illustrates a meridional projection of the LP turbine mesh employing 
smooth endwalls (upper plot) and the same configuration where the shrouded rotor 
seal cavities are discretely modeled (lower plot). An illustration of the piedicted 
axisymmetric-averaged Mach number contours for the EEE LP turbine with shrouded 
rotor cavity model is presented in Figure 9.15. The influence of the cavities would 
appear to be limited to local regions along the case near the inflow/outflow openings 
of the cavity. The resulting effect on the predicted spanwise profiles at the exit of 
the turbine are also illustrated on the plots on Figure 9.16. The characteristics of the 
spanwise profiles were not significantly altered due to the addition of the shrouded 
rotor cavity model; however, significant changes in the overall turbine peifoimance 
parameters were detected. These changes are discussed in more detail in the following 
separate subsections. 

9.2.7 Summary of Variations in Turbine Parameters on De- 
sign Point Performance 

A summarv of the overall performance characteristics due to the variations described 
above is given in Table 9.5. In terms of overall performance, variations in inlet 
profile did not appear to have a significant effect on the predicted mass flow rate, 
exit total pressure, total temperature, or efficiency for the smooth endwall model. 
In the cavity endwall model calculations, the differences due to inlet profile were 
more pronounced. Variations in first vane reset primarily affected the predicted mass 
flow rate. The 1 degree (open) reset of the first stage vane resulted in an increase 
in flow of 0.78% for the smooth endwall test case, and an increase of 1.59% for 
the cavity endwall model test case. Finally, in terms of the effects of variations in 
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Figure 9.14: Comparison of axisymmetric projection of mesh systems for the EEE LP 
turbine with smooth endwalls (upper) and with modeled shrouded rotor seal cavities 
(lower). 
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Figure 9.15: Illustration of predicted axisymmetric-averaged Mach number contours 
for the EEE LP turbine with shrouded rotor cavity endwall model. 


the endwall model, the most prominent characteristics were reductions in predicted 
mass flow rate and adiabatic efficiency due to the cavity endwall flow model. The 
reduction in efficiency was quite dramatic - on the order of 3%-5% depending on 
the test case. One problem encountered during this evaluation was an inability to 
consistently maintain a constant mass flow from blade row to blade row in the cavity 
endwall solutions. Typical variations in mass flow from blade row to blade low in the 
multistage simulations using the smooth endwall model was 0.3%, while the cavity 
endwall model resulted in blade row to blade row variations as high as 2.0%. The 
large variation in the cavity flow model was a result of the complicated mixing-plane 
arrangement employed to numerically couple the 2-D cavity passage openings with 
both the upstream and downstream neighboring blade row 3-D mesh systems. Given 
this large level of mass flow variation, the large predicted efficiency reduction due to 
the addition of the shrouded rotor cavities should be interpreted qualitatively, not 
necessarily quantitatively. 


9.2.8 Off-Design Analysis 

A summary of the off-design component performance validation efforts for the EEE 
LP turbine are presented in this section. ADPAC solutions for the LP turbine were 
compared with GE scaled test rig Block II, Configuration 5 experimental data [101]. 
ADPAC was employed to generate several operating point solutions near the design 
blade-jet speed ratio ( ujC a = 0.412 where u is the turbine inlet mean axial velocity 
and C 0 is the turbine tip speed) for the 2.4 million point LP mesh. The mesh included 
2-D shrouded rotor seal geometries. A constant blade-jet speed ratio was set by fixing 
the inlet-to-exit pressure ratio and solving for the necessary shaft rotational speed. 
Pressure ratios of 2.0, 3.0 and 4.76 were used for computations and ADPAC data was 
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ADPAC EEE LP Turbine Analysis 


5th Stage Exit Spanwise Total Pressure Variation 



ADPAC EEE LP Turbine Analysis 



Figure 9.16: Comparison of predicted and experimental spanwise variation in fifth 
stage exit total temperature distributions for EEE LP turbine analyses with variations 
in inlet profile. 
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Corrected Speed (N/(T t39 ) 

Figure 9.17: Comparison of predicted (ADPAC) and measured equivalent energy 
extraction for the Energy Efficient Engine (EEE) LP turbine. 


reduced to enable comparison of equivalent energy extraction, inlet flow function, 
total-to-total efficiency and total-to-static efficiency. The comparisons are displayed 
in Figures 9.17-9.20. 

The predicted trends for equivalent energy extraction and inlet flow function com- 
pare well with the scaled rig test data. The absolute levels of these performance 
parameters is also predicted reasonably well, in spite of the numerous uncertainties 
concerning the test vehicle and the data reduction procedures. The predicted trends 
in efficiency were also captured reasonably well; however, the predicted efficiencies 
are consistently 2%-4% low. This difference was due, in part, to the modeling of 
shrouded rotor seal flow, which caused a 3%-5% drop in adiabatic efficiency when 
compared to the smooth endwall prediction. The discrepancy in efficiency varied 
considerably based on the numerical method used to compute the efficiency (total 
temperature, angular momentum change, mass averaging versus area averaging, etc.). 
The large number of unpublished features of the test rig operation, and the uncer- 
tainties associated with the numerical cavity model prohibited timely investigation of 
this discrepency. 
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Total-to-Total Efficiency 



Figure 9.19: Comparison of predicted (ADPAC) and measured total to total adiabatic 
efficiency for the Energy Efficient Engine (EEE) LP turbine. 
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Total-to-Static Efficiency 



Pressure Ratio (P T42 /P S55 ) 


Figure 9.20: Comparison of predicted (ADPAC) and measured total to static adia 
batic efficiency for the Energy Efficient Engine (EEE) LP turbine. 
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9.3 EEE Core Compressor Analysis 

9.3.1 Description of the Design 

The core (HP) compressor system for the EEE is a 10-stage axial flow compressor 
with a design pressure ratio of 22.6:1 [83]. The design corrected tip speed is 456 
m/s (1495 ft/s), resulting in a corrected airflow of 53.5 kg/s (118 lb/s). The goal 
adiabatic and polytropic efficiencies for the EEE core compressor were 86.1% and 
90.6%, respectively. The design was the result of an extensive compressor optimiza- 
tion study to identify desirable compressor design features for a subsonic transport 
engine. This optimization included analyses of the effects of aspect ratio, solidity, 
inlet specific flow, exit Mach number, reaction ratio, inlet radius ratio, exit radius 
ratio, and number of stages. The effects of each parameter were examined based on 
efficiency, weight, cost, aircraft direct operating cost (DOC) and fuel consumption. 
Two engine configurations were considered during this early study: an engine having 
a core compressor total pressure ratio of 14 with booster stages on the LP spool, and 
an unboosted core compressor with a total pressure ratio of 23. It was determined 
that best compressor performance was achieved using medium values of aspect ratio, 
solidity, and reaction ratio, and low values of inlet radius ratio, inlet specific flow, and 
exit Mach number. The 10-stage configuration offered the best overall combination 
of desirable features: compactness, low cost, high efficiency, low DOC and low fuel 
usage. Design parameters for the EEE core compressor are tabulated in Table 9.6. 


9.3.2 Mesh System 

Simulation of the core compressor was intended primarily as a check on the mesh 
generation system developed for the EEE IGES component definitions, and the fi- 
delity of the EEE core compressor geometry database. As such, only a design point 
simulation was performed on the baseline core compressor geometry (variable stator 
schedules in their “design” setting). The resulting mesh system was composed of 21 
mesh blocks (1 per blade row for IGV and 10 stages) and is illustrated in an axisym- 
metric projection in Figure 9.21. The total number of computational cells in this 
mesh is 3,553,000. A typical block size is 97x33x49 (axial, radial, tangential). This 
mesh density is typical of design analysis calculations. 


9.3.3 Design Point Analysis 

A design point analysis was performed for the EEE core compressor. The simulation 
was performed primarily as a check of the geometry database and the mesh generation 
and solution procedures. The solution was performed on a Silicon Graphics Power 
Challenge XL computer, employing 12 of the 16 available processors. This resulted 
in no more than two blade rows per processor. Subsequent calculations employed 21 
processors with one blade passage per blade row. The resulting decrease in CPU time 
was nearly linear with the number of processors. 
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EEE Core Compressor Aerodynamic Design Parameters 

Engine Design Point: Max Climb (M=0.8, 35,000 ft ISA) 


Corrected Tip Speed 

(m/s) 

(ft/s) 

456 

1495 

Inlet Radius Ratio 


0.503 

Flow/Annulus Area 

(kg/s/m 2 ) 

(lbm/s/ft 2 ) 

185.5 

38.0 

Rotor 10 Exit Hub Speed 

(m/s) 

(ft/s) 

352.7 

1157 

Rotor 10 Exit Radius Ratio 


0.93 

Outlet Guide Vane Exit Mach Number 

0.30 

Number of Rotors and Stators 

1672 

Average Aspect Ratio 


1.48 

Average Pitch Solidity 


1.36 

Adiabatic Efficiency 


85.7% 

Stall Margin Potential 


25% 


Table 9.6: EEE core compressor aerodynamic design parameters. 



Figure 9.21: Axisymmetric projection of EEE core compressor component validation 
mesh system. 
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Including solution initialization, which was performed on coarser meshes using 
the ADPAC full multigrid initialization routine (essentially a combined grid sequenc- 
ing/multigrid solution strategy), the complete design point simulation using the 12- 
processor configuration was obtained in 15 hours. This clearly indicates that a com- 
plete operating map could be derived in essentially one day given enough available pro- 
cessors. For example, using two 21-processor machines of current computing power, 
a single constant speed operating line (6 different pressure ratios at a constant speed) 
could be evaluated for this compressor in about 14 hours using the ADPAC solution 
strategy. This estimate includes the reduction in total solution time afforded by the 

ability to restart from previous solutions. 

Mass-averaged estimates of performance from the HP compressor design point 
simulation were used to predict a mass flow rate, total pressure iatio, and adiabatic 
efficiency of 120.8 lbm/s, 22.37, and 86.6%, respectively. These estimates compare 
very well with the corresponding design values listed above. The prediction also 
demonstrated good agreement based on the measured data [83] for this compiessor, 
although there were obvious potential sources of error such as bleed flows, clearance 
changes, etc. 

A comparison of predicted and experimental spanwise total pressure distributions 
aft of the sixth stage rotor and aft of the tenth stage stator aie piesented on Fig- 
ures 9.22 and 9.23, respectively. Test data from both the compressor rig test and the 
engine core test are included on both figures. Since the prediction and expeiiment 
represent operation at slightly different total pressure ratios, the absolute levels of 
total pressure are slightly mismatched, but it is clear from Figure 9.22 that the span- 
wise character of the flow is very accurately predicted. There is a large difference 
between the rig and core engine test data for the tenth stage stator exit data plotted 
on Figure 9.23. This large difference was attributed to a rather large difference in 
rotor tip clearance which degraded the performance of the outei endwall flow foi the 
compressor rig test. 
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EEE HP Compressor (IGV+10 Stages) 


Rotor 6 Exit Spanwise Total Pressure Survey 



Total Pressure Ratio 


Figure 9.22: Comparison of predicted and experimental spanwise total pressure ratio 
distribution aft of the sixth stage rotor for the EEE HP compressor (design point operation). 
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EEE HP Compressor (IGV+10 Stages) 


Diffuser Exit Spanwise Total Pressure Survey 



Figure 9.23: Comparison of predicted and experimental spanwise total pressure ratio 
distribution aft of the tenth stage stator for the EEE HP compressor (design point operation). 


NAS A/CR— 1998-206597 


66 




Parameter (Units) 


Max 

Climb 


Inlet Temperature 

(deg. K) 
(deg. R) 

1588 

2858 

Energy 
Ah n 

(J/kg/de. K) 
(But/lbm/deg. R) 

353.4 

0.0844 

Speed 

ufyr 

(rad/sec//jbeg. K) 
(rpm^deg. R) 

33.19 

236.2 

Corrected Flow 
W/y^T/P 

((yfiep. K/sec/Pa) 
(Ibrry^deg. R/sec/psi) 

0.8648 

17.65 

Loading 

Ah/2U 2 


0.635 

Efficiency 

(n.%) 


91.9 


Max 

Max 

Max 

Sea Level 

Climb+4% 

Climb 

Cruise 

Takeoff+27 F 

1591 

1557 

1515 

1618 

2863 

2802 

2728 

2913 

353.4 

355.5 

353.4 

354.6 

0.0844 

0.0849 

0.0844 

0.0847 

33.78 

33.56 

33.68 

34.22 

240.4 

238.9 

239.7 

243.6 

0.8913 

0.8643 

0.8638 

0.8628 

18.19 

17.64 

17.63 

17.61 

0.625 

0.624 

0.616 

0.599 

91.9 

92.4 

92.4 

92.1 


Table 9.7: EEE HP turbine critical operating data. 


9.4 EEE HP Turbine Analysis 

9.4.1 Description of the Design 

The EEE HP turbine design [84] evolved from overall engine integration and system 
studies performed at General Electric Corporation during the development of the 
EEE engine test vehicles. The design point for the HP turbine was operation at a 
Mach number of 0.8, at 35,000 ft. ISA. The efficiency goal was 92.4%. A summary 
of the EEE HP turbine critical operating data is given in Table 9.7. 

The final HP turbine configuration was the result of detailed studies aimed at as- 
sessing the potential benefits of geometric alterations about a baseline design which 
resulted from the early cycle studies. The alterations considered were: 

• number of stages 

• outer diamater 

• annulus height 

• stage work distribution 

The EEE HP turbine consists of a 2-stage design with moderately loaded airfoils. 
A summary of stage aerodynamic parameters for the EEE HP turbine is given in 
Table 9.8. The individual airfoil blade aerodynamic geometry parameters are listed 
in Table 9.9. 


9.4.2 Mesh System 

The mesh system for the EEE HP turbine consisted of 4 mesh blocks with a total of 
627,396 points. A typical mesh block size for each blade row was 97x33x49 (axial, 
radial, tangential). An axisymmetric projection of the EEE HP turbine component 
validation mesh system is given in Figure 9.24. 
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Parameter 


Stage 


1 2 


Pressure Ratio 

2.18 

2.18 

Dh/2U 2 

0.69 

0.56 

Tip Speed (Take-off) (m/s) 

513.9 

535.2 

(ft/s) 

Cooling and Leakage (%) 

1686.0 

- 1 « D _ 

1 756 



Exit Mach Number 

0.34 

0.43 

Reaction 

0.38 

0.35 

Swirl, degrees 

15 

1 

Number of Vanes 

46 

48 

Number of Blades 

76 

70 

Radius Ratio 

0.88 

0.82 

%Tip Clearance 

1.0 

0.6 


Table 9.8: EEE HP turbine stage aerodynamic parameters. 


Parameter 

Stage 1 
Vanes 

Stage 1 
Vanes 

Stage 2 
Blades 

Stage 2 
Blades 

Number 

46 

48 

76 

70 

Solidity 

AW/t 

0.71 

1.07 

0.96 

1.06 

Zweifel Number 

0.67 

0.79 

1.08 

1.03 

% Trailing Edge 
Blockage 

7.2 

6.6 

8.1 

7.4 

Aspect Ratio 
AR=h/do 

3.3 

4.4 

3.8 

4.6 

Unguided Turn 
DB S 

8.4 

11.0 

13.0 

15.5 


Table 9.9: EEE HP turbine stage blade aerodynamic geometry. 
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Figure 9.24: Axisymmetric projection of EEE HP turbine multi-block H-type mesh 
system. 

9.4.3 Design Point Analysis 

A design point numerical simulation of the EEE High Pressure (HP) turbine was 
performed to permit comparison with test data from the full scale warm air rig test 
performed during the EEE engine development cycle [100]. Although the warm air 
test rig included cooling flow, no attempt was made in the present numerical analysis 
to account for the effects (both aerodynamic and thermal) of the cooling flow system. 
The results must therefore be interpreted with this limitation in mind. The analysis 
was performed on a Silicon Graphics Power Challenge XL multiprocessor computer 
with 2 GB of main memory. Converged solutions were obtained in a total of 3 hours 
(wall clock time) using four processors. 

Predicted turbine surface static pressure contours are illustrated in Figure 9.25. 
The orientation of the stages and nature of the blading is evident in this picture. 

A comparison of predicted and experimental spanwise variation of second stage 
exit total pressure and total temperature profiles is given in Figure 9.26. The correla- 
tion between rig test and calculation displays a consistent deviation in both temper- 
ature and pressure across the entire span. This deviation (albeit small) is believed to 
be due to the fact that cooling air injection present in the warm air turbine rig test 
was not modeled in the numerical simulation. Since the only purpose of this simula- 
tion was to validate the geometry and solution procedure, the original solution was 
deemed sufficient for this purpose, and no further effort to identify this discrepency 
was attempted. 
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Figure 9.25: Predicted surface static pressure contours for EEE HP turbine. 
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Energy Efficient Engine High Pressure Turbine 


Turbine Exit Total Pressure Ratio 



Total Pressure Ratio (Tt, cxit/TUnlet) 


Energy Efficient Engine High Pressure Turbine 


Turbine Exit Total Temperature Ratio 



Figure 9.26: Comparison of predicted and experimental spanwise distributions of sec- 
ond stage exit total pressure and total temperature profiles for the EEE HP turbine. 
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9.5 Lobed Exhaust Mixer Analysis 

9.5.1 Description of the Design 

Static scale model tests were conducted to evaluate exhaust system mixers for a 
high bypass ratio engine as part of the NASA sponsored Energy Efficient Engine pro- 
gram [36]. Gross thrust coefficients were measured for a series of mixer configurations 
which included variations in the number of mixer lobes, tailpipe length, mixer penetra- 
tion, and length. All of these parameters have a significant impact on exhaust system 
performance. In addition, flow visualization pictures and pressure/ temperature tra- 
verses were obtained for selected configurations. Parametric performance trends were 
defined based on these results. Mixer configuration variables included lobe number, 
penetration and perimeter, as well as several cutback mixer geometries. Mixing effec- 
tiveness and mixer pressure loss were determined using measured thrust and nozzle 
exit total pressure and temperature surveys. These scaled results provided a data 
base to aid the analysis and design/development of the EEE mixed-flow exhaust sys- 
tem. The final EEE Flight Propulsion System (FPS) lobed exhaust mixer employed 
a scalloped, 12-lobe design based on the results of the extensive rig testing. 


9.5.2 Mesh System 

The mesh system for the EEE lobed exhaust mixer represented one of the more 
challenging aspects of this project. Since this geometry is dissimilar to the bladed 
flowpath geometries of the fan, compressor and turbine sections, mesh generation 
was performed essentially by hand using the GRIDGEN mesh generation program. 
A partial geometry database was constructed by NASA during this study and was 
employed for the EEE LPS simulations described in this section and the following 
chapter. The geometry is at least representative of the final design, but there remains 
some uncertainty as to the complete accuracy of the lobed surfaces. In addition, the 
actual test article employed scallops on the lobes to enhance mixing. Since no detailed 
information on scallop configuration was available, the cut-outs were not modeled in 
this study. An illustration of the modeled surfaces of the EEE lobed exhaust mixer is 
given in Figure 9.27. The EEE lobed exhaust mixer mesh system along the lobe plane 
of symmetry is given in Figure 9.28. A total of 9 mesh blocks were employed to define 
the coannular engine flow streams and the external flow stream. An illustration of 
the mesh system at the mixer plane is given in Figure 9.29. 

The final mesh block sizes and total number of computational cells for the lobed 
exhaust mixer component validation study are tabulated in Table 9.10. 

9.5.3 Design Point Analysis 

A design flow analysis was performed for the EEE lobed exhaust mixer using the 
ADPAC code. Results from the analysis were integrated and qualitatively compared 
to the test data from the rig test study [36]. Only a qualitative compaiison was 
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Figure 9.27: Illustration of EEE lobed exhaust mixer geometric surfaces modeled 
during the component validation study. 


Number of blocks: 9 


Block 

ISize 

J Size 

K Size 

1 

65 

45 

81 

2 

81 

9 

81 

3 

81 

49 

81 

4 

81 

49 

81 

5 

73 

49 

81 

6 

73 

49 

81 

7 

49 

49 

61 

8 

65 

49 

81 

9 

49 

49 

61 


Total Number of Computational Cells: 2,275,992 


Table 9.10: Tabulation of EEE lobed exhaust mixer mesh block sizes and total number 
of computational cells employed during the component validation study. 
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Figure 9.28: Illustration of EEE lobed exhaust mixer symmetry plane mesh surfaces 
employed during the component validation study. 
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Figure 9.29: Illustration of EEE lobed exhaust mixer exit plane mesh system em- 
ployed during the component validation study (analysis employs one lobe and assumes 
periodicity from lobe to lobe). 
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possible due to uncertainty between the modeled mixer and the geometries described 
in the rig tests. 

Figure 9.30 illustrates the predicted surface static pressure contours from the de- 
sign point analysis for the EEE lobed exhaust mixer. Total temperature contours 
on an axial plane roughly one diameter downstream of the nozzle exit plane are also 
illustrated on this figure. The symmetric horseshoe-shaped total temperature regions 
result from the secondary flow vorticies which develop as a result of the lobed mixer. 
Figure 9.31 illustrates a series of iso-surfaces defining boundaries of constant total 
temperature for temperature ratios varying from 1.8 to 1.1. The high temperature 
surfaces are confined within the mixer as the large temperature differences between 
the two streams are initially reduced rather rapidly. At lower temperature ratios, and 
consequently farther downstream, the lobed mixer flow patterns control the shapes of 
the constant temperature surfaces. The vortical nature of the flow displays a bifurca- 
tion of the iso-surface (24 segments as opposed to 12) for the iso-surface defined by a 
total temperature ratio of 1.3. The iso-surface returns to a 12-segment configuration 
for lower temperature ratios. 

Spanwise total temperature profiles at the mixer/nozzle exit are illustrated in 
Figure 9.32. Predicted and experimental total temperature ratios are plotted against 
a normalized nozzle area distribution along several circumferentially spaced arrays 
spanning a single half-lobe of the mixer. The test data was derived from a study [36] 
of mixer configurations of varying penetration, area ratio, etc. To validate the mixer 
predictions, test data was derived from an essentially equivalent mixer (Configura- 
tion F3, 12 lobes, 39% penetration) which was tested under the referenced study. In 
general, the spanwise characteristics of the mixer are qualitatively captured, partic- 
ularly along the lobe radial peak (Station A on Figure 9.32 survey. There is some 
noticable disagreement between prediction and test at survey Stations D and E. This 
discrepencv is likely due to the fact that the numerical and test mixer geometries 
were not exactly similar, and also due to the generally accepted observation that the 
algebraic turbulence model employed in the present analyses is not well suited for 
temperature mixing problems of this sort. The algebraic model does not promote 
turbulent mixing at the shear layer between the two streams, and the general con- 
sequence is that predicted temperatures tend to display more abrupt profile changes 
than the test data. 
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Figure 9.30: Predicted surface static pressure contours and axial plane total temper- 
ature contours (one diameter aft of nozzle exit) for the EEE lobed exhaust mixer. 
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Figure 9.31: Predicted iso- temperature surfaces for EEE lobed exhaust mixer simu- 
lation illustrate temperature distribution patterns due to mixing. 
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Nozzle Area Distribution (A/A.total) Nozzle Area Distribution (A/AtotaJ) 




Station D Station E Station F 



Total Temperature Ratio Total Temperature Ratio Total Temperature Ratio 


Figure 9.32: Comparison of predicted and experimental radial total temperature 
surveys for the EEE lobed exhaust mixer. 
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Chapter 10 

EEE/LP Subsystem Analysis 


This chapter deals with the results of numerical modeling of the Low Pressure (LP) 
Subsystem of the General Electric (GE) Energy Efficient Engine (EEE). The LP 
Subsystem model was developed following the component validation studies described 
in the previous chapter. The LP Subsystem analysis employed both fixed and variable 
core boundary specifications based on the results of an engine cycle model for the High 
Pressure (HP) core behavior. 


10.1 LP Subsystem Mesh Construction 

Grid generation for the EEE LP subsystem analysis was based essentially on collect- 
ing the individual meshes for the major subcomponents (fan, HP/LP turbines and 
lobed mixer) employed during the component validation study. The existing fan, 
quarter-height booster stage, HP turbine, LP turbine, and lobed mixer subsystem 
component meshes were assembled for this purpose. In addition, new meshes were 
generated using GRID GEN to model those regions which were not discretized by any 
of the component validation models. These new regions included the forward-most 
flow in the inlet, external flow about the nacelle, and the bypass duct flow between 
the fan section bypass vane and the lobed exhaust mixer. For computational sim- 
plicity, these new regions were modeled in a two-dimensional fashion (the analysis is 
certainly not limited in this respect), and were computational coupled to the three- 
dimensional domains using the ADPAC mixing plane strategy (see e.g. Figure 10.1). 
It should be emphasized that all primary components (blade rows, for example) were 
still modeled with 3-D mesh systems. The collection and assembly of these meshes 
resulted in a numerical model of the entire EEE (minus the engine core compres- 
sor and combustor). It should be noted that although the high pressure compressor 
and combustor were not discretely modeled, the influences of these components were 
approximated by equivalent inflow and outflow boundary conditions. Figure 10.1 
illustrates axisymmetric projections of the resulting EEE mesh/geometry model. 

The resulting primary mesh for the EEE LP analysis consisted of 74 separate 
blocks and approximately 6.7 million grid points. The meshes and corresponding 
boundary data file were sequentially “coarsened” by removing every other grid point 
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Figure 10.1: Axisymmetric projection of Energy Efficient Engine (EEE) Low Pressure 
(LP) Subsystem analysis component layout and mesh system. 
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in each comput ational coordinate direction to yield a mesh system of of 900,000 points. 
This “coarse” mesh was employed to debug the overall solution specifications, and to 
ease the difficulties encountered with numerical initialization of the solution caused 
by the high pressure ratios encountered in a complete engine solution. It was also 
found that the solution could be effectively initialized by employing the individual 
component solution data obtained from the component validation studies. This type 
of data would not normally be available for the analysis of a new engine design, and 
it was therefore deemed important to be able to demonstrate that the solution could 
be initialized from an arbitrary initial specification. 


10.2 EEE LPS Processing Strategy 

All calculations for the EEE LP Subsystem were performed on parallel computing 
systems. Four such systems were described in detail in Chapter 6. At the very end 
of this program, timing comparisons were also obtained on a 128-processor Silicon 
Graphics Origin 2000 computer. All parallel calculations employed the MPI [114] 
message passing specification, with the primary programming sublayer coded using 
the the APPL [112] message passing library. Conversion between the APPL specifi- 
cation and the MPI specification was handled by a conversion library (APPLMPI). 
This layered coding structure is outlined in more detail in the ADPAC reference 
manuals [107], [142]. Inter-processor communication based on the MPI programming 
specification was handled using two different MPI libraries. Initially, a public domain 
MPI library referred to as MPICH was employed as this package was self-contained, 
had an automatic configuration script, and was available for a very wide variety of 
computing platforms. During the latter stages of this project, Silicon Graphics Corp. 
developed a proprietary MPI implementation (SGI MPI 3.0) which was also employed 
during this study. Timing comparisons for the various computational platforms and 
communication libraries employed in this study for the EEE LP Subsystem analysis 
are provided in Table 10.1. 

Several interesting conclusions can be drawn from this comparison. Solution times 
varied widely based on computing platform. Overall computation time (wall clock 
time in this instance) is governed by essentially three factors: system load, processor 
load, and communication load. For each of the times presented in Table 10.1, every 
effort was made to perform the timing study on an unloaded system. That is, the 
system, while not necessarily dedicated, was essentially unloaded when the timing 
comparison was performed. This was assumed to eliminate the system load factor 
as a significant contributor to the overall time. The remaining time was therefore 
essentially a function of CPU load and communication load. CPU load was con- 
trolled through the block/processor assignment algorithm employed by the ADPAC 
analysis. The ADPAC code performs parallel computations via a domain decompo- 
sition coarse grained computing strategy. The division of the computational effort 
is accomplished by assigning one or more blocks of the multiple block mesh to spe- 
cific processors. This assignment can be directly specified by the user, or through 
the code predefined assignment strategy. The overall processing load for a given 
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Wall Clock Time Summary 

(100 Iterations of EEE/LP Model) 



Coarse Mesh 
Number of Processors 

Fine Mesh 

Number of Processors 


8 

16 

32 

8 

16 

32 

LACE 

communication 

5380 

2139 

2707 

N/A 

23063 

_ 

total solver time 

7762 

4846 

4198 

N/A 

77427 

— 

Babbage 

communication 

952 

403 

735 

— 

— 

8763 

total solver time 

2673 

1418 

1089 

— 

— 

17518 

Davinci 

communication 

- 

- 

- 

4617 

- 

- 

total solver time 

- 

- 

— 

18122 

— 

— 

Allison SGI 
Power Challenge 

communication 

585 

182 

N/A 

N/A 

N/A 

N/A 

total solver time 

1278 

673 

N/A 

N/A 

N/A 

N/A 

SGI 

Origin 2000 

communication 

268 

153 

264 



3105 

total solver time 

781 

403 

327 

- 

- 

5528 


N/A - not applicable (machine resources insufficient to performing the operation) 

LACE: NASA Lewis IBM RS-6000 cluster. 

Babbage: NASA Ames IBM SP2 cluster 

Davinci: NASA Ames SGI cluster 


Table 10.1: Tabulation of parallel computing CPU time estimates for platforms em- 
ployed for the EEE LP Subsystem analysis (all times given are wall clock time on 
non-dedicated systems with precautions taken to eliminate outside loading factors). 
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processor is based on the total number of mesh cells contained within the blocks as- 
signed to that processor. Thus, for a mesh system with widely varying mesh sizes, 
optimizing the block/processor assignment to balance the processor computational 
load can be a difficult task. This was exactly the case with the EEE LPS analysis. 
The smallest single mesh block in the EEE LPS mesh system contained 25 compu- 
tational cells while the largest single mesh block contained 448,497 computational 
cells. This disparity in block size and the flexibility in the number of processors and 
the block/processor assignment strategy makes balancing the overall computational 
load in the parallel computing environment a very difficult proposition. Table 10.1 
illustrates CPU time estimates based on three different numbers of processors (8, 16, 
32). In each case, the block processor assignment strategy was to attempt to balance 
the computational load. Xo real bias was devoted to incorporating the communica- 
tion overhead in the block/processor load balancing strategy. The block/processor 
assignment was not necessarily considered optimal, but should be reasonable in terms 
of providing a good estimate of the type of parallel computing performance which can 
be achieved in a production environment. The final factor controlling overall CPU 
time is the inter-processor communication load. The communication load, in turn, 
is governed by many factors including system hardware, communication library, and 
block/processor assignment. This area is often the limiting factor in determining 
the total number of processors which can be effectively applied to a large-scale CFD 
simulation. As more processors are added to attack a given problem, the individual 
CPU load goes down, while the communication requirements go up. The experi- 
ence gained in this study suggests that for the current status of computer equipment 
(processor power, communication speed) a near optimal arrangement for multistage 
turbomachinery calculations was achieved when 1-2 processors was assigned for each 
blade row in the machine. For the complete EEE LPS simulation, this level was not 
acheived on every system tested as there were 15+ blade rows in every simulation, 
and several systems were limited to a maximum of 16 processors. 

Overall, the following comments can be made concerning the parallel performance 
studies: 

• Peak processing speed was acheived on a Silicon Graphics Origin 2000 using the 
SGI MPI 3.0 communication library. 

• Estimated turnaround time for a single operating point was estimated to be 10 
hours on the SGI Origin 2000 system using 32 processors. 

• Load balance was non-optimal for the present mesh configuration. It seems 
entirely possible that significant improvements in parallel computing efficiency 
might be achieved through a more structured specification of mesh block di- 
mensions in the overall problem. 

• For the faster systems, parallel computing efficiency was still nearly linear with 
the addition of more processors. This implies that the problem could still be 
effectively acclerated if systems with larger numbers of processors (> 100) were 
available. 
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10.3 EEE LPS Design Point Simulation 

Preliminary solutions for the EEE LP Subsystem model focused on demonstrating 
the solution convergence behavior and accuracy for engine design point analysis. For 
this simulation, the HP compressor and combustor were modeled through the ap- 
propriate boundary specifications for the engine HP compressor inlet and combustor 
exit planes. The boundary specifications were based on a design point engine cycle 
analysis derived from results from the XEPP computer code. Note that for this set 
of results, the HP turbine (normally considered a core, or HP subsystem component) 
was employed in the CFD model to permit a more reasonable specification of the 
spanwise flow profiles entering the LP turbine. Subsequent large-scale simulations of 
the LP Subsystem did not employ the CFD representation of the HP turbine as it was 
ultimately demonstrated that the LP turbine performance is relatively insensitive to 
inlet flow profile. 


10.4 EEE LPS Shaft Power Balance 

An important aspect of engine simulation, compared to component simulation, is 
that the mating of components often involves both aerodynamic and mechanical cou- 
plings. This concept is illustrated for both single-spool and twin-spool gas turbine 
engines in Figure 10.2. This concept is commonly employed in cycle deck analyses 
(e.g. NEPP ) for components connected by a common shaft. The same concept can 
be applied to larger-scale simulations by providing the appropriate aerodynamic con- 
sistency between components (mass flow, etc.) as well as equating the overall power 
requirements for common shaft-mounted components. This balance was iteratively 
achieved in the present simulation through an iterative procedure which employed 
shaft rotational speed as the means of achieving the desired shaft power balance. 

A series of solutions for the EEE/LP Subsystem was obtained for fixed shaft 
rotational speeds. For each shaft speed, computed power and torque for the rotating 
components were integrated for the rotating components of both the LP turbine 
and fan/booster-stage assemblies. Differences between the computed power/torque 
requirements for the fan and LP turbine assemblies were then employed to estimate 
a new shaft speed for the subsequent solution. Simple physical reasoning suggests 
that if there is power excess, then the shaft speed should increase, and if there is a 
power deficit, then the shaft speed should decrease. A simple linear interpolation was 
employed to estimate the updated shaft speed based on the integrated results from 
two previous solutions. 

A portion of the iterative history of the ADPAC EEE LP shaft power balance is 
given in Table 10.2. As the shaft speed was reduced, the power required by the fan 
was reduced, while the power provided by the LP turbine increased. Eventually, these 
two power levels were essentially identical. The balance was deemed converged when 
the power balance was within 1%. Note that in spite of the changes to the LP system, 
the HP turbine power was relatively constant. This is essentially a result of the fact 
that the core performance was fixed during the shaft power balance procedure. The 
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Aero-Mechanical Coupling in a Single Spool Gas Turbine 



Aero-Mechanical Coupling in a Dual Spool Gas Turbine 

LP Compressor LP Turbine 



LP Shaft 

Mechanical Coupling via Shafts (Separately for LP/HP) 

Aerodynamic Coupling via Flowpath (Joint LP/HP) 

Figure 10.2: Illustration of aerodynamic/mechanical balance required for single-spool 
and twin-spool gas turbine engines. 
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ADPAC Solution (Fixed RPM/Fixed Core, ft-lbf/sec) 


Shaft RPM Fan 


LP Turbine HP Turbine 


3507 

3407 

3250 

3200 


8 , 622,000 

8 , 028,000 

7 , 522,300 

7 , 243,100 


6 , 947,200 

7 , 002,600 

7 , 301,400 

7 , 322,700 


12 . 634.000 

12 . 619.000 

12 . 557.000 

12 . 547.000 


NEPP Solution (Design Point) 

Shaft RPM Fan LP Turbine HP Turbine 

3538.5 8 , 182,000 8 , 182,000 11 , 625,300 

Table 10.2: Tabulation of coarse mesh EEE LP Subsystem shaft power balance iter- 
ative results. 


absolute power levels must be interpreted with the limitations of the CFD analysis in 
mind. The analysis was performed with a constant specific heat, when in fact, given 
the range of temperatures in the machine, the specific heat actually varies up to 5%. 
In addition, parasitic losses in the compressor (endwall leakages, cavity flows, etc.) 
have not been included in the analysis. The shaft power balance also assumes a 100% 
transmission efficiency (no bearing losses). These solutions were typically not run to 
full convergence as only an indicator of the level of power balance was required for the 
intermediate solutions during the balancing procedure. The primary intent at this 
point was to validate the convergence of the shaft power balance iterative process, 
and not necessarily isolate all of the individual features of the problem. 

Also tabulated on Figure 10.2 are the corresponding power estimates from the 
NEPP cycle deck analysis for the design point. The most glaring discrepancy between 
the ADPAC and XEPP results is that the power generated by the HP turbine is lower 
and the power generated by the LP turbine is higher than the corresponding ADPAC 
predictions. It is clear that the predicted LP turbine power output is rather low 
compared to the XEPP cycle analysis data. This is believed to be due to the shrouded 
rotor cavity model applied for the LP turbine in the EEE/LP Subsystem analysis. 
It was demonstrated in the component validation study that the present shrouded 
rotor cavity model can result in a 3%-5% reduction in turbine efficiency, which would 
explain much of the noted discrepancy. In addition, it should be noted that the data 
obtained from the NEPP cycle model represents operation along constant operating 
lines , while the ADPAC simulation is obtained along fixed speed lines. This subtle 
difference is illustrated in Figure 10.3. 

Other potential influences may be due to cooling flows or specific heat ratio vari- 
ations. 

The convergence of the shaft power balance is something of a milestone effort, rep- 
resenting solutions which are both aerodynamically consistent (within the limitations 
of the CFD model, of course), and mechanically consistent. 
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Chapter 11 

ADPAC/NEPP Engine Analysis 

11.1 ADPAC/NEPP Coupling Procedure 

Following the completion of the effort to develop an LP Subsystem shaft power balance 
computational procedure, the logical next step in the LP Subsystem analysis was 
to couple the 3-D ADPAC predictions with a lower order (cycle deck) analysis of 
the core component performance. This coupling is consistent with the “zooming” 
philosophy inherent in the XPSS system architecture. In the present application, 
the core cycle model was based on predictions from the NEPP code. In order to 
incorporate the NEPP results in a systematic fashion, the various interactions between 
the XEPP core model and the ADPAC LP Subsystem model must be addressed. 
One interpretation of these interactions is outlined schematically in Figure 11.1. The 
specifications required from XEPP for the ADPAC analysis are an estimate of the 
core compressor inlet flow (represented initially by a static pressure which is used to 
set the flow in the ADPAC solution), and a specification of the HP turbine inlet total 
pressure and total temperature profiles describing the flow out of the EEE combustor. 
The specifications required from the ADPAC analysis for the XEPP analysis include 
the core compressor inlet total pressures, temperature and velocities (which result 
from the CFD analysis of the fan section). Intertwined in this cross specification is 
the fact that the LP shaft RPM may change as the overall solution evolves, and the 
level and frequency by which the exchanged boundary data between the two analyses 
occurs may be critical. It would also be useful if the ADPAC LP Subsystem analysis 
solution could be initiated based on cycle predictions from the XEPP code. This 
would essentially eliminate the complex solution initialization process required for 
ADPAC analysis for this complex problem. This procedure was, unfortunately, not 
available for the current set of calculations. 

The computational system resulting from the combined NEPP/ ADPAC compu- 
tational procedure is illustrated graphically in Figure 11.2. Since this procedure was 
designed for demonstration purposes, the coupling between the ADPAC and NEPP 
analyses was controlled by a UXIX shell script which sequentially applied the analyses 
in an iterative fashion. Following the application of each analyses, the appropriate 
flow information was extracted from output files by hardwired programs developed 


NASA/CR— 1998-206597 


91 




ADPAC/NEPP Interaction 


A. Use NEFP output to create ADPAC Boundata and Restart files 

1 . No direct link between NEPP and APAC 

2. NEPP output file updates EXIT and INLET conditions and RPM boundata file 

3. NEPP output file creates ADPAC restart file to aid in "start-up" process 

B. Create coding changes in NEPP to perform A 

C. Use SYSTEM call in ADPAC boundata file to run NEPP from ADPAC 

D. Use front-end script to control NEPP- ADPAC interations 


Figure 11.1: Coupled ADPAC/XEPP analysis schematic data flow representation. 
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LP Subsystem Model 
(Color Components/External Flow) — , 
(ADPAC CFD Analysis) 



Dual Spool 
Concentric 
Shaft Arrangement 


_ Core Components 
(Grey Color) 
NEPP Cycle Analysis 


Predicted EEE LP Subsystem Surface Static Pressure 



Figure 11.2: Illustration of coupled ADPAC/XEPP prediction for the EEE LP Sub- 
system (color contours indicate predicted static pressure ratio: red- 10.0, blue-0.36, 
grey scale components are represented by the XEPP cycle analysis). 
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specifically for these two codes, and based specifically on the format of the output for 
each codes. This was, unfortunately an inflexible system, but did have the advantage 
that it could be assembled rather quickly to demonstrate the overall concept. 

A solution for the EEE/LP Subsystem using the coupled ADPAC/ NEPP solution 
strategy was obtained for the design operating point. Problems encountered during 
the initial tests of the solution procedure were traced to excessive variations in the 
boundary specifications during the initial phases of the calculations. These excursions 
were modulated using a simple under-relaxation procedure. The behavior of the 
overall solution procedure was then relatively stable, albeit very slow. Individual 
ADPAC solutions acquired during the iterative cycle can take up to 8 hours on a 
parallel system, with some 10-20 iterations required to achieve complete coupling 
between the ADPAC and XEPP analyses. 

It should be noted that the present demonstration did not employ the LP shaft 
power balance procedure which would be essential to complete the coupled solution 
procedure. At this point, a demonstration of the concept was considered of primary 
importance. The capability demonstrated through this exercise validates the XPSS 
primary objective of “ zooming ' and can hopefully lead to further research in em- 
ploying this type of analysis for future gas turbine engine studies. 
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Chapter 12 

CONCLUSIONS 


Computational Fluid Dynamics (CFD) analysis of the complete Low Pressure (LP) 
Subsystem of the General Electric Energy Efficient Engine (EEE) was demonstrated. 
This study identified several important topical areas to consider in the planning and 
execution of large-scale simulations of complete gas turbine engine propulsion sys- 
tems. The topical areas include geometry manipulation, mesh generation, solution 
initialization, application of parallel computing, full-scale engine simulation, and in- 
terpretation of computational results. Each area is discussed in the sections below. 


12.1 Geometry Manipulation 

This study served to identify both the strengths and the weaknesses of the consistent 
geometry database representation strategy. The NASA EEE Master Engine Geome- 
try Database package evaluated and updated during this study consists of a collection 
of IGES curve-based and surface-based entities which define the primary flowpaths 
and bladed elements in the EEE engine. The component validation study served to 
validate many of these geometry representations through both the mesh generation 
process and the comparison of CFD prediction with test data. The contributions to 
the geometry database elements derived from this study included the following: 


• A one degree (open) reset was applied to the LP turbine first stage vane, based 
on both computational results and discussions with engineers from General 
Electric Corporation familiar with the design. 

• The construction of an axisymmetric representation of the LP turbine shrouded 
rotor endwall seal cavities was also performed. These geometry elements were 
const ructed by hand based on int erpret at ions of published drawings of the act ual 
test rig hardware. 

• Some of the original geometry elements in the Master Engine Geometry Database 
package were mislabeled (rotors and stators swapped). These corrections were 
subsequently reported to NASA so that the database could be updated. 
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Shortcomings of the IGES database geometry representation identified during this 
study included the following: 


• A lack of variable geometry capability hinders the analysis of HP compression 
systems. Specifically, the rotational axes of variable setting stators should be 
defined to permit adjustment representative of the actual compreasion system 
geometry under off-design operating conditions. 

• Some identification of engine secondary flow systems such as customer bleed 
and cooling flow bleed from the HP compression system, coolant flow injection 
in the HP turbine, and drum/cavitv purge flows such as found in both HP and 
LP turbine should be considered. 

• Variations of the geometry with temperature should be identified in some fash- 
ion. At present, it is not well understood whether the geometry elements aie 
the hot or cold (as manufactured) representations. 

• The ability to model statistical variations in the geometry such as might occur 
due to manufacturing tolerances or erosion would also be useful. 

All of these elements point to the need for a flexible geometry manipulation tool 
which could act in concert with the mesh generation/CFD procedures described in 
this report. 


12.2 Mesh Generation 

A relatively simple mesh generation procedure was established during this project 
which employed the geometry database described in Chapter 3. Based on the IGES 
entity blade and flowpath definitions, sheared H-type mesh systems could be rapidly 
generated for multistage compressor or turbine flows. It has been acknowledged that 
the mesh systems are non-optimal in the sense that orthogonality and mesh aspect 
ratio are somewhat compromised for the convenience and simplicity of this neaily 
automated procedure. The analyses described in this report are in no way limited to 
this type of mesh system. Complete automation of the mesh generation pioceduie 
for arbitrary engine configurations would be a significant accomplishment, and was 
beyond the scope of this study. Some estimates of the sensitivity of the solution to 
the mesh density was afforded through the examination of results from “coarsened” 
mesh systems derived by eliminating every other point from an original “fine” mesh. 
It is unlikely that the solutions presented for the LP Subsystem simulation were 
mesh independent as the number of grid points was typically minimized to reduce 
the overall solution computation time. A full IGES-compatible parser could also be 
added to the SEARCH program (the current version of the program is limited to 
specific IGES entity designations). 
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12.3 Solution Initialization 


Some general comments regarding solution initialization seem appropriate given the 
magnitude of the computational effort and the numerical difficulties presented during 
this large-scale simulation. The ideal solution initialization procedure would be to 
have a collection of 3-D CFD isolated component analyses derived with appropriate 
boundary conditions such that an overall representation of the flow could be gen- 
erated by simply assembling the individual component data. This procedure was, 
in effect, verified following the component validation study. Another reasonable ap^ 
proach would be to employ 2-D engine simulations from which an axisvmmetric 3-D 
solution initialization field could be established. Further down the ladder defined by 
this type of hierarchy would be extending 1-D or cycle deck simulations to the 3-D 
space through some sort of interpolation system. Naturally this interplation is some- 
what arbitrary as the extension of lower order data to the higher order system has 
many possible solutions. Finally, least desirable is initialization of the solution from 
an essentially meaningless initial condition (uniform flow). It was demonstrated that 
it is possible to generate the LP Subsystem simulations from the lowest order ini- 
tialization routine, but that this process required a “stair-stepping” of the boundary 
piessuie and tempeiatuie specifications in order to avoid overwhelming the simula- 
tion with nonrealistic pressure or temperature ratios. The development of automated 
couplings between the 3-D and lower order analyses of the types described above (2-D 
and 1- D/cycle analyses) would afford a great simplification in the solution initializa- 
tion pioceduie, and also accelerate the generation of results when evaluating a new 
operating condition. 


12.4 Application of Parallel Computing 

Parallel computing constructs were used extensively during this project, and included 
architectures based on multiprocessor shared-memory computers, to distributed mem- 
ory, network-connected workstation clusters. The analysis was demonstrated under 
foui diffeient parallel computing environments of both NASA and industry origin. 
Oveiall, the best performance was achieved by assigning approximately one blade 
row per processor in the parallel computing environment, as this afforded the best 
compromise between processor load and communication overhead. Parallel comput- 
ing efficiencies on the order of 75% were achieved during this study for the large-scale 
simulations. Load balancing ultimately became the issue which was the greatest ob- 
stacle for improving performance. Careful planning of grid distributions would help to 
i educe this problem, although the large variations in length scales for the components 
in a complete engine simulation make grid distribution a difficult task. 


12.5 Full-Scale Engine Simulation 

Full-scale engine simulations of the GE EEE engine were demonstrated during this 
study based on the 3-D CFD LP Subsystem simulation coupled with the XEPP cycle 
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deck core simulation. Multidisciplinary coupling of common spool components was 
achieved through the application of a shaft power balance. This procedure demon- 
strated both aerodynamic and mechanical coupling of the LP shaft components by 
varying the shaft rotational speed until a mechanical power balance was achieved. 
This requirement is nearly always overlooked in other reported “full-scale” engine 
simulations. The iterative process to achieve the shaft power balance was accelerated 
through the use of the A D PA C mesh sequencing and multigrid capability, which per- 
mits early iterations of the shaft power balance on coarser meshes before proceeding 
with the process on the more computationally expensive finer mesh. 

Direct coupled analyses of the engine operation employing both the 3-D ADPAC 
LP simulation and the NEPP HP/combustor simulation were also successfully de- 
veloped. This procedure demonstrated a two-way coupling between the analyses to 
derive the engine operating condition. This type of analysis features the XPSS zoom- 
ing” concept whereby a portion of the engine is simulated on one level of fidelity, while 
other components are simulated at a different fidelity level. 


12.6 Interpretation of Computational Results 

By and large, all of the computational result developed in this study exhibited qual- 
itative agreement with available data, and did not display any gross violations of 
the expected physical behavior for a given component or collection of components. 
Some comments on the order of accuracy of the analysis are, however, appropriate at 
this time. The mesh sensitivity of the results was discussed earlier, and will not be 
repeated here. Certain aspects of the calculations possessed known errors such as the 
use of non- varying specific heat ratio in the simulations. Other factois contributing 
to errors include omission of windage and mechanical drag in the shaft power balance, 
slight inconsistencies in mass flow due to the 2-D LP turbine shrouded rotor cavity 
model, and the application of adiabatic wall boundary conditions throughout the ma- 
chine. Given these known inconsistencies, detailed evaluation of engine performance 
parameters (i.e., specific fuel consumption (SFC), thrust, etc.) makes little sense, and 
therefore no attempt was made to correlate predicted and experimental values of this 
nature. 


12.7 Recommendations for Future Study 

Given the demonstrated capability of complete engine simulation, a natural extension 
to this work would be to remove the inconsistencies described above, and to proceed 
with the effort to validate the engine operation prediction capabilities of this scheme 
with actual engine data. Large amounts of data exist for simpler, single-shaft en- 
gine configurations which could be used to evaluate the predicted capabilities of this 

scheme. . 

Enhancing the capabilities of the geometry database is also considered a priority. 
Several details related to this effort are described above under the geometry section. 
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Finally, additional studies directed at the evaluation of secondary flow systems 
effects on primary gas flowpath performance would serve to define what level of fi- 
delity is required to incorporate the overall engine performance effect of these “real 
world" components. Simulations of this type have been performed for isolated com- 
pressor airfoils employing inner banded stators [143]. Similar analyses for multistage 
compressor and turbine systems incorporating secondary flow systems models could 
be performed to assess the relative impact of these gas path features. 
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Appendix A 

ADPA C07 Input Files for EEE/LP 
Simulation 


ADPAC Input File for EEE LP Analysis 


*• 


t EEE LP Analysis 

•YARNAME - 

VARIABLE VALUE 

COMMENT 

CASES AKE 

EEELP 

Casa naaa usad for fila naming 

FBCWARN 

1.0 

Turn on warning for B.C. arrors 

RMACH 

0.800000 

Rafarsnca Mach Numbar 

FINYVI 

1.000000 

Viscous triggar (On-1) 

CAMMA 

1.400000 

Spacific hast ratio 

PREF 

769.052800 

Rafaranca Total Prassur* (lbf/ft“2) 

TREE 

444.319200 

Rafaranca Total Tamparatura (dag. A) 

AGAS 

1716.350700 

Cas constant 

DIAH 

0.083333 

Rsfaranca langth to convart grid to faat 

EPSX - 

1.000000 

Rasidual smoothing multipliar 

EPSY 

0.600000 

Rasidual smoothing multipliar 

EPSZ 

1.000000 

Rasidual smoothing multipliar 

VIS2 

0.500000 

2nd ordar dissipation coafficiant 

VIS4 

0.015625 

4th ordar dissipation coafficiant 

CF1 

-6.000000 

Tima stap multipliar (CFL numbar) 

FNCHAX 

300.00000 

Numbar of itarations on fina grid 

FITCHK 

6.000000 

Chackpoint rastart itaration intarval 

FTIMEI 

1.000000 

Tima stap calculation intarval 

FTUABI 

1.000000 

Turbulanca modal updata intarval 

FTUABB 

20.000000 

Turbulanca modal initiation itaration numbar 

PANO 

0.700000 

Prandtl numbar 

PRTNO 

0.900000 

Turbulant Prandtl Numbar 

FSOLYE 

1.000000 

Solution schama triggar (1 • 4-staga schama) 

FRESID 

1.000000 

Rasidual smoothing triggar (0n-l) 

FREST 

1.000000 

Rastart flag (-1, rastart using casa .rastart .old) 

P3DPRT 

1.0 

PL0T3D fila output triggar 

FUN I NT 

99999.000000 

Unstaady output itaration intarvml 

FMULTI 

2.000000 

Numbar of srultigrid lavals, -l,no multigrid 

FSUBIT 

3.000000 

Numbar of subitarations during multigrid 

FFULMC 

0.000000 

Full multigrid triggar (start on coarsar mash-1) 

FC0AC1 

2.000000 

Full multigrid starting mash laval 

FC0AC2 

2.000000 

Full multigrid anding mash laval 

FITFMC 

200.000000 

Full multigrid itarations (on coarsa mashas) 

VISCG2 

0.126000 

Coarsa mash dissipation coafficiant 

FCRAFIX 

6.000000 

Intaractiva graphics triggar (0-off) 

FCRAFINT 

1.000000 

Intaractiva graphics updata intarval 

FXMCSAV 

0.000000 

Intaractiva graphics scraansava triggar 

FIHCINT 

99999.000000 

Intaractiva graphics scraasava intarval 

FVTSFAC 

7.000000 

Viscous tima stap factor 

FTOTSM 

1.000000 

Hultigrid smoothing triggar (on-1.0) 

EPSTOT 

0.100000 

Multigrid smoothing coafficiant 

FWALLF 

1.000000 

Wall Function triggar (on-l) 

CFKAX 

2.200000 

Raf. CFL # for implicit rai. smoothing 


I > Sat tha block rotational spaads hara 


• Upstraam outar nacalla 

• 

RPM(l) - 0.000 

APH(2) - 0.000 


• Fan + Quartar staga boo* tar 
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RPM(3) 


3507.000 

RPN(4) 


3607.000 

RPM(5) 


0.000 

RPM(6) 


3607.000 

RPM (7) 


0.000 

RPH(8) 


0.000 

RPM(9) 


0.000 

RPN(lO) 


0.000 

« 



0 Bypass 

duct 


I 



RPM(ll) 

- 

0.000 

0 



• HP turbine 


C 



RPH(12) 

- 

0.000 

RPM( 13) 

- 

12627.316 

RPM(14) 

• 

0.000 

RPM(16) 

- 

12627.316 

• 



• LP turbine 


0 



RPM(IC) 


0.000 

RPM(17) 


3607.000 

RPM(lfl) 


0.000 

RPM(19) 


3607.000 

RPM (20) 


0.000 

RPN(21) 


3607.000 

RPM (22) 


0.000 

RPM (23) 


3607.000 

RPM (24) 


0.000 

RPH(25) 


3607.000 

• 



t LP turbine shrouded rotor cavities 

« 



RPM (26) 


0.000 

RPHC27) 


0.000 

RPH(26) 


0.000 

RPH(29) 


0.000 

RPH(30) 


0.000 

RPM(31) 


0.000 

RPM (32) 


0.000 

RPM (33) 


0.000 

RPM (34) 


0.000 

RPM (36) 


0.000 

RPM (36) 


0.000 

RPM (37) 


0.000 

RPM (38) 


0.000 

RPM (39) 


0.000 

RPM (40) 


0.000 

RPH(41) 


0.000 

RPM (42) 


0.000 

RPM(43) 


0.000 

RPM (44) 


0.000 

RPM (46) 


0.000 

RPM (46) 


0.000 

RPM(47) 


0.000 

RPM(48) 


0.000 

RPM (49) 


0.000 

RPM (60) 


0.000 

RPM(61) 


0.000 

RPM (62) 


0.000 

RPM (63) 


0.000 

RPM (64) 


0.000 

RPM (66) 


0.000 

RPM (66) 


0.000 

RPM(67) 


0.000 

RPM (68) 


0.000 

RPM(69) 


0.000 

RPM (60) 


0.000 

RPM(61) 


0.000 

RPM (62) 


0.000 

RPM (63) 


0.000 

RPM (64) 


0.000 

RPM (66) 


0.000 

• 



0 Lobed 

exhaust 

ssixer and downs treaa 

0 



RPM (66) 


0.000 

RPM(67) 


0.000 

RPM(68) 


0.000 

RPM(69) 


0.000 

RPM (70) 


0.000 

RPM(71) 


0.000 

RPM (72) 


0.000 

RPM(73) 


0.000 
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RPM(74) 


0.000 


« 

i > Set the number of blades for each block where possible 

* 


0 

0 Upstream outer nacelle 


* 



NBLD(l) 

- 

1.000 

NBLD (2) 

“ 

1.000 

0 



• Fan ♦ 

quarter 

stage booster 

0 



NBLD (3) 


32.000 

NBLD (4) 


32.000 

NBLD (6) 


60.000 

NBLD(6) 


66.000 

NBLD(7) 


64.000 

NBLD(B) 


64.000 

NBLD(9) 


34.000 

NBLD(IO) 


34.000 

0 



0 Bypass 

duct 


0 



NBLD (11) 

“ 

1.000 

0 



0 HP turbine 


0 



NBLD (12) 

- 

46.000 

NBLD(13) 

- 

76.000 

NBLD(14) 

- 

48.000 

NBLD(15) 

* 

70.000 

0 



0 LP turbine 


0 



NBLD(16) 


72.000 

NBLD (17) 


120.000 

NBLD (18) 


102.000 

NBLD (19) 


122.000 

NBLD (20) 


96.000 

NBLD (21) 


122.000 

NBLD (22) 


114.000 

NBLD (23) 


156.000 

NBLD (24) 


120.000 

NBLD (25) 


110.000 


* 

0 LP turbine shrouded rotor cavities 


NBLD (26) 
NBLD (27) 
NBLD (26) 
NBLD (29) 
NBLD (30) 
NBLD(31) 
NBLD(32) 
NBLD (33) 
NBLD (34) 
NBLD (35) 
NBLD (36) 
NBLD (37) 
NBLD (38) 
NBLD (39) 
NBLD (40) 
NBLD (41) 
NBLD (42) 
NBLD (43) 
NBLD (44) 
NBLD (46) 
NBLD (46) 
NBLD (47) 
NBLD (46) 
NBLD (49) 
NBLD (60) 
NBLD(61) 
NBLD (62) 
NBLD (63) 
NBLD (64) 
NBLD (66) 
NBLD (66) 
NBLD (67) 
NBLD (68) 
NBLD (69) 
NBLD (60) 
NBLD (61) 
NBLD (62) 
NBLD (63) 


1.000 

1.000 

1.000 

1.000 


1.000 


1.000 


1.000 

1.000 

1.000 

1.000 

1.000 


1.000 


1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 


NASA/CR— 1998-206597 


115 



NBLDC64) 

NBLD(65) 


1.000 

1.000 


i Lob«d exhaust mixer and downstream 


I 

NBLD(66) 

NBLDC67) 

m0 (66) 

mo (69) 
NBLDC70) 
mD(71) 
mD(72) 
NBLD(73) 
KBLD(74) 


12.000 

12.000 

12.000 

12.000 

12.000 

12.000 

12.000 

12.000 

12.000 
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ADPAC Boundary Data File for EEE LP Analysis 

» 

i Block dascription* for boundary condition* 

a 


t BLOCX 1 - upstraam and outar nacalla 

mash 

(2D) 

• BLOCX 2 - upstraam and outar nacalla 

mash 

(2D) 

t BLOCK 3 - quart ar stag a and fan 

• 

f 

mash 

(3D) 

* 

i BLOCX 10 - quart ar stag a and fan 

mash 

(3D) 

• BLOCK 11 - bypass duct 

mash 

(2D) 

• BLOCK 12 - high prassura t urbina 

« 

§ 

mash 

(3D) 

1 

V BLOCK 15 - high prassura turbina 

mash 

(3D) 

• BLOCK 16 - low prassura turbina 
» 

mash 

(3D) 

» 

• BLOCK 66 - low prassura turbina 

mash 

(3D) 

# BLOCK 66 - mixar and downstraam 

• 

mash 

(3D) 

S 

V BLOCK 74 - mixar and downstraam 

• 

mash 

(3D) 


9 UPSTREAM AND OUTER NACELLE CRIDDINC 

PATCH 1 2JJPMIK 1 101 1 66 1 2 1 65 1 2 

PATCH 2 1JJKPIK 101 1 166 1 2 1 66 1 2 

• Uppar part of inlat 

FREE 1 1IIPPJX 1 1 1 49 1 2 1 49 1 2 

PTOT TTOT EM INF ALPHA 

1.0 1-0 0.8 0.0 

• Nacalla top 

SSVI 1 1 J J P P I X 1 1 66 129 1 2 66 129 1 2 

RPMVALL TV ALL 

0.0 0.0 

• 2D linkup w/ 3D mixar (abova nacalla) 

KBCAVC 1 70 I I M P J X 129 1 1 49 1 2 1 49 1 81 

MSECS 
1 

LBLDCX2B LFACE2B LDIR2B L2LIMB M2LIM1B K2LIM2B N2LIM1B N2LIM2B 

70 I P 1 1 49 1 81 

KBCAVC 70 1 I I P M J X 1 129 1 49 1 81 1 49 1 2 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIMB M2LIM1B M2LIM2B N2LIM1B N2LIM2B 

1 I M 129 1 49 1 2 

• Top nacalla farfiald 

FREE 1 1JJHMIK 49 49 1 129 1 2 1 129 1 2 

PTOT TTOT EMINF ALPHA 

1.0 1.0 0.8 0.0 

t Lovar part of inlat 

FREE 2 2IIPPJK 1 1 1101 1 2 1101 1 2 

PTOT TTOT EMINF ALPHA 

1.0 1.0 0.8 0.0 

f 2D linkup s/ 3D fan (insida inlat) (lovar saction) 

KBCAVC 2 31 IKPJX93 1 181 1 2 181 149 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIMB K2LIM1B K2LIM2B N2LIK1B N2LIK2B 

3 I P 1 1 81 1 49 

KBCAVC 3 2IIPMJX 1 93 1 81 1 49 1 81 1 2 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIK3 H2LIM1B M2LIK2B N2LIM1B N2LIM2B 

2 I H 93 1 81 12 

f 2D linkup w/ 3D fan (insida inlat) (uppar saction) 

KBCAVC 2 4 I I M P J X 93 1 81 101 1 2 1 21 1 49 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIMB M2LIH1B H2LIM2B N2LIM1B N2LIM2B 

4 I P 1 1 21 1 49 

KBCAVC 4 2 I I P M J X 1 93 1 21 1 49 81 101 1 2 

NSECS 
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1 

LBL0CX23 LFACE2B LDIR2B L2LIKB M2LIM1B M2LIM2B N2LIH1B N2LIH2B 
2 I H 93 81 101 1 2 


i Axis of symmetry 
SSIN 2 2 J J 

P P 

I X 

1 1 

1 

93 

1 2 

1 

93 

1 2 

RPHVALL 

TV ALL 










0.0 

0.0 










9 Inside nacelle 

SSVI 2 2 J J 

M M 

I K 

101 101 

65 

93 

1 2 

66 

93 

1 2 

APMVALL 

TV ALL 










0.0 

0.0 











! 

m 

STACE AND 

FAN 















PATCH 

3 

4 

J 

J 

K 

P 

I 

X 

81 

1 

1 

67 

1 

49 

1 

67 

1 

49 

PATCH 

4 

3 

J 

J 

P 

M 

I 

X 

1 

81 

1 

67 

1 

49 

1 

67 

1 

49 

PATCH 

3 

4 

J 

J 

M 

P 

I 

X 

81 

l 

81 

113 

1 

49 

81 

113 

1 

49 

PATCH 

4 

3 

J 

J 

P 

M 

I 

X 

1 

81 

81 

113 

1 

49 

81 

113 

1 

49 

PATCH 

3 

3 

X 

X 

P 

M 

I 

J 

1 

49 

1 

33 

1 

81 

1 

33 

1 

81 

PATCH 

3 

3 

X 

X 

M 

P 

I 

J 

49 

1 

1 

33 

1 

81 

1 

33 

1 

81 

PATCH 

3 

3 

X 

X 

P 

K 

I 

J 

1 

49 

97 

113 

1 

81 

97 

113 

1 

81 

PATCH 

3 

3 

X 

X 

M 

P 

I 

J 

49 

1 

97 

113 

1 

61 

97 

113 

1 

81 

PATCH 

4 

4 

X 

X 

P 

M 

I 

J 

1 

49 

1 

33 

1 

21 

1 

33 

1 

21 

PATCH 

4 

4 

X 

X 

M 

P 

I 

J 

49 

1 

1 

33 

1 

21 

1 

33 

1 

21 

PATCH 

4 

4 

X 

X 

P 

M 

I 

J 

1 

49 

97 

113 

1 

21 

97 

113 

1 

21 

PATCH 

4 

4 

K 

X 

M 

P 

I 

J 

49 

1 

97 

113 

1 

21 

97 

113 

1 

21 

MBCAYC 

NSECS 

3 

6 

I 

I 

M 

P 

J 

K 

113 

1 

1 

65 

1 

49 

1 

65 

1 

33 

1 

LBL0CK23 

LFACE23 

LDIR23 

L2LIMB 

M2LIM1B 

M2LIK2B 

N2LIH1B 

N2LIM2B 



5 


I 


P 



1 


1 


65 


1 


33 



MBCAVG 

NSECS 

5 

3 

I 

I 

P 

M 

J 

X 

1 

113 

1 

65 

1 

33 

1 

66 

1 

49 

1 

LBL0CX2B 

L FACE 23 

LDIR2B 

L2LIHB 

M2LIM1B 

M2LIM23 

N2LIM1B 

N2LIK2B 



3 


I 


M 


113 


1 


65 


1 


49 



MBCAVG 

NSECS 

3 

10 

I 

I 

M 

P 

J 

X 

113 

1 

65 

81 

1 

49 

1 

17 

1 

49 

1 

LBL0CK2B 

LFACE2B 

LDIR2B 

L2LIN3 

H2LIM1B 

M2LIM2B 

N2LIM1B 

N2LIK2B 



10 


I 


P 



1 



1 

17 


1 


49 



MBCAVG 

NSECS 

10 

3 

I 

I 

P 

M 

J 

X 

1 

113 

1 

17 

1 

49 

65 

81 

1 

49 

1 

LBL0CX2B 

LFACE2B 

LDIA2B 

L2LIKB 

K2LIM1B 

K2LIM2B 

N2LIM1B 

N2LIM2B 



3 


I 


M 



113 


65 

61 


1 


49 




f before spinner 

SSIN 3 3JJPPIK 1 1 1 9 149 1 9 14 9 

APMVALL TV ALL 

0.0 0.0 

• spinner surface 

SSVI 3 3 J J P P I X 1 1 9 113 1 49 9 113 1 49 

APMVALL TV ALL 

3507.000 0.0 

SSVI 3 3 J J H N I K 81 81 57 81 1 49 67 61 1 49 

APMVALL TV ALL 

3507.000 0.0 

SSVI 3 3 K K P P I J 1 1 33 97 1 81 33 97 1 81 

APMVALL TV ALL 

3507.000 0.0 

SSVI 3 3 K X M M I J 49 4 9 33 97 1 81 33 97 1 81 

APMVALL TV ALL 

3507.000 0.0 

M3CAVC 4 10 I I M P J X 113 1 1 21 1 49 17 37 1 49 

NSECS 
1 

LBL0CX2B LFACE2B LDXA2B L2LIMB M2LIM1B M2LIM2B N2LIM1B N2LIH2B 

10 I P 1 17 37 1 49 

MBCAYC 10 4 I I P K J K 1 113 17 37 1 49 1 21 1 49 

NSECS 
1 

LBL0CX2B LFACE2B LDIA2B L2LIMB M2LIM1B M2LIK2B N2LIM1B N2LIM2B 
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4 I M 113 1 21 1 49 

SSVI 4 4 J J P P I K 1 1 67 81 1 49 67 81 1 49 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 4 4JJHMIX2121 1 113 1 49 1 113 1 49 

RPMVALL TV ALL 

0.0 0.0 


SSVI 

4 4 K 

X P 

P 

I 

J 

1 

1 

33 

97 

1 

21 

33 

97 

1 

21 

RPMVALL 

TV ALL 















3607.000 

0.0 














SSVI 

4 4 K 

X M 

H 

I 

J 

49 

49 

33 

97 

1 

21 

33 

97 

1 

21 

RPMVALL 

TV ALL 















3607.000 

0.0 














PATCH 

6 5 X 

K P 

M 

I 

J 

1 

33 

1 

17 

1 

65 

1 

17 

1 

65 

PATCH 

5 6 X 

K M 

P 

I 

J 

33 

1 

1 

17 

1 

65 

1 

17 

1 

65 

PATCH 

5 5 X 

X P 

M 

I 

J 

1 

33 

81 

97 

1 

65 

81 

97 

1 

65 

PATCH 

6 5 X 

X M 

P 

I 

J 

33 

1 

81 

97 

1 

65 

81 

97 

1 

65 

MBCAVC 

5 6 1 

I M 

P 

J 

X 

97 

1 

1 

66 

1 

33 

1 

66 

1 

33 

NSECS 
















LBL0CK2B 

L FACE 2 3 

LDIR2B 

L2LIMB 

M2LIM1B 

K2LIM2B 

N2LIH1B 

N2LIM2B 



6 

I 

P 


1 


1 


65 



1 

33 



MBCAVC 

6 5 1 

I P 

M 

J 

X 

1 

97 

1 

66 

1 

33 

1 

65 

1 

33 

NSECS 
















LBL0CK2B 

LFACE2B 

LDIR2B 

L2LIM3 

H2LIM1B 

M2LIK2B 

N2LIM1B 

N2LIH2B 



6 

I 

M 


97 


1 


65 



1 

33 



SSVI 

5 5 J 

J P 

P 

I 

K 

1 

1 

1 

97 

1 

33 

1 

97 

1 

33 

RPMVALL 

TV ALL 















0.0 

0 

.0 














SSVI 

6 6 J 

J M 

M 

I 

X 

65 

66 

1 

97 

1 

33 

1 

97 

1 

33 

RPMVALL 

TV ALL 















0.0 

0 

.0 














SSVI 

6 6 K 

X P 

P 

I 

J 

1 

1 

17 

81 

1 

65 

17 

81 

1 

66 

RPMVALL 

TV ALL 















0.0 

0 

.0 














SSVI 

6 6 X 

K M 

H 

I 

J 

33 

33 

17 

81 

1 

66 

17 

81 

1 

65 

RPMVALL 

TV ALL 















0.0 

0 

.0 














PATCH 

6 6 X 

X P 

M 

I 

J 

1 

33 

1 

17 

1 

65 

1 

17 

1 

65 

PATCH 

6 6 X 

K M 

P 

I 

J 

33 

1 

1 

17 

1 

66 

1 

17 

1 

65 

PATCH 

6 6 X 

K P 

M 

I 

J 

1 

33 

81 

97 

1 

65 

81 

97 

1 

65 

PATCH 

6 6 X 

K M 

P 

I 

J 

33 

1 

81 

97 

1 

65 

81 

97 

1 

65 

MBCAVC 

6 7 1 

I M 

P 

J 

X 

97 

1 

1 

33 

1 

33 

1 

33 

1 

33 


NSECS 

l 


LBLDCX23 LFACE2B LDIR2B L2LIMB M2LIM1B M2LIK2B K2LIM1B N2LIM2B 

7 IP1 133 133 

KBCAVC 7 6 I I P M J X 1 97 133 1 33 1 33 133 

NSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIMB M2LIM1B K2LIH2B N2LIM1B N2LIH2B 

6 I M 97 1 33 1 33 

MBCAVC 6 8IIMPJX97 1 33 65 1 33 1 33 1 33 

MSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIMB K2LIM1B K2LIK2B N2LIM1B N2LIK2B 

8 I P 1 1 33 1 33 

MBCAVC 8 6IIPMJX 1 97 1 33 1 33 33 66 1 33 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIMB M2LIM1B H2LIK2B N2LIM1B N2LIM2B 

6 I H 97 33 65 1 33 

SSVI 6 6JJPPIK 1 1 1 97 1 33 1 97 1 33 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 6 6 J J M H I K 65 65 1 97 1 33 1 97 1 33 

RPMVALL TV ALL 

0.0 0.0 

SSVI 6 6XKPPIJ 1 1 17 81 1 66 17 81 166 

RPMVALL TV ALL 
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3607.000 


0.0 


SSVI 

6 

6 

X 

X 

N 

H 

I 

J 

33 

33 

17 

81 

1 

65 

17 

81 

1 

66 

RPMVALL 


TWALL 
















3607. 

.000 



0 . 

0 














PATCH 

7 

8 

J 

J 

H 

P 

I 

X 

33 

1 

1 

17 

1 

33 

1 

17 

1 

33 

PATCH 

8 

7 

J 

J 

P 

H 

I 

X 

1 

33 

1 

17 

1 

33 

1 

17 

1 

33 

PATCH 

7 

7 

X 

K 

P 

M 

I 

J 

1 

33 

1 

26 

1 

33 

1 

25 

1 

33 

PATCH 

7 

7 

X 

K 

M 

P 

I 

J 

33 

1 

1 

25 

1 

33 

1 

26 

1 

33 

PATCH 

7 

7 

X 

X 

P 

H 

I 

J 

1 

33 

89 

105 

1 

33 

89 

106 

1 

33 

PATCH 

7 

7 

X 

X 

M 

P 

I 

J 

33 

1 

89 

105 

1 

33 

89 

106 

1 

33 

PATCH 

8 

8 

X 

X 

P 

H 

I 

J 

1 

33 

1 

29 

1 

33 

1 

29 

1 

33 

PATCH 

8 

8 

X 

X 

M 

P 

I 

J 

33 

1 

1 

29 

1 

33 

1 

29 

1 

33 


• Exit fro* byp*s» (inlet to cooproesor) 

EXITC 7 7 I I N M L L 106 106 1 33 1 33 1 33 1 33 

PEXIT 

1.6242046 1.076 0.02 


SSVI 

7 

7 J 

J 

p 

P 

I 

X 

1 

1 

1 

106 

1 

33 

1 

105 

1 

33 

RPHWALL 


TVALL 
















0.0 


0 . 

0 















SSVI 

7 

7 J 

J 

M 

H 

I 

X 

33 

33 

17 

105 

1 

33 

17 

105 

1 

33 

RPMVALL 


TVALL 
















0.0 


0 . 

0 















SSVI 

7 

7 X 

X 

P 

P 

I 

J 

1 

1 

26 

89 

1 

33 

25 

69 

1 

33 

RPHWALL 


TVALL 
















0.0 


0 . 

0 















SSVI 

7 

7 X 

X 

M 

H 

I 

J 

33 

33 

26 

89 

1 

33 

25 

89 

1 

33 

RPHWALL 


TWALL 
















0.0 


0.0 















HBCAVC 

MSECS 

1 

LB LOCK 23 

8 

9 I 

I 

M 

P 

J 

X 

29 

1 

1 

33 

1 

33 

1 

33 

1 

49 

LFACE2B LDIR2B 

L2LIKB 

H2LIM1B 

M2LIM2B 

N2LIM1B 

N2LIN2B 



9 


I 

p 



1 


1 


33 


1 

49 



HBCAVC 

MSECS 

9 

8 I 

I 

P 

H 

J 

X 

1 

29 

1 

33 

1 

49 

1 

33 

1 
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0.0 


0.0 


SSVI 10 10 J J H M I X 37 37 1 137 1 49 1 137 1 49 

RPMVALL TV ALL 

0.0 0.0 

SSVI 10 10 X X P P I J 1 1 49 121 1 37 49 121 1 37 

RPMVALL TV ALL 

0.0 0.0 

SSVI 10 10 X X H M I J 49 49 49 121 1 37 49 121 1 37 

RPMVALL TV ALL 

0.0 0.0 

• Lovar section 30 qtr-staga to 20 bypass duct 

KBCAVC 9UIIMPJX89 1 133 1 49 1 17 1 2 

MSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIM3 M2LIM1B K2LIM2B N2LIM1B N2LIM2B 

11 I P 1 1 17 12 

MBCAVC 11 9IIPMJK 1 89 1 17 1 2 133 1 49 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIKB M2LIM1B H2LIK2B N2LIM1B N2LIM2B 

9 I M 89 1 33 1 49 

• Uppar saction 30 qtr-stag* to 20 bypass duct 

MBCAVC 10 11 I I M P J X 137 1 1 37 1 49 17 49 1 2 

MSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIMB M2LIM1B M2LIM2B M2LIM1B N2LIH2B 

11 I P 1 17 49 1 2 

MBCAVC 11 10 I I P M J K 1 137 17 49 1 2 1 37 1 49 

MSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIKB H2LIM1B M2LIM2B N2LIM1B M2LIM2B 

10 I M 137 1 37 1 49 


I BYPASS DUCT 

SSVI 11 11 J J P P I K 1 1 1 97 1 2 1 97 1 2 

RPMVALL TV ALL 

0.0 0.0 

SSVI 11 11 J J M H I K 49 49 1 97 1 2 1 97 1 2 

RPMVALL TV ALL 

0.0 0.0 

SSVI 11 11 I I M M J K 97 97 1 49 1 2 1 49 1 2 

RPMVALL TWALL 

0.0 0.0 

* 20 bypass duct to 30 mixar linkup 

MBCAVC 11 73 I I M P J K 97 1 1 49 1 2 1 49 1 81 

MSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIMB M2LIM1B H2LIM2B N2LIM1B N2LIM2B 

73 I P 1 1 49 1 81 

MBCAVC 73 11 I I P M J X 1 97 1 49 1 81 1 49 1 2 

NSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIKB M2LIM1B K2LIM2B N2LIM1B N2LIM2B 

11 I M 97 1 49 1 2 


• HICK PRESSURE TURBINE (HPT) 
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•• Inlet 
•INLETC 

to KPT 
12 12 I 

I P 

P J X 1 

1 

1 33 

1 49 

1 33 

1 49 

* PTOT 
< 35.66 

TTOT 

6.632 








« Van* 1 

•SSVI 

12 12 J 

J P 

P I X 1 

1 

1 97 

1 49 

1 97 

1 49 


• RPHVALL TVALL 

* 0.0 0.0 

•SSVI 12 12 J J H H I X 33 33 1 07 1 49 1 97 1 49 

* RPHVALL TV ALL 

** 0.0 0.0 

•SSVI 12 12 X X P P I J 1 1 17 81 1 33 17 81 1 33 

• RPHVALL TWALL 

« 0.0 0.0 

•SSVI 12 12 X X H H I J 49 49 17 81 1 33 17 61 1 33 

* RPHVALL TWALL 

• 0.0 0.0 

KBCAVG 12 13 I I H P J X 97 1 1 33 1 49 1 33 1 49 

NSECS 
1 

LBL0CX23 LFACE2B LDIR2B L2LIKB H2LIH1B H2LIH2B N2LIM1B N2LIK2B 

13 I P 1 1 33 1 49 

KBCAVC 13 12IIPHJX 1 97 1 33 1 49 1 3 3 1 49 

NSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIMB H2LIH1B M2LIH23 N2LIH1B N2LIM2B 

12 I H 97 1 33 1 49 

• Blade 1 

•SSVI 13 13 J J P P I X 1 1 1 97 1 49 1 97 1 49 

• RPHVALL TWALL 

» 12627.315 0.0 

•SSVI 13 13 J J H H I X 33 33 1 97 1 49 1 97 1 49 

I RPHVALL TWALL 

• 0.0 0.0 


•*»•*• HERE BE THE TEST - ADD HTPERSPACE TIP CLEARANCE 

•SSVI 13 13 X X P P I J 1 1 17 81 1 33 17 81 1 33 

•• RPHVALL TV ALL 
•* 12627.315 0.0 

•• 

•SSVI 13 13 X X H H I J 49 49 17 81 1 33 17 81 1 33 

« RPHVALL TWALL 

• 12627.315 0.0 

• 

•SSVI 13 13 X X P P I J 1 1 17 81 1 29 17 81 1 29 

• RPHVALL TV ALL 

• 12627.316 0.0 

•SSVI 13 13 X X H K I J 49 49 17 81 1 29 17 81 1 29 

• RPHVALL TV ALL 

• 12627.315 0.0 

• 

•PATCH 13 13 X X P H I J 1 49 17 81 29 33 17 81 29 33 

•PATCH 13 13 X X H P I J 49 1 17 81 29 33 17 81 29 33 


HBCAVC 13 14 I I H P J X 97 1 133 1 49 1 33 1 49 

NSECS 
1 
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LBL0CK2B LFACE2B LDIR2B L2LIKB K2LIM1B M2LIM2B N2LIM1B N2LIM2B 

14 I P 1 1 33 1 49 

• 

MBCAVC 14 13 I I P M J X 1 97 1 33 1 49 1 33 1 49 

NSECS 
1 

LBL0CX2B LFACE2S LDIR2B L 21 I MB M2LIM1B K2LIK2B N2LIM1B N2LIM2B 

13 I H 97 1 33 1 49 

• 

»» Van* 2 

» 

•iSSVI 14 14 J J P P I X 1 1 1 97 1 49 1 97 1 49 

* RPHWALL TV ALL 

« 0.0 0.0 

• 

•SSVI 14 14 J J M M I X 33 33 1 97 1 49 1 97 1 49 

• RPHWALL TV ALL 

i 0.0 0.0 

• 

#SSVI 14 14 X K P P I J 1 1 17 81 1 33 17 81 1 33 

t RPHWALL TV ALL 

• o.o o.o 

•SSVI 14 14 K K M M I J 49 49 17 81 1 33 17 81 1 33 

• RPMVALL TWALL 

• 0.0 0.0 

KBCAVC 14 15 I I M P J X 97 1 1 33 1 49 1 33 1 49 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIKB K2LIM1B K2LIM2B N2LIH1B N2LIM2B 

IB I P 1 1 33 1 49 

* 

MBCAVC 15 14 I I P M J X 1 97 1 33 1 49 1 33 1 49 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIMB M2LIM1B M2LIK2B N2LIM1B N2LIM23 

14 I M 97 1 33 1 49 

« 

*• 3 lad* 2 

• 

i»SSVI 15 15 J J P P I X 1 1 1 97 1 49 1 97 1 49 

• RPMVALL TWALL 

• 12627.315 0.0 

•SSVI 15 15 J J M M I X 33 33 1 97 1 49 1 97 1 49 

» RPMVALL TWALL 

» 0.0 0.0 


••*»• HERE BE THE TEST - ADD HTPERSPACE TIP CLEARANCE 

•SSVI 15 15 X X P P I J 1 1 17 81 1 33 17 81 1 33 

« RPMVALL TWALL 

• 12627.315 0.0 

« 

•SSVI 15 15 X X M M I J 49 49 17 81 1 33 17 81 1 33 

» RPMVALL TWALL 

» 12627.315 0.0 


•SSVI 16 15 K K P P I J 1 1 17 81 1 29 17 81 1 29 

• RPMVALL TWALL 

• 12627.315 0.0 

•SSVI 15 15 X X M M I J 49 49 17 81 1 29 17 81 1 29 

« RPHWALL TWALL 

• 12627.316 0.0 

•PATCH 16 15 K X P M I J 1 49 17 81 29 33 17 81 29 33 

•PATCH 16 15 K K M P I J 49 1 17 81 29 33 17 81 29 33 


• HPT to LPT linkup 

•MBCAVC 15 16 I I M P J X 97 1 1 33 1 49 1 49 1 33 

•NSECS 

• 1 

•LBL0CX2B LFACE2B LDIR2B L2LIMB M2LIM1B K2LIM2B N2LIM1B N2LIM2B 

• 16 I PI 1 49 133 

* 

•MBCAVC 16 15 I I P M J X 1 97 1 49 1 33 1 33 1 49 

•NSECS 

• 1 

•LBL0CX2B LFACE2B LDIR2B L2LIKB K2LIM1B K2LIM2B N2LIM1B N2LIM2B 

*15 I M 97 1 33 1 49 

* 

* > This is th* fak* KPT *xit 

* 

•EXITC 15 15 I I M H J X 97 97 1 33 1 49 1 33 1 49 

•PEXIT 

•9.2424 


NASA/CR— 1998-206597 


123 



0 

0 > Shut off the HP turbine which is in the grid fll* 

0 (Essentially remove this from the LP analysis) 

0 

KILL 12 12 I I P P J X 1 1 1 33 1 49 

LSTART LEND 

1 97 

KILL 13 13 I I p p j K 1 1 1 33 1 49 

LSTART LEND 
1 97 

KILL 14 14 I I P P J K 1 1 1 33 1 49 

LSTART LEND 
1 97 

KILL IB 15 I I P P J X 1 1 1 33 1 49 

LSTART LEND 
1 97 


1 33 1 49 

1 33 1 49 

1 33 1 49 

1 33 1 49 


0 

0 > LPT INLET SPECIFICATION 

0 

INLETC 16 16IIPPJK 1 1 1 49 1 33 1 49 1 33 

PTOT TTOT 
9.2424 B. 01216 
0.6 4.46 


0 LOW PRESSURE TURBINE (LPT) 

0 Blocks 16-26 are alternating stator/rotor in 3d 
4 Blocks 26-65 are seal cavity geometries in 2d 

PATCH 16 16 X X P H I J 1 33 1 17 1 49 1 17 1 49 

PATCH 16 16 X X M P I J 33 1 1 17 1 49 1 17 1 49 

PATCH 16 16 X X P M I J 1 33 81 97 1 49 81 97 1 49 

PATCH 16 16 X X H P I J 33 1 81 97 1 49 81 97 1 49 

PATCH 17 17 K X P M I J 1 49 1 17 1 49 1 17 1 49 

PATCH 17 17 X X H P I J 49 1 1 17 1 49 1 17 1 49 

PATCH 17 17 X X P H I J 1 49 81 97 1 49 81 97 1 49 

PATCH 17 17 X X M P I J 49 1 61 97 1 49 81 97 1 49 

PATCH 18 18 X X P N I J 1 49 1 17 1 49 1 17 1 49 

PATCH 18 18 X X H P I J 49 1 1 17 1 49 1 17 1 49 

PATCH 18 18 X X P M I J 1 49 81 97 1 49 81 97 1 49 

PATCH 18 18 X K H P I J 49 1 81 97 1 49 81 97 1 49 

PATCH 19 19 X X P H I J 1 49 1 17 1 49 1 17 1 49 

PATCH 19 19 X X M P I J 49 1 1 17 1 49 1 17 1 49 

PATCH 19 19 X X P N I J 1 49 81 97 1 49 81 97 1 49 

PATCH 19 19 X X H P I J 49 1 81 97 1 49 81 97 1 49 

PATCH 20 20XKPMIJ 1 49 1 17 1 49 1 17 1 49 

PATCH 20 20KKMPIJ 49 1 1 17 1 49 1 17 1 49 

PATCH 20 20 X X P H I J 1 49 81 97 1 49 81 97 1 49 

PATCH 20 20 X X M P I J 49 1 81 97 1 49 81 97 1 49 

PATCH 21 21 K X P H I J 1 49 1 17 1 49 1 17 1 49 

PATCH 21 21 X X H P I J 49 1 1 17 1 49 1 17 1 49 

PATCH 21 21 X X P H I J 1 49 81 97 1 49 81 97 1 49 

PATCH 21 21 X X H P I J 49 1 81 97 1 49 81 97 1 49 

PATCH 22 22KKPMIJ 1 49 1 17 1 49 1 17 1 49 

PATCH 22 22KKMPIJ 49 1 1 17 1 49 1 17 1 49 

PATCH 22 22 X X P M I J 1 49 81 97 1 49 81 97 1 49 

PATCH 22 22 X X H P I J 49 1 81 97 1 49 81 97 1 49 

PATCH 23 23XXPMIJ 1 49 1 17 1 49 1 17 1 49 

PATCH 23 23XKMPIJ 49 1 1 17 1 49 1 17 1 49 

PATCH 23 23 X K P M I J 1 49 81 97 1 49 81 97 1 49 

PATCH 23 23 X X H P I J 49 1 81 97 1 49 81 97 1 49 

PATCH 24 24KKPMIJ 1 49 1 17 1 49 1 17 1 49 

PATCH 24 24 K X M P I J 49 1 1 17 1 49 1 17 1 49 

PATCH 24 24 X X P M I J 1 49 81 97 1 49 81 97 1 49 

PATCH 24 24 X X M P I J 49 1 81 97 1 49 81 97 1 49 

PATCH 25 25 X X P H I J 1 49 1 17 1 49 1 17 1 49 
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rotor 2 
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* Blocks 

42-49 compos* 

s*al cavity 

for 

rotor 3 
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* Blocks 

60-57 compos* 

seal cavity 

for 

rotor 4 
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f Blocks 

58-6! 

> compos* 

seal cavity 

for 

rotor 

5 

flow 

ord»r 68 

-59- 

-61- 

60-62 
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• Stator 

1 


















SSVI 

16 

16 

J 

J 

P 

P 

I 

X 

1 

1 

1 

97 

1 

33 

1 

97 

1 

33 


RPMVALL Ttf ALL 

0.0 0.0 

* Old tip «/ no s**ls 

#SSVI 16 16 J J M H I X 49 49 1 97 1 49 1 97 1 49 

* RPMVALL TV ALL 

* 0.0 0.0 


t Tip v/ s*al 

SSVI 16 16 J J M H I X 49 49 1 81 1 33 1 81 1 33 

RPMVALL TV ALL 

0.0 0.0 

SSVI 16 16 K X P P I I 1 1 17 81 1 49 17 81 1 49 

RPMVALL TV ALL 

0.0 0.0 

SSVI 16 16 K X M M I J 33 33 17 81 1 49 17 81 1 49 

RPMVALL TV ALL 

0.0 0.0 

MBCAYC 16 17 I I M P J K 97 1 1 49 1 33 1 49 1 49 

KSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIK3 M2LIM1B M2LIH2B N2LIM1B N2LIM2B 

17 I P 1 1 49 1 49 

MBCAVC 17 16 I I P M J X 1 97 1 49 1 49 1 49 1 33 

NSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIMB M2LIM1B K2LIM2B N2LIM1S N2LIH2B 

16 I M 97 1 49 1 33 


• Rotor 1 


SSVI 17 17 J J P P I X 1 1 1 97 1 49 1 97 1 49 

RPMVALL TV ALL 

3607.000 0.0 

SSVI 17 17 J J M M I X 49 49 1 97 1 49 1 97 1 49 

RPMVALL TV ALL 

3607.000 0.0 

SSVI 17 17 X X P P I J 1 1 17 81 1 49 17 81 1 49 

RPMVALL TV ALL 

3507.000 0.0 


SSVI 17 17 X X M M I J 49 49 17 81 1 49 17 81 1 49 

RPMVALL TV ALL 

3507.000 0.0 


MBCAVC 17 18 I I M P J X 97 1 1 49 1 49 1 49 1 49 

KSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIMB M2LIH1B H2LIM2B N2LIM1B N2LIM2B 
16 I P 1 1 49 1 49 


MBCAVC 18 17 I I P M J X 1 97 1 49 1 49 1 49 1 49 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIHB M2LIH1B M2LIM2B N2LIM1B N2LIM2B 
17 I M 97 1 49 1 49 
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# Stator 2 


SSVI 18 18 J J P P I K 1 1 1 97 1 49 1 97 

RPMVALL TV ALL 

0.0 0.0 

• Old tip v/ no sa&l 

•SSVI 18 18 J J M M I K 49 49 1 97 1 49 1 97 

• RPKWALL TV ALL 

• 0.0 0.0 


• Maw tip w/ aaal 

SSVI 18 18 J J M M I K 49 49 17 81 1 49 17 61 

RPMVALL TV ALL 

0.0 0.0 

SSVI 18 18 X K P P I J 1 1 17 81 1 49 17 81 

RPMVALL TVALL 

0.0 0.0 

SSVI 18 18 X X M M I J 49 49 17 81 1 49 17 81 

RPMVALL TVALL 

0.0 0.0 

MBCAVC 18 19 I I H P J X 97 1 1 49 1 49 1 49 

MSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIMB M2LIM1B M2LIM2B N2LIM1B N2LIH2B 

19 I P 1 1 49 1 49 

MBCAVC 19 18 I I P M J X 1 97 1 49 1 49 1 49 

NSECS 
1 

LBL0CX2B L FACE 23 LDIR2B L2LIKB M2LIM13 M2LIM2B V2LIM1B N2LIM2B 

18 I M 97 1 49 1 49 


V Rotor 2 


SSVI 19 19 J J P P I X 1 1 1 97 1 49 1 97 

RPMVALL TVALL 

3607.000 0.0 

SSVI 19 19 J J M M I K 49 49 1 97 1 49 1 97 

RPMVALL TVALL 

3607.000 0.0 

SSVI 19 19 X K P P I J 1 1 17 61 1 49 17 81 

RPMVALL TVALL 

3507.000 0.0 

SSVI 19 19 K X M M I J 49 49 17 81 1 49 17 81 

RPMVALL TVALL 

3607.000 0.0 

MBCAVC 19 20 I I M P J K 97 1 1 49 1 49 1 49 

NSECS 

1 

LBL0CK2B LFACE2B LDIR2B L2LIMB M2LIM1B M2LIK2B N2LIM1B N2LIM2B 

20 I P 1 1 49 1 49 

MBCAVC 20 19 I I P M J K 1 97 1 49 1 49 1 49 

NSECS 

1 

LBL0CX2B L FACE 23 LDIR2B L2LIKB K2LIM1B M2LIM2B N2LIM1B N2LIM2B 

19 I M 97 1 49 1 49 


• Stator 3 

SSVI 2020JJPPIX 1 1 1 97 1 49 1 97 1 

RPKWALL TVALL 

0.0 0.0 

• Old tip v/ no saal 

•SSVI 20 20 J J H M I X 49 49 1 97 1 49 1 97 

» RPMVALL TVALL 

• 0.0 0.0 

* Maw tip o/ aaal 

SSVI 20 20 J J M M I X 49 49 17 81 1 49 17 81 1 

RPMVALL TVALL 

0.0 0.0 


SSVI 20 20 X X P P I J 1 1 17 81 1 49 17 81 1 

RPMVALL TWALL 

0.0 0.0 


SSVI 20 20 X X M M I J 49 49 17 81 1 49 17 81 1 

RPMVALL TVALL 

0.0 0.0 


1 49 


1 49 


1 49 


1 49 


1 49 


1 49 


1 49 


49 


49 


49 


49 


49 


49 


49 


1 49 


49 


49 


49 
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t Rotor 3 













SSVI 21 

21 J 
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1 97 

1 49 

1 97 
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49 

RPMVALL 

TV ALL 












3807.000 


0.0 











SSVI 21 

21 J 

J H 
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I 

X 

49 

49 

1 97 

1 49 

1 97 
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49 

RPMVALL 

TV ALL 












3607.000 
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SSVI 21 

21 X 
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17 81 

1 49 

17 81 
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49 

RPMVALL 

TV ALL 












3607.000 
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SSVI 21 
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X N 
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J 

49 

49 

17 81 

1 49 
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49 

RPMVALL 

TV ALL 












3607.000 
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MBCAVC 21 
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1 49 

1 49 

1 49 

1 

49 

MSECS 
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49 
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M2LIM2B 

N2LIM1B 

N2LIH2B 



21 

I 
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97 


1 


49 

1 

49 




• Stator 4 


SSVI 22 22JJPPIX 1 1 1 97 1 49 1 97 1 49 

HP MV ALL TV ALL 

0.0 0.0 


• Old tip «/ no so&l 
iSSYI 22 22 J J M 

» R? MV ALL TV ALL 

» 0.0 0.0 

i Now tip w/ «*al 
SSVI 22 22 J J K 

RPHWALL TV ALL 

0.0 0.0 


SSVI 22 

22 X 

X P 

RPMVALL 

TV ALL 


0.0 

0 

.0 

SSVI 22 

22 X 

X M 

RPMVALL 

TV ALL 


0.0 

0 

.0 

MBCAVC 22 

23 I 

I M 

MSECS 



1 



LBL0CK2B LFACE2B 

LDIR2B 

23 

I 

P 

MBCAVC 23 

22 I 

I P 

MSECS 



1 



LBL0CK2B LFACE2B 

LDIR2B 

22 

I 

H 


t Rotor ' 

4 



SSVI 

23 

23 J 

J P 

RPMVALL 


TV ALL 


3507 

.000 


0.0 

SSVI 

23 

23 J 

J M 

RPMVALL 


TV ALL 



3607.000 0.0 


M I K 49 49 1 97 

M I K 49 49 17 B1 

P I J 1 1 17 81 

M I J 49 49 17 81 

P J K 97 1 1 49 

L2LIKB M2LIM1B H2LXK2B 
1 1 49 

K J K 1 97 1 49 

L2LIMB N2LXM1B M2LXH2B 
97 1 49 

P I K 1 1 1 97 

M I X 49 49 1 97 


1 49 1 97 1 49 

1 49 17 81 1 49 

1 49 17 81 1 49 

1 49 17 81 1 49 

1 49 1 49 1 49 

N2LIM1B N2LIH2B 
1 49 

1 49 1 49 1 49 

M2LIM1B N2LIM2B 
1 49 

1 49 1 97 1 49 

1 49 1 97 1 49 
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X P P 


I J 1 1 17 81 1 49 17 81 1 49 


SSVI 23 23 K 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 23 23 X X M H I J 49 49 17 81 1 49 17 81 

RPMVALL TV ALL 

3507.000 0.0 

KBCAVC 23 24 I I M P J K 97 1 1 49 1 49 1 49 

NSECS 
1 

LBL0CK2B LFACE23 LDIR2B L2LIKB M2LIM13 M2LIM2B N2LIH1B N2LIM2B 

24 I P 1 1 49 1 49 

HBCAVC 24 23 I I P M J K 1 97 1 49 1 49 1 49 

NSEGS 
1 

LBLQCX2B LFACE2B LDIR2B L2LIHB M2LIM1B M2LIH2B N2LIM1B N2LIH2B 

23 I M 97 1 49 1 49 

• Stator 5 

SSVI 24 24 J J P P I K 1 1 1 97 1 49 1 97 

RPMVALL TV ALL 

0.0 0.0 


§ Old tip v/ no saal 

•SSVI 24 24 J J M M I K 49 49 1 97 1 49 1 97 

• RPMVALL TWALL 

• 0.0 0.0 

f Nav tip v/ laal 

SSVI 24 24 J J M H I X 49 49 17 81 1 49 17 81 

RPMVALL TV ALL 

0.0 0.0 

SSVI 24 24 X K P P I J 1 1 17 81 1 49 17 81 

RPMVALL TV ALL 

0.0 0.0 

SSVI 24 24 K X H M I J 49 49 17 61 1 49 17 81 

RPMVALL TV ALL 

0.0 0.0 

KBCAVC 24 25 I I H P J X 97 1 1 49 1 49 1 49 

NSEGS 
1 

LBL0CK2B LFACE2B LDIR23 L2LIMB M2LIM1B M2LIM2B N2LIM1B N2LIH2B 

25 I P 1 1 49 1 49 

KBCAVC 25 24 I I P M J X 1 97 1 49 1 49 1 49 

NSEGS 
l 

LBL0CK2B LFACE23 LDIR2B L2LIK3 M2LIM1B M2LIM2B N2LIM1B N2LIM23 

24 I M 97 1 49 1 49 


• Rotor 5 

SSVI 25 26 J J P P I X 1 1 1 129 1 49 1 129 

RPMVALL TV ALL 

3507.000 0.0 

• Split to accapt tip aaal 

SSVI 26 25 J J M M I X 49 49 1 97 1 49 1 97 

RPMVALL TV ALL 

3507.000 0.0 


# Split to accapt tip aaal 

SSVI 25 25 J J M M I X 49 49 113 129 1 49 113 129 

RPMVALL TWALL 

3607.000 0.0 

SSVI 25 26 X X P P I J 1 1 17 81 1 49 17 81 

RPMVALL TWALL 

3607.000 0.0 

SSVI 25 26 X X M M I J 49 49 17 81 1 49 17 81 

RPMVALL TWALL 

3607.000 0.0 

• LPT to mixar linkup 

KBCAVC 25 66 I I M P J X 129 1 1 49 1 49 1 45 1 

NSEGS 

1 

LBL0CK2B LFACE2B LDIR2B L2LIKB M2LIM1B K2LIH2B N2LIH1B N2LIH2B 
66 I P 1 1 46 1 81 

HBCAVC 66 25 I I P M J X 1 129 1 46 1 81 1 49 1 

NSEGS 
1 


1 49 


1 49 


1 49 


1 49 


1 49 


1 49 


1 49 


49 


49 


49 


49 


49 


49 


49 


49 


81 


49 
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LBLDCX23 LFACE2B LDIR2B L2LIK3 M2LIM1B M2LIH2B N2LIM1B N2LIM2B 
25 I H 129 1 49 1 49 


• Rotor 1 **al cavity 2d/3d mixing plan* linkup and out*r boundari** 
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HBCAYC 16 

NSECS 

26 J 

J M 

P 

I 

X 

49 

1 

81 

97 

1 

33 

1 

17 
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97 
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33 



SSVI 27 
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K 

17 

17 
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17 

1 

2 
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17 
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0.0 

0.0 














SSVI 25 

28 J 

J P 

P 

I 

X 

1 

1 

1 

5 

1 

2 

1 

6 

1 

2 

RPMVALL 

TWALL 















3507.000 


0.0 














SSVI 28 

28 J 

J M 

M 

I 

X 

5 

5 

1 

5 

1 

2 

1 

5 

1 

2 

RPMVALL 

TWALL 
















0.0 0.0 


SSVI 29 29IIPPJK 1 1 113 1 2 113 1 2 

R? MV AIL TWALL 

3507.000 0.0 

SSVI 29 29JJPPIX 1 1 121 1 2 121 1 2 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 29 29JJMMIK17 17 121 1 2 121 1 2 

RPMVALL TV AIL 

0.0 0.0 


SSVI 3030IIPPJK 1 1 17 33 1 2 17 33 1 2 

RPMVALL TV ALL 

0.0 0.0 

SSVI 30 30 I I M M J K 21 21 129 1 2 1 29 1 2 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 30 30JJPPIX 1 1 121 1 2 121 1 2 

RPMVALL TWALL 

3507.000 0.0 

SSVI 30 30 J J M M I K 33 33 1 21 1 2 1 21 1 2 

RPMVALL TWALL 

0.0 0.0 

SSVI 3131JJPPIK 1 1 1 6 1 2 1 6 1 2 

RPMVALL TWALL 

3507.000 0.0 

SSVI 3131JJMHIK 6 5 1 5 1 2 1 5 1 2 

RPMVALL TWALL 

0.0 0.0 

SSVI 32 32 IIPPJK 1 1 113 1 2 113 1 2 

RPMVALL TWALL 

3507.000 0.0 

SSVI 32 32 JJPPIK 1 1 117 1 2 117 1 2 

RPMVALL TWALL 
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3507.000 


0.0 


SSVI 32 

32 J 

J 

M 

M 

I 

K 

17 

17 

1 

17 

1 

2 

1 

17 

1 

2 

RPMVALL 

TV ALL 
















0.0 

0 

.0 















SSVI 33 

33 I 

I 

P 

P 

J 

K 

1 

1 

1 

13 

1 

2 

1 

13 

1 

2 

RPMVALL 

TV ALL 
















3507.000 


0 . 

0 














SSVI 33 

33 I 

I 

H 

M 

J 

X 

17 

17 

1 

29 

1 

2 

1 

29 

1 

2 

RPMVALL 

TV ALL 
















0.0 

0 

.0 















SSVI 33 

33 J 

J 

H 

M 

I 

X 

29 

29 

1 

17 

1 

2 

1 

17 

1 

2 

RPMVALL 

TV ALL 
















0.0 

0 

.0 















HBCAVC 18 

33 J 

J 

H 

P 

I 

X 

49 

1 

1 

17 

1 

49 

1 

17 

1 

2 

KSEGS 

















LBL0CK2B LFACE2S LDIR2B 

L2LIMB 

M2LIM1B 

M2LIM2B 

N2LIM1B 

N2LIM2B 



33 

J 

P 



1 


1 


17 



1 


2 



KBCAVC 33 

18 J 

J 

P 

H 

I 

X 

1 

49 

1 

17 

1 

2 

1 

17 

1 

49 


NSEGS 

1 

LBL0CX2B LFACE2B LDIR2B L2LIM3 M2LIM1B M2LIM2B M2LIH1B N2LIM2B 
18 J M 49 1 17 1 49 


• Rotor 2 • cavity 2d/3d mixing plan* linkup and out*r boundari** 


SSVI 34 

34 I 

I P 

p 

J 

X 

1 

1 

1 

29 

1 

2 

1 

29 

1 

2 

RPMVALL 

TV ALL 















0.0 

0.0 














SSVI 34 

34 I 

I H 

M 

J 

X 

17 

17 

1 

9 

1 

2 

1 

9 

1 

2 

RPMVALL 

TV ALL 















3507.000 


0.0 














SSVI 34 

34 J 

J M 

M 

I 

X 

29 

29 

1 

17 

1 

2 

1 

17 

1 

2 

RPMVALL 

TV ALL 















0.0 

0.0 














MBCAVC 18 

34 J 

J M 

P 

I 

X 

49 

1 

81 

97 

1 

49 

1 

17 

1 

2 

MSECS 
















1 
















LBL0CX2B LFACE2B LDIR2B 

L2LIMB 

M2LIM13 

M2LIM2B 

N2LIM1B 

N2LIM2B 



34 

J 

P 


1 


1 


17 



1 


2 



MBCAVC 34 

18 J 

J P 

M 

I 

X 

1 

49 

1 

17 

1 

2 

81 

97 

1 

49 

MSECS 
















1 
















LBL0CX2B LFACE2B LDIR2B 

L2LIKB 

M2LIH1B 

M2LIM2B 

N2LIM1B 

N2LIM2B 



18 

J 

M 


49 


81 


97 



1 

49 



SSVI 36 

36 I 

I M 

M 

J 

X 

17 

17 

1 

17 

1 

2 

1 

17 

1 

2 

RPMVALL 

TV ALL 















3607.000 


0.0 














SSVI 36 

36 J 

J P 

P 

I 

X 

1 

1 

1 

17 

1 

2 

1 

17 

1 

2 

RPMVALL 

TV ALL 















3507.000 


0.0 














SSVI 36 

36 J 

J H 

M 

I 

X 

21 

21 

1 

17 

1 

2 

1 

17 

1 

2 

RPMVALL 

TV ALL 















0.0 

0 . 

.0 














SSVI 36 

36 J 

J p 

P 

I 

X 

1 

1 

1 

6 

1 

2 

1 

5 

1 

2 

RPMVALL 

TV ALL 















3507.000 


0.0 














SSVI 36 

36 J 

J M 

M 

I 

X 

6 

5 

1 

6 

1 

2 

1 

6 

1 

2 

RPMVALL 

TV ALL 















0.0 

0 . 

,0 














SSVI 37 

37 I 

I p 

P 

J 

X 

1 

1 

1 

13 

1 

2 

1 

13 

1 

2 

RPMVALL 

TV ALL 















3507.000 


0.0 














SSVI 37 

37 J 

J p 

P 

I 

X 

1 

1 

1 

21 

1 

2 

1 

21 

1 

2 

RPMVALL 

TV ALL 















3507.000 


0.0 














SSVI 37 

37 J 

J M 

M 

I 

X 

17 

17 

1 

21 

1 

2 

1 

21 

1 

2 

RPMVALL 

TV ALL 















0.0 

0 . 

.0 














SSVI 38 

38 I 

I p 

P 

J 

K 

1 

1 

17 

33 

1 

2 

17 

33 

1 

2 


RPMVALL TV ALL 
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0.0 


0.0 


SSVI 3838IIMMJK2121 129 1 2 1 29 1 2 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 38 38JJPPIX 1 1 121 1 2 121 1 2 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 38 38JJMHIX 33 33 1 21 1 2 1 21 1 2 

RPMVALL TV ALL 

0.0 0.0 

SSVI 39 39 JJPPIK1 1 1 5 121 5 1 2 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 39 39 JJMMIK 5 6 1 5 1 2 1 6 1 2 

RPMVALL TV ALL 

0.0 0.0 


SSVI 40 40IIPPJK 1 1 113 1 2 113 1 2 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 40 40JJPPIX 1 1 117 1 2 117 1 2 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 40 40JJMMIX17 17 117 1 2 117 1 2 

RPMVALL TV ALL 

0.0 0.0 

SSVI 4141IIPPJX 1 1 113 1 2 113 1 2 

RPMVALL TV ALL 

3607.000 0.0 

SSVI 41 41 I I M H J X 17 17 129 1 2 129 1 2 

RPMVALL TV ALL 

0.0 0.0 

SSVI 41 41 J JMMIX2929 1 17 1 2 1 17 1 2 

RPMVALL TV ALL 

0.0 0.0 

MBCAVC 20 41 JJMPIX 49 1 1 17 1 49 1 17 1 2 

NSEGS 

1 

LB LOCK 2B LFACE2B LDIR2B L2LIMB M2LIH1B H2LIM2B N2LIM1B K2LIM2B 

41 J P 1 1 17 12 

MBCAVC 41 20 J J P M I X 1 49 1 17 1 2 1 17 1 49 

KSECS 

1 

LBL0CX2B LFACE2B LDIR2B L2LIKB M2LIM1B K2LIM2B N2LIN1B K2LIM2B 

20 J M 49 1 17 1 49 


« Rotor 3 a*al cavity 2d/3d mixing plan* linkup and out*r boundari** 


SSVI 42 

42 I 

I P 

P 

J 

X 

1 

1 

1 

29 

1 

2 

1 29 

1 

2 

RPMVALL 

TV ALL 














0.0 

0.0 













SSVI 42 

42 I 

I M 

H 

J 

X 

17 

17 

1 

9 

1 

2 

1 9 

1 

2 

RPMVALL 

TV ALL 














3607.000 


0.0 













SSVI 42 

42 J 

J H 

M 

I 

X 

29 

29 

1 

17 

1 

2 

1 17 

1 

2 

RPMVALL 

TV ALL 














0.0 

( 

5.0 













MBCAVC 20 

42 J 

J M 

P 

I 

X 

49 

1 

81 

97 

1 

49 

1 17 

1 

2 

MSECS 















LBL0CK2B L FACETS 

LDIR2B 

L2LIHB 

M2LIN1B 

M2LIK2B 

N2LIM1B 

N2LIK23 



42 

J 

P 


1 


1 


17 



1 

2 



MBCAVC 42 

20 3 

J P 

M 

I 

X 

1 

49 

1 

17 

1 

2 

81 97 

1 

49 

KSECS 















L3L0CK2B LFACE2B 

LDIR2B 

L2LIHB 

H2LIM1B 

M2LIK2B 

K2LIM1B 

N2LIM2B 



20 

J 

M 


49 


81 


97 



1 

49 



SSVI 43 

43 I 

I M 

M 

J 

X 

17 

17 

1 

17 

1 

2 

1 17 

1 

2 

RPMVALL 

TV ALL 














3507.000 


0.0 













SSVI 43 

43 J 

J P 

P 

I 

X 

1 

1 

1 

17 

1 

2 

1 17 

1 

2 


RPMVALL TV ALL 

3507.000 0.0 
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SSVI 43 43 J J M M I K 21 21 1 17 1 2 1 17 1 2 

RPMWALL TWALL 

0.0 0.0 


SSVI 44 44JJPPIK 1 1 1 5 1 2 1 6 1 2 

RPMWALL TV ALL 

3607.000 0.0 


SSVI 4444J J H M I X 6 

RPMWALL TV ALL 

0.0 0.0 

SSVI 45 46 I I P P J X 1 

RPMWALL TWALL 

3607.000 0.0 

SSVI 4S45JJPPIK 1 

RPMWALL TWALL 

3607.000 0.0 


616121612 


1 1 13 1 2 1 13 1 2 


1 1 21 1 2 1 21 1 2 


SSVI 46 46 J J M M I K 17 17 1 21 1 2 1 21 1 2 

RPMWALL TWALL 

0.0 0.0 


SSVI 46461 IPPJK 1 1 17 33 1 2 17 33 1 2 

RPMWALL TWALL 

0.0 0.0 

SSVI 4646IIMMJX2121 129 1 2 1 29 1 2 

RPMWALL TWALL 

3607.000 0.0 

SSVI 46 46 JJPPIK 1 1 121 1 2 121 1 2 

RPMWALL TWALL 

3607.000 0.0 

SSVI 46 46 J J M M I X 33 33 1 21 1 2 1 21 1 2 

RPMWALL TWALL 

0.0 0.0 

SSVI 47 47 JJPPIX 1 1 1 6 1 2 1 6 1 2 

RPMWALL TWALL 

3607.000 0.0 

SSVI 47 47 JJMMIK 6 6 1 6 1 2 1 6 1 2 

RPMWALL TWALL 

0.0 0.0 


SSVI 48 48 IIPPJX 1 1 113 1 2 113 1 2 

RPMWALL TWALL 

3507.000 0.0 

SSVI 48 48 J J P P I X 1 1 1 17 1 2 1 17 1 2 

RPMWALL TWALL 

3607.000 0.0 

SSVI 46 48 J J M M I K 17 17 1 17 1 2 1 17 1 2 

RPMWALL TWALL 

0.0 0.0 


SSVI 49 49 I I P P J K 1 1 1 13 1 2 1 13 1 2 

RPMWALL TWALL 

3607.000 0.0 

SSVI 49 49 I I M M J X 17 17 1 29 1 2 1 29 1 2 

RPMWALL TWALL 

0.0 0.0 


SSVI 49 49 J J M H I X 29 29 1 17 1 2 1 17 1 2 

RPMWALL TWALL 

0.0 0.0 


MBCAVC 22 49 J J M P I K 49 1 1 17 1 49 1 17 1 2 

MSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIKB M2LIM1B K2LIM2B N2LIM1B N2LIM2B 

49 J P 1 1 17 12 

MBCAVC 49 22JJPMIX 149 117 1 2 117 149 

NSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIMB M2LIM1B M2LIM2B N2LIH1B N2LIM2B 

22 J M 49 1 17 1 49 


• Rotor 4 seal cavity 2d/3d mixing plan* linkup and outer boundaries 

SSVI 50 50IIPPJX 1 1 129 1 2 1 29 1 2 

RPMWALL TWALL 

0.0 0.0 
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SSVI 50 50 I I M M J K 17 17 1 9 1 2 1 9 1 2 

RPMVALL TV ALL 

3507.000 0.0 

SSVI 50 50JJMMIK29 29 1 17 1 2 1 17 1 2 

RPMVALL TVALL 

0.0 0.0 

KBCAVC 22 50 J J M P I K 49 1 81 97 1 49 1 17 1 2 

NSECS 

1 

LBL0CK2B LFACE2B LDIR2B L2LIMB M2LIM1B H2LIM2B N2LIM1B N2LIM2B 

50 J P 1 1 17 12 

KBCAVC 50 22JJPMIX 149 117 1 2B197 149 

NSECS 

1 

L3LDCX2B LFACE2B LDIR2B L2LIMB K2LIM1B H2LIM2B N2LIM1B N2LIH2B 

22 J N 49 81 97 1 49 

SSVI 61 51 I I M M J K 17 17 1 17 1 2 1 17 1 2 

RPMVALL TVALL 

3507.000 0.0 

SSVI 5151JJPPIK 1 1 117 1 2 117 1 2 

RPMVALL TWALL 

3507.000 0.0 

SSVI 51 51 J J M M I K 21 21 1 17 1 2 1 17 1 2 

RPMVALL TVALL 

0.0 0.0 

SSVI 62 52 JJPPIX 1 1 1 5 1 2 1 5 1 2 

RPMVALL TVALL 

3507.000 0.0 

SSVI 52 82 JJMHIK 55 1 5 1 2 1 5 1 2 

RPMVALL TVALL 

0.0 0.0 


SSVI 63 53 IIPPJX 1 1 113 1 2 113 1 2 

RPMVALL TVALL 

3507.000 0.0 

SSVI 63 53 JJPPIK 1 1 121 1 2 121 1 2 

RPMVALL TVALL 

3507.000 0.0 

SSVI 53 53 J J M M I K 17 17 1 21 1 2 1 21 1 2 

RPMVALL TVALL 

0.0 0.0 


SSVI 54 64IIPPJK 1 1 17 29 1 2 17 29 1 2 

RPMVALL TVALL 

0.0 0.0 

SSVI 64 54IIMHJX2121 1 25 1 2 1 25 1 2 

RPMVALL TVALL 

3507.000 0.0 

SSVI 6454JJPPIK 1 1 121 1 2 121 1 2 

RPMVALL TVALL 

3507.000 0.0 

SSVI 64 54JJHMIX2929 1 21 1 2 1 21 1 2 

RPMVALL TVALL 

0.0 0.0 


SSVI 65 56 JJPPIK 1 1 1 5 1 2 1 6 1 2 

RPMVALL TVALL 

3607.000 0.0 

SSVI 65 66 JJMMIK 5 5 1 5 1 2 1 6 1 2 

RPMVALL TVALL 

0.0 0.0 

SSVI 56 55IIPPJX 1 1 117 1 2 117 1 2 

RPMVALL TVALL 

3507.000 0.0 

SSVI 55 66JJPPIX 1 1 117 1 2 117 1 2 

RPMVALL TVALL 

3507.000 0.0 

SSVI 5656JJMMIK2121 117 1 2 117 1 2 

RPMVALL TVALL 

0.0 0.0 

SSVI 67 57 I I P P J X 1 1 1 13 1 2 1 13 1 2 

RPMVALL TVALL 

3507.000 0.0 
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M H J K 17 17 1 33 1 2 1 33 i 2 


SSVI 57 57 I I 

RPMWALL TWALL 

0.0 0.0 


SSVI 57 57 J J M M I X 33 33 1 17 1 2 1 17 

RPMWALL TV ALL 

0.0 0.0 


HBCAVC 24 57 J J M P I K 49 1 1 17 1 49 1 17 

NS ECS 
1 

LBL0CK2B LFACE2B LDIH2B L2LIMB M2LIM1B M2LIM2B N2LIM1B N2LIM2B 

67 J P 1 1 17 12 

HBCAVC 57 24 J J P M I K 1 49 1 17 1 2 1 17 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIKB M2LIM1B M2LIK2B N2LIM1B N2LIK2B 

24 J M 49 1 17 1 49 


• Rotor 5 i*bl cavity 2d/3d nixing plan* linkup and outar boundariaa 

SSVI 68 661 IPPJX 1 1 129 1 2 129 

RPMWALL TWALL 

0.0 0.0 


SSVI 58 6SIIHHJK17 17 1 9 1 2 1 9 1 

RPMWALL TWALL 

3507.000 0.0 


SSVI 68 68 JJMMIX29 29 1 17 1 2 1 17 1 

RPMWALL TWALL 

0.0 0.0 


HBCAVC 24 68 J J M P I K 49 1 81 97 1 49 1 17 1 

NSECS 
1 

L3L0CX23 LFACE2B LDIR2B L2LIKB K2LIM1B M2LIM2B N2LIM1B N2LIM2B 

68 J P 1 1 17 12 

HBCAVC 58 24 J J P M I X 1 49 1 17 1 2 81 97 1 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIMB M2LIM1B M2LIM23 N2LIM1B N2LIK2B 


24 

J 

N 


49 


81 


97 



1 

49 


SSVI 69 

59 I 

I M 

M 

J 

X 

17 

17 

1 

17 

1 

2 

1 17 

1 

RPMWALL 

TWALL 













3607.000 


0.0 












SSVI 59 

59 J 

J P 

P 

I 

X 

1 

1 

1 

17 

1 

2 

1 17 

1 

RPMWALL 

TWALL 













3607.000 


0.0 












SSVI 59 

69 J 

J M 

M 

I 

X 

21 

21 

1 

17 

1 

2 

1 17 

1 


RPMWALL TWALL 

0.0 0.0 


i Ramambar block* 45 and 46 ara awitchad from logical saquanca 
SSVI 60 60IIPPJX 1 1 117 1 2 117 1 

RPMWALL TWALL 

3507.000 0.0 


SSVI 6060JJPPIX 1 1 121 1 2 121 1 

RPMWALL TWALL 

3607.000 0.0 


SSVI 60 60 J J M M I X 21 21 1 21 1 2 1 21 1 

RPMWALL TWALL 

0.0 0.0 


SSVI 6161JJPPIX 1 1 1 6 1 2 1 5 1 

RPMWALL TWALL 

3607.000 0.0 


SSVI 61 61 J J M M I X 6 6 1 5 1 2 1 5 1 

RPMWALL TWALL 

0.0 0.0 


SSVI 62 62 I I M M J X 21 21 1 17 1 2 1 17 1 

RPMWALL TWALL 

3607.000 0.0 


SSVI 62 62 JJPPIX 1 1 121 1 2 121 1 

RPMWALL TWALL 

3607.000 0.0 


SSVI 62 62 J J M H I X 21 21 1 21 1 2 1 21 1 

RPMWALL TWALL 

0.0 0.0 


1 2 


2 


49 


2 


2 


2 


2 


49 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 
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SSVI 63 63 J J p p I K 

RPMVALL TWALL 

3607 .000 0.0 

SSVI 63 63J JMHIX 

RPHUALL TWALL 

0.0 0.0 

SSVI 64 64 I I p p j K 

RPMVALL TWILL 

3607.000 0.0 

SSVI 64 64JJPPIK 

RPMVALL TWALL 

3607.000 0.0 

SSVI 64 64 J J M M I K 

RPMVALL TWALL 

0.0 0.0 

SSVI 66 66 IIPPJ* 

RPMVALL TWALL 

3507.000 0.0 

SSVI 65 65IIHNJX 

RPMVALL TWALL 

0.0 0.0 

SSVI 66 65 J J M M I K 

RPMVALL TWALL 

0.0 0.0 

HBCAVC 25 65JJMPIX 

NSECS 
1 

LBL0CX2B LFACE2B LDIR2B L2LIKB 

66 J P 1 

HBCAVC 66 25 J J P M I K 

NSECS 
1 

LBL0CK2B LFACE2B LDIR2B L2LIK3 

25 J M 49 


111612161 


661512161 


1 1 1 13 1 2 1 13 1 


1 1 1 17 1 2 1 17 1 


17 17 1 17 1 2 1 17 1 


1119 12 19 1 


17 17 1 26 1 2 1 26 1 


25 26 1 17 1 2 1 17 1 


49 1 97 113 1 49 1 17 


M2LIM1B K2LIM2B N2LIM1B N2LIK2B 
1 17 12 

1 49 1 17 1 2 97 113 


M2LIM1B H2LIM23 N2LIM1B N2LIM2B 
97 113 1 49 


• MIXER AND EXIT CRIDDXNC 
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1 2 
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1 

2 

PATCH 

66 72 
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M 
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72 66 
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J 

PATCH 

66 66 
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66 66 
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P 

H 

I 
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67 68 
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J 

H 
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I 
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68 67 
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J 
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PATCH 

67 67 

X 

X 

M 
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PATCH 
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X 

X 

P 

M 

I 

PATCH 

68 72 

I 

I 

P 

M 

J 

PATCH 

72 68 

I 

I 

M 

P 

J 

PATCH 

68 69 

J 

J 

M 

P 

I 

PATCH 

69 68 

J 

J 

P 

M 

I 

PATCH 

68 68 

X 

X 

H 

P 

I 

PATCH 

68 68 

K 

X 

P 

M 

I 

PATCH 

69 74 

I 

I 

P 

M 

J 

PATCH 

74 69 

I 

I 

M 

P 

J 


L L M 
12 1 
L L L 
I I I 
M M M 


M N N 

1 1 1 

L L L 

1 I I 

M M M 

2 12 


M H N 
2 2 2 
L L L 
I I I 
M K H 
1 2 1 


K 

K 

J 

J 

K 

K 

J 

J 

K 

K 

X 

X 

J 

J 

X 

K 


65 

1 

1 

45 

1 

65 

1 

45 

81 

1 

1 

66 

1 

81 

1 

65 

9 

1 

1 

81 

1 

9 

1 

81 

81 

1 

1 

81 

1 

81 

1 

81 

1 

49 

1 

49 

49 

1 

1 

49 

49 

1 

1 

81 

1 

49 

1 

81 

81 

1 

1 

81 

1 

81 

1 

81 

1 

49 

1 

49 

49 

1 

1 

49 


1 

81 

1 

46 

1 

81 

1 
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0 . OOOOOOE+OO O.OOOOOOE+OO 


SSVI 73 73 J J M H I x 49 49 1 65 1 81 1 65 1 81 

RPKVALL TVALL 

O.OOOOOOE+OO O.OOOOOOE+OO 

SSVI 74 74 J J N M I X 49 49 1 49 1 81 1 49 1 81 

RPMVALL TV ALL 

O.OOOOOOE+OO O.OOOOOOE+OO 

SSVI 74 74 X X P P I J 1 1 33 49 45 49 33 49 46 49 

RPMVALL TW ALL 
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Appendix B 

SEARCH Program Source Code 


The source code for the SEARCH mesh generation utility program developed during 
this study is printd below for reference. 


program saarch 

paramatar(mxdim-U00,idia-150, jdtm«100,kdim-49, nblks-50) 

raal t(-5:mxdim) ,s(-5:axdim) ,v(0 :mxdim,0;mxdim) , 

. x(0;mxdim,0:axdim) , y(0:mxdi»,0 :mxdim) , 
*(0:mxdim,0:mxdim) ,u(0:l),v(0:l), 
sla(axdim) ,»ta(axdim) ,tla(mxdia) ,tta(mxdim) , 
xt« (mxdia) , yta {mxdia) ,xta(mxdim) .rta(mxdia) , 
xl« (mxdia) ,yla (axdia) ,zla (axdla) ,rla (mxdia) 

raal csla (axdim) , cat a (mxdia) .cxla(axdim) ,cxta (axdim) , 
cyla (mxdia) ,cyta (axdia) ,czla(axdim) ,c 2 ta(mxdim) 

raal xn(idim, jdim,kdim) , 
yn(idim, jdim.kdim) , 
xnfidia, jdim.kdim) , 
x2d(idia» jdia) , r2d(idia, jdim) , 
that ( idim .jdim , kdim) , rtha t (kdim) 

intagar il(nblks) ,jl(nblks) ,kl(nblki) 


character *60 filnm 


common/ surf data/ 
c oamon/par ami ia/ 
c osmon/nurb s 1 z a / 
common/ax imash/ 
common/ful Lma *h/ 
c omaon/aa s hi ia/ 
c oamon/ xyzadga a / 

common/ < t _a dga a / 
c oamon/ xrpo i n t / 
common/ttval/ 
common/* sval/ 
common/ if lags/ 


>diV|X 

u,v 

kl,k2 

x2d,r2d 
xn , yn , xn 
il,jl,kl,mg 

xla.xta ,yla,yta ,zla,zta ,rla,rta, 
cxla , cxta.cyla, cyta,czla,czta 
*la , sta.tla ,tta ,c*la ,c«ta 
xx.rr 
tt 

• ft 

isaarch, Ida bug, lin, lout , lgrid 


♦** Xnititializa I/O, flag* and constants ♦♦♦ 
a 

lin - 6 

lout ■ 6 

lgrid - 12 

idabug - 0 

pi ■ 4 . 0*atan(l .0) 


*** Aaad in ICES axcarpt fila *a* 
a 

writ* (lout ,*) 'Entar NASA ICES filanama:' 
raad(lin, * (a) ’) filna 
vrita (lout ,*) 'Filanama: filna 

call raadiga s (filna) 

* 

*** Find radial distributions of Xmin and Umax *♦* 

• 

writa(lout,*) ’Entar nuabar of »plina point*: {- to dabug NURBS) * 
raad(lin,*) nsp 

rrita{lout , *) ’Numbar of splina points:', nsp 
if (nsp. It. 0) idabug » 1 
nsp ■ abs(nsp) 
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call lataf ind(nsp) 


writadout,*) 'Maridional landing and trailing *dg** vrittan. ' 
writadout,*) ’Continua? (1-y)’ 
raadClin,*) inns 
if (ians.na.l) atop 

*** Calculata aplina coafficianta *** 
call a pc oaf (tla, ala ,c ala ,nap) 
call apcoaf (tta.ata.cata.nap) 

call spcoaf (rla,xla,cxla,nsp) 
call spcoaf (rta.xta.cxta, nap) 

call spcoaf (rla.yla ,cyl* ,nsp) 
call apcoaf (rta.yta, cyta, nap) 

call apcoaf (rla.sla.czla , nap) 
call apcoaf (rta.xta.cxta, nap) 

* 

#** Intarrogata NURBS aurfaca uaing aaarch routinaa (idabug - 1) *** 

* 

if (idabug. na . 0) than 

5 writadout ,*) 'Entar a,t valua: (a - 999 to quit)* 

raadClin,*) aa.tt 
if (aa.na.999.) than 

call gataurf xrt (aa ,tt ,xa ,ra, ta) 
vrita(lout ,*) 'xrt -> ' ,xa ,ra ,ts*180 ./pi 
go to 5 
and if 

7 writadout,*) ’Entar x,r valua to match: (x-999 to quit) 1 

raadClin,*) xx.rr 
if (xx.na.999) than 
iaurf ■ 1 

aa - u(0)+0.25*(u(l)-u(0}) 
tt - v(0)+0.50*(v(l)-v(0)) 

writadout,*) ’Uaing altamatlng aacant aaarch...’ 
iaaarch ” 1 

call f indat (nap, iaurf ) 

call gataurf xrt (a* ,tt ,xa ,ra , ta) 

writadout,*) '«,t *> ’, aa.tt 

writa (lout ,*) *xrt -> ’ ,xs ,rs ,ts*l&0. /pi 

writ*(lout ,*) 'Uaing undar-ralaxad gradiant aaarch. . . ’ 

iaaarch « 2 

call f indat (nap, iaurf ) 

call gataurf xrt ( aa, tt ,xa ,ra,ta) 

writa (lout ,*) ’i,t -> '.aa.tt 

writa (lout , *) ’xrt ■> * ,xa ,ra ,ta*180. /pi 

go to 7 

alaa 

atop 
and if 
and if 

• 

♦** Raad in 2-D PLQT3D math for maridional point diatribution *** 

♦ 

writadout ,*) 'Entar 2-D axis yin. PL0T3D binary maah fllanam*;’ 

raaddin, ’ (a30) 1 ) filnm 

writa (lout,*) 'Filanam* : filnm 

call raadaxiaaah(f ilna) 

writadout,*) ’Entar landing adga indica:’ 
raadClin,*) ila 

writa (lout,*) ’Landing adga: *,ila 

writadout,*) ’Entar trailing adga indica:' 
raad(lin,*) ita 

writadout,*) ’Trailing adga: *,ita 

writadout,*) 'Entar numbar of bladaa in whaal: * 
raadClin,*) nbladas 

writadout,*) 'Blada count: * , nbladas 

pitch • 2-*pi/nbladaa 

writa(lout,*) ’Pitch: ’ ,pitch*180. /pi 

writadout,*) 'Entar numbar of pointa across passago:’ 
raadClin,*) kl(l) 
writadout,*) ’Kpts: ’ ,kl(l) 

writadout,*) ’Entar naar blada spacing for first point’ 
writadout,*) ’off tha bind* aurfaca:' 
raadClin,*) drthl 

writadout.*) 'Naar-wall spacing: '.drthl 

drth2-drthl 

a 

*♦* assign landing and trailing adga valuas from splinas *** 

* 

1 » 1 
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do 60 j-l.jl(l) 

ytosnp - spaval(r2d(ila,j) ,rla,yla,cyla,nsp) 
ztamp - spaval(r2d(ila,j) ,rla,zla,cyla,nsp) 
that (ils, j , 1) « atan2(ztamp, ytamp) 
that(ila,j ,kl(l)) - that (ila ,j , 1) 

ytamp - spaval(r2d(ita,j) ,rta,yta,cyta,n#p) 
xtamp - spaval(r2d(ita,j) ,rta,zta,cyta,nsp) 
that ( ita , j , 1) * atan2(ztsmp, ytamp) 
that (Ha, j,kl(l)) - that (ita, j ,1) 

60 continua 

♦ 

*** find thata valuas for intarior x,r pairs ♦** 
a 

writa Clout,*) 'Entsr saarch aathod: ' 
writ* (lout, •) * 1. Alternating sacant saarch* 

writs (lout , *) ’ 2. Undar-ralaxad gradisnt saarch* 

raad(lin, a) isaarch 

writa(lout ,*) * Saarch mathod; ’,isaarch 
do 100 1 « 1 , mg 
kl(l) - kl(l) 
do 100 isurf - 1,2 
writa(lout.a) 'Surfaca isurf 
do 100 i-ila+l,lta-l 

writs (lout , *) 1-ila+l,* of '.ita-ila+1 
idsbug — 0 
do 100 j-l,jl(l) 
xx - x2d ( i , j ) 
rr - r2d ( i , J ) 
call f indst (nsp, isurf ) 
call gatsurf xrt (as ,tt,xs,rs,ts) 
if (isurf .aq. 1) thst(i, j,l)-ts 
if (isurf. aq. 2) that(i, j ,kl(l))-ts 
100 continua 

♦ 

aa* Sat thatas upstraam and dovnstraam of blada aaa 


if (ila.na.l) than 
** upstraam ♦* 

do 200 i-l.ila-l 
do 200 j-l,jl(l) 

that (i , j , 1) - that (ila, j ,1) 
that (i, j ,kl(l)) - that(ila,j ,kl(l)) 
200 continua 

and if 

if (ita .na .il(l)) than 
aa dowcstrsam aa 

do 210 i-ita+1 , il(l) 
do 210 j-l,jl(l) 

that(i, j , 1) - that (ita, j ,1) 
that (i, j ,kl (1» - that (ita ,j ,kl(l)) 
210 continua 

and if 


♦** Distrlbuta points across passaga *** 


*** switch surfacas 142? **• 
ia ■ lot (0.5*(ila+ita) ) 
jm - int(0.5*jl(l)) 

thatl • that(im, jm, 1) 
that 2 ■ that (im, jm,kl(l))+pitch 
dal that - abs (thatl - that2) 
writadout,*) 'dalthat - * ,dalthat*180. /pi 
writa(lout ,a) 'pitch • ’ ,pitch*180 . /pi 
if (dalthat .gt .pitch) than 
*** swap surfacas 

writs (lout, *) ’swapping thata surfacas...' 

1-1 

do 360 i-l.il(l) 
do 360 j«l,jl(l) 
tamp-that (i, j ,1) 

that(i,j,l)-that(i,j,kl(l)) 

that (i, j ,kl(l) )-tsmp 
350 continua 
and if 

1-1 

do 400 i«l,il(l) 
do 400 j-l.jl(l) 

*** shift surfaca by pitch 

that(i , j ,kl(l>) - that (i,j ,kl(l))+pitch 

kasa - 0 

iarr - 0 
kpts - kl(l) 

rtl-r2d(i, j )*that (i, j ,1) 
rt2*r2d(i, j) *that (i , j ,kl(l)) 
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eqspc - (rt2 - rtl) / real(kpts-l) 
if (abs(eqspc).gt.drthl) then 

us* Vinokur clustered spacing *** 

call VINOKUR (rthet ,kpts ,rtl ,rt 2 ,drthl ,drth2,KASE, IERR) 
if (ierr.ne.O) th*n 

writedout,*) 'Error in VINOKUR routin*. . .stopping’ 
stop 
•nd if 

•Is* 

••* us* *qual rth*t spacing *** 

do 390 k»2,kl(l)-l 

rthet (k)-rtl+eqspc*real(k-l) 

390 continu* 

•nd if 

do 400 k«2.kl(l)-l 

thet(i , J,k)-rthet(k)/r2d(i , j) 

400 cont inu* 

* 

*** Anglo upstream and downstream grid points to natch blad* angl* **+ 

write(lout ,*) 'Enter circumferential angl* for upstream: (d*g) ' 

r* ad(lin,*) ucang 

vrit*(lout ,*) 'Upstream: », ucang 

if (ucang. n*.0) th*n 
ucang— ucang*pi/ 180 . 

* 

**** shift extension block *** 

* 

1-1 

do 800 i-ile-1.1,-1 
do 800 j-ljl(l) 
do 800 k-l.kl(l) 

* tangential location 

d* " xn(i, j ,k)-xn(ile, j ,k) 
dt - dx/yn(i,l.k)*atan (ucang) 

xn(i,j,k)-zn(ile,j,k)+dt 

300 continu* 

•nd if 

vrit* (lout ,*) 'Enter circumferential angle for downstream; (deg)’ 
re*d(lin,*) dcang 

wrlte(lout,*) 'Downstream: dcang 

if (dcang. ne.0) then 

dcang-dcang*pi/180, 

* 

**** shift extension block *** 

* 

1-1 

do 810 i-lte+1 ,11(1) 
do 810 j-l,jl(l) 
do 810 k-l.kl(l) 

* tangential location 

dx * xn(i , j ,k)-xn(ite, j ,k) 
dt - dx/yn(i,l,k)*at an (dcang) 
xn(i,j ,k)-xn(ite , j ,k)+dt 
810 continue 
end if 


*** Convert back to cartesian coordinates *** 

e 

writedout,*) 'Converting to cartesian coordinates...’ 
1-1 

writ*(lout,*) 'Block: * ,1 
do 700 i-l t ll(l) 
do 700 j-l.jl(l) 
do 700 k-1 ,kl(l) 
xn(i,J,k)-x2d(i,j) 

yn(i » j ,k)-r2d(i , j ) *cos (thet (l , j ,k) ) 

*n(i , j ,k)-r2d(i, j)*sin(th*t(i,J,k)) 

700 continu* 


*** Output PL0T3D file *•* 
e 

vrit«(lout,*) * — OUTPUT ’ 

c writ*(lout,*) ’Enter 3-D PLQT3D binary filename:’ 

c read(lin, ' (aBO) ’ ) filna 
f ilna - 'blade. mesh' 
writedout,*) 'Filename: ’.filnm 

call qdopendgrid.f ilnm, je) 
call qdputi (lgrld.mg, Je) 
writedout,*) ’Number of grids: ’ ,mg 

do 899 1 - 1, mg 

call qdputi (lgrid.il(l) , je) 
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call qdputi(lgrid, jl(l) . j*) 

call qdputi (Igrid, kl(l) , j*) 

899 continue 
1-1 

vrite(lout,*) 'Block: ’ ,1,U(1) , jl(l) ,kl (1) 
do 911 k - 1, kl(l) 
do 911 j - 1, jl(l) 
length - i 1 ( 1 > 

call qdpuea(lgrid,xn(l, j ,k) , length, je) 

911 continue 

do 912 k - 1. kl(l) 

do 912 j - 1, jl(l) 
length - il(l) 

call qdpueaflgrid ,yn (1 , j ,k) , length, Je) 

912 continue 

do 913 k - 1, kl(l) 

do 913 j - 1, jl(l) 
length - 11(1) 

call qdpuea(lgrid,zn(l ,j ,k) .length, je) 

913 continue 

call qdcloe (Igrid, Je) 

• top 
end 


SUBROUTINES 


**** 


•ubroutine readaxime*h (filnm) 
parameter (idim-150,j dim- 100, nblk*-50) 

real x2d(idim, jdim) ,r2d(idim, jdim) 
integer il (nblk«) , jl(nblk») ,kl (nblk*) 

character*80 filnm 

co®mon/axime*h/ x2d,r2d 
common /me* hi ia/ il,jl,kl,mg 

common/ if lag*/ isearch, i debug, lin , lout , Igrid 

call qdopenf Igrid, filnm, je) 

call qdgeti( Igrid, mg, je) 

write (lout ,*) 'Number of grid*: ’ ,mg 

if (ag.gt.l) then 

write(lout,*) 'Warning!! Multiple grid* in 2-D file...’ 
vrite (lout ,*) ’Only u*ing grid block il . ’ 
end if 

do 9 1 - 1 , mg 

call qdget i( Igrid, i 1(1) , je) 
call qdgeti ( Igrid, j 1(1) , je) 
call qdgeti (Igrid, kl(l) , je) 

9 continue 

♦** only single me*h now (1-1) ***• 

1-1 
mg - 1 

write(lout , +) ’Block: * ,1,11(1) p jl(l) 
do ill j - l, jid) 
length - il(l) 

call qdgeea(lgrid,x2d(l, j), length, je) 

111 continue 

do 112 j - 1, jl(l) 
length • 11(1) 

call qdgeea(lgrid,r2d(l, j) , length, je) 

112 continue 

call qdclo* (Igrid, je) 

return 

end 


•ubroutine f ind*t (n*p, i*urf) 
parameter (mxdi*-l 100) 

real t(-6:mxdim) .•(-6:mxdim) ,w(0:mxdim,0:mxdim) , 
x(0:mxdim,0:ttxdim) ,y(0:mxdim,0:mxdim) , 
x(0 : mxdim, 0; mxdim) ,u(0: 1) ,v(0: 1) , 

•le (mxdia) ,*te (mxdim) ,tle (mxdim) , t te (mxdia ) t 
cele(mxdim) ,c*te(mxdis) 
real p(2),xi(2,2) 


common/ surf data/ *,t,w,x,y,z 
coannon/paramlim/ u , v 
comaon/nurbsize/ kl,k2 

common/ st .edges/ *le, *te,tle,tte .csle.cste 
common/ xrpoint/ xx,rr 
common/* «val/ ** 

common/ ttval/ tt 

common/ if lag*/ isearch, idebug, lin, lout , Igrid 
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• •* initial gusss point and sat bounds for (a,t) 

* 

tt • v(0)+0.6*(v(l)-v(0) ) 

tain - v(0) 
taax ■ v(l) 

sain * spaval(tt ,tla ,sla ,csla ,nsp) 
saax » §psval(tt ,tta , sta .csta.nsp) 

if (sain. gt .saax) than 
stamp “ sain 
sain “ sou 
saax ■ stsap 
•nd if 

if (isurf.aq.2) thsn 
*** swap sain and saax *** 
stsap • sain 
sain ■ saax 
saax • stsap+u(l) 
and if 

ss “ 0 . 5* (sain+smax) 

if (issarch.sq. 1) thsn 


*♦* Altsrnating sscant ssarch (inaw - 0) *** 


ftol « 0.0005 
tol - 0.1 

ipass ■ 1 

5 if Ctol.lt. 0.0001) tol - 0.0001 

♦ 

♦*« SVEEP T at CONST S ♦•* 

* 

15 told - 0 . 95*tt 

if {told. It. 0.1) thsn 
tt - 0.1 
told - tt+0.5 
and if 

dslt - tt - told 

if (Idsbug. ns .0) writa(lout ,*) *tt .told, dslt » ’ ,tt .told, dslt 
call gstsurfxrt (ss .told, xs ,rs ,ts) 

c writa(3, *) xs.rs 

frold- rr - rs 

itsr - 1 

if (idsbug. ns .0) writs (lout ,*) ’rr.frold - ' .rr.frold 

20 call gstsurfxrt (ss ,tt ,xs ,rs ,ts) 
c writs(3,*) xs,rs 

frnsv • rr - rs 

if ( idsbug . ns . 0) writa(lout ,*) 'ss.tt, it, frnsw - ’ ,i»,tt p rr,fniM 
dslt ** -fmsw / (frnsw - frold)*dslt 
ttsap ■ tt + dslt 

If (idsbug. ns. 0) writs(lout ,*) ’ttsap, dslt, tain, taax 
t t sap, dslt .tain, taax 

if (ttsaqj.lt .tain) thsn 
dslt - tain - tt 
ttsap “ tain 
and if 

if (ttsap. gt .taax) thsn 
dslt * taax - tt 
ttsap ■ taax 
and if 

tt ■ ttsap 

if (idsbug. ns. 0) writs (lout ,♦) ’itsr, tt, dslt - itsr ,tt , dslt 

ss* try until tolsrancs is ast or 10 tiass sxcssdsd *** 
if (abs(dslt) .gt.tol) thsn 
frold ■ frnsw 
itsr ■ itsr + 1 
if (itsr.lt.il) go to 20 
and if 


ss* SWEEP S at CONST T *** 
* 

7 itsr - 1 


NASA/CR— 1998-206597 


144 



•min - spaval (tt , tla , si* , cal* , nsp) 
smax ” spaval (tt,tta,sta,csta,nsp) 

if (idabug.na .0) writa(lout,*) ’laS, taS: sain, smax 

if (smin. gt . smax) than 
stamp - sain 
smin * smax 
ixax - stamp 

if (idabug.na.O) writa (lout ,*) 'swappad: min, max; * , smin, smax 
and if 

if (isurf.aq.2) than 
•** swap smin and smax •** 
stamp ■ smin 
smin ■ smax 
•max * staap+u(l) 
and if 

sold ■ 0.95+ss 

if (sold. la . smin) sold ■ smin +0.1 
if (ss.aq.smax) sold * smax - 0.1 
dais ■ ss - sold 

if (idabug.na.O) writa (lout , ♦) ’ss, sold, dais - ’ , ss , sold, dais 

call gatsurfxrt(sold,tt,xs,rs,ts) 
c writa(3,+) xs,rs 

fxold* xx - xs 

if (idabug.no . 0) vrita(lout ,♦) ’smin, smax * smin, smax 

10 call gatsurfxrt (ss, tt ,xs ,rs ,ts) 

c writa(3,*) xs,rs 

fxnaw ■ xx - xs 

dais » -fxnaw / (fxnaw - fxold)+dals 

stamp ■ ss + dais 

if (stamp. It . smin) than 
dais ■ smin - ss 
stamp - smin 
and if 

if (stamp. gt .smax) than 

dais * smax - ss 
stamp - smax 
and if 

ss * stamp 

*** try until tolaranca is mat or 10 timas axcaadad *aa 
if (abs(dals) .gt . tol) than 
fxold ■ fxnaw 
itar ■ itar + 1 
if (itar.lt.il) go to 10 
and if 

* 

+++ calculate distanca »++ 

a 

25 call gatsurfxrt (ss, tt ,xs,rs,ts) 
c writa(3,+) xs,rs 

dist » sqrt( (xx-xs)+»2 + (rr-rs)a*2) 

ipass ” ipass 4 1 

if (idabug.na.O) than 

writa (lout , •) ’PASS i’, ipass 
writa(lout,*) ’DIST.s.t - ’,dist,ss,tt 
writa (lout , +) ’smin, smax* smin, smax 

writa (lout ,a) ’ dx , dr - ’ ,xx-xs ,rr-rs 

and if 

if (ipass .gt .7) than 

writa (lout ,*) 'giving up ratry ipass 

writa(lout ,*) 'DIST.s.t - ',dist,ss,tt 
writa (lout , +) ’ smin, smax* * , smin, smax 

writa (lout , +) *dx,dr * ’ , xx-xs , rr-rs 
rat urn 
and if 

if (dist .It . f tol) than 
*+* solution found aaa 
raturn 

also 

**+ try to gat closar (tight an tolaranca) aaa 
tol - tol/10. 

if (i da bug . na . 0) writa (lout , *) ' > 

if (idabug.na.O) writa(lout,+) ’NEW tol - *,tol 
go to 5 
and if 

also 
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* sispl* derivative ••arch 


it*r - 1 

call g«t*urfxrt(*»,tt ,x« ,r* , t*) 

23 call d*riv*(«*,tt ,dxd*,drd*,dxdt .drdt ,n*p,i*urf ) 

dx - xx - x» 
dr - rr - r* 

If (dxd*.D«.0) then 
delil - dx / dxds 

•1m 

d«ltl - 0.0 
•nd If 

if (drd».n«.0) than 
d« 1*2 - dr / drdt 

•Is* 

d«la2 - 0.0 
•nd If 

d* “ d«lsl + d«l»2 

if (dxdt.n*.0) than 
d*ltl - dx / dxdt 

• ls« 

d*ltl - 0.0 
•nd If 

if (drdt .n« .0) th«n 
d«lt2 - dr / drdt 

•li« 

dtU2 - 0.0 
•nd if 

dt - d«ltl + d«lt2 

•♦* und«r-r*lax valu* • ** 
omega - 0.1 
sn«v ■ •• + d**om*ga 
tn«v - tt ♦ dt*om*ga 

if (tn«v. It . v(0) ) tn«v ■ v(0) 
if (tn*w.gt.v(l)) tn*w • v(l) - 0.0001 

tt - tn*v 

tain - *p*val(tt , tl«,»l* ,ctl« ,n*p) 

■ «p«val(tt ,tt«.st* ,c»t*,ntp) 

if (amin.gt . smax) than 
•t*mp * tain 
•min - ssax 
•max • tt*np 
•nd if 

if (isurf.*q.2) then 
*** swap tain and «sax **• 
t temp « tain 
•min “ tmax 
•max ■ *t#mp+u(l) 

•nd if 

if (sn«w.l«.*min) »n«w - »min ♦ 0.01 
if (sn«w.g*. smax) anew * »max * 0.01 
•• - ■□•v 

call g*t*urfxri(»*,tt ,x* ,rs,t«) 
ditt - •qrt((xx-x*)**2 + (rr-r»)**2) 

if (it«r.g«. 100) then 

writ*(lout,*) 'D*riv ••arch FAILED!!!' 
writ* (lout,*) ’ iter,di«t ■',it*r,di*t 
return 
•nd if 

if (dlst.lt.ftol) then 
•** solution found **♦ 
return 
• lM 

lt*r - it*r * 1 
go to 23 
•nd if 


•nd if 
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• ad 


subroutine derive (ss ,tt ,dxds ,drds ,dxdt ,drdt , nip, isurf) 
parameter (axdia*l 100) 
r*al u(0: 1) ,v(0: 1) , 

sle(mxdia) ,ste(axdim) ,tle(axdim) ,tte(axdim) , 

. csl« (nut dim) .cite (mxdim) 

comaon/paraalia/ u,v 

common/ st.edges/ sle,ste,tle,tte,csle,cste 

c if (isurf.eq.2) write(lout.e) * derive — ’ 

c if (isurf.eq.2) write (lout ,*) 'ss,tt - ’,ss,tt 

sain ■ speval(tt ,tle,sle ,csle,nsp) 
saajt - spevaKtt.tte.ste.cste.nsp) 

c if (isurf.eq.2) write (lout,*) ’u,v - ' ,u(0) ,u(l) ,v(0) , v(l) 

if (sain.gt . saax) then 
•temp “ sain 
sain ■ saax 
saax ■ stemp 
end if 

if (isurf.eq.2) then 
*** swap sain and saax *** 
stemp - sain 
sain “ saax 
saax ■ stemp+u(l) 
end if 

c if (isurf.eq.2) vrite(lout ,*) 'sain, saax ■ sain, saax 

deism * 0.001 
delta “ 0.001 
delsp • 0.001 
deltp - 0.001 

•1 " ss-delsa 

if (sl.le.smin) then 
sl“sain 

deism * ss - si 

end if 

t2 - tt+deltp 
if (t2.ge . v(l) ) then 
t2 - v(l) 
deltp - t2 - tt 
end if 

•3 ■ ss + delsp 

if (s3.ge.saax) then 
s3 B saax 

delsp - s3 - ss 

end if 

t4 ■ tt-delta 
if (t4.1e.v(0)) then 
t4 - v(0) 
deltp ■ tt - t4 
end if 

call getsurfxrt(sl,tt,xl,rl,tl) 
call getsurf xrt (is , t2 ,x2,r2, t2) 
call getsurf xrt (s3, tt ,x3 ,r3,t3) 
call ge t surf xrt (ss , t4,x4 ,r4,t4) 

dels “ delsp + deism 

c if (isurf.eq.2) write (lout , ♦) 'dels .dels .delsp, deism 

dxds - 0 . Be(x3 - xl) /dels 
drds - 0 . 6* (r3~rl ) /dels 

delt ■ deltp + delta 

c if (isurf.eq.2) write(lout ,♦) ’delt -’,delt 

dxdt • 0 . 5* (x2-x4) / de It 
drdt - 0. B*(r2-r4)/delt 

return 

end 


subroutine getsurf xrt (ss,tt,xs,rs,ts) 

c parameter (axdim- 1100) 

c real t(-6:axdia) ,s(-6;axdim) * 

c w(0:axdia,0 : axdim) , 
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c x(0:mxdi»,0:mxdim),y(0:mxdim,0:mxdim). 

C x(0:mxdia,0 :mxdim) ,u(0:l) ,v(0: 1) 

c common/surf data/ s.t.w.x.y.z 

c c omagra /par ami im/ u, v 

c comaion/nurb*!*#/ kt,k2 

call gat surf xyz(ss ,tt ,x» .ys.zs) 
ra ■ »qrt(y***2 + zs**2) 
ts “ atan2(x» ,ys) 

rat urn 
and 


subrout lna gatsurf xyx(ssin, tt ,xs .ys, zs) 
a 

♦** calculata point along tha NURBS curva **♦ 
a 

paramatar (mxdim-1100) 
raal t(-5:mxdim) ,s(-6 : nix dim) , 
v(0:mxdim,0:mxdim) , 

x (0 :mxdin , 0 : mxdim) , y (0: mxdim, 0: mxdim) , 
s { 0 : mxdim , 0 : mxdim) ,u(0:t) ,v(0:l) , 
rN0(0:mxdim,2) ,rNl (0:mxdim,2) , 
rN2(0:mxdim,2) ,rN3(0:mxdim,2) 

common /surf data/ s.t.V.X.y.z 
common/ par ami im/ u , v 
common/nurbsixa/ kl,k2 

c ommon/ if lags/ i saarch , idabug , 1 in , lout , lgrid 

aaa account for parlodicity in s *** 

ss - sain 

if (ssin.lt .0. ) ss - u(l) + ssin 
if (ssin.ga.u(l)) ss - ssin - u(l) 

if (tt.lt.T (0) ) tt - v(0) 
if (tt.ga.v(l)) tt * v( 1) -0 . 0001 

aaa calculata basis function: bi(s) *** 
do 1810 ii-0,kl 
*• initializa rN3 to 0 *♦ 
rN3(il,l) - 0.0 
rN2(ii,l) - 0.0 
rNl(ii,l) - 0.0 
rNO ( i i , 1 ) - 0.0 

if (ss.ga .s(ii) . and.ss .It .s(ii+l)) than 
rNO(ii.l) -1-0 
imid - ii 
and if 
1810 continua 


do 1811 i i- imid- 1, imid 
do 1811 ii-O.kl 
tarml - 0.0 

if (rN0(ii,l) .na. 0.0. and. 

(s(ii+l)-s(li)).na.0.0) 

tarml - (ss-s(ii>) / (s(ii+l)-s (ii)) * rNO(ii.l) 
tarm2 * 0.0 

if (rN0(ii+l,i) .na. 0.0. and. 

(s(ii+2)-s(U+l)>.na.O.) 

tarm2 - (s(ii+2)-ss) / (s(U*2)-s(ii+l) ) * rN0(ii+l,l) 

rNl(ii.l) - tarml + tarm2 

1811 continua 

do 1812 ii-iaid-2,imid 
do 1812 ii-O.kl 
tarml - 0.0 

if (rNl (ii, 1) .na .0.0. and. 

(s(ii+2)-s(ii)).na.0.) 

tarml - (ss-s(ii)) / (s(ii*2)-s(ii)) * rNl(ii.l) 
tarx2 - 0.0 

if (rNl(ii+l,l) .na. 0.0. and. 

(s(ii+3)-s(ii+l)).na.0.) 

tarm2 • (s(ii+3)-ss) / (s(ii+3)-s(ii+l) ) ♦ rNl(ii*l,l) 

rN2(ii,l) - tarml + tarm2 

1812 continua 

do 1813 ii“imid-3 , imid 
c do 1813 ii-O.kl 

tarml - 0.0 

if (rN2(ii,l) .na .0 .0. and. 

(s(ii+3)-s(ii)>.na.0.) 

tarml - (ss-s(ii)) / (*<ii+3)-.<il) ) * rN2(ii,l) 
tarm2 • 0.0 
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if (rN2(ii+l,l) .ne .O.O.and. 

(s(ii+4)-s(ii+l)) . ne .0 . ) 

t*rm2 - (s(ii+4)-ss) / (s(ii+4)-s(ii+l) ) * rN2(il+l,l) 

rN3(ii,l) ■ terml + term2 
1813 continue 


*♦* calculate basis function: bj(t) 
do 1820 j j *0 , k2 
rN3(J j ,2) - 0.0 
rN2(j j ,2) - 0.0 
rNl(jj ,2) - 0.0 
rN0(jj,2) - 0.0 

if (tt.ge.t(jj) . and .tt.lt,t(jj+l)) then 
rN0(jj,2) - 1.0 
jaid - jj 
end if 

1820 continue 

c if (jmid.lt. ijdel) jmid - ijdel 

c if (jmid.gt.(k2-ijdel)) jmid - k2-ijdel 

do 1821 jj-jmid-l.jmid 
c do 1821 jj-0,k2 

terml - 0.0 

if (rN0(jj, 2). ne.0.0. and. 

(t(jj+l)-t(jj)).ne.O.) 

terml - (tt-t(jj)) / (t (j j+l)-t ( j j) ) * rN0(jj,2) 
term2 “ 0.0 

if (rN0(jj+l,2) .ne.0.0. and . 

(t(jj+2)-t(jj+l)).ne,0.) 

term2 - (t(jj+2)-tt) / (t (j j+2)-t(j j+1) ) * rN0(jj+l,2) 

rNKjJ»2) ■ terml + term2 

1821 continue 

do 1822 j j*Jmid-2, jmid 
c do 1822 jj-0,k2 

terml ■ 0.0 

if (rfTUjj, 2). ne.0.0. and. 

(t(jj+2)-t(jj)).ne.O.) 

terml - (tt-t(jj)) / (t{j j+2)-t(j j ) ) ♦ rNl(jj,2) 
term2 * 0.0 

if (rNl (J j+1 ,2) .ne.0.0. and. 

(t(jj+3)-t(jj+l)).ne.O.) 

term2 - (t(jj+3)-tt) / (t ( j j+3)-t (j j+1) ) * rNl(jj+l,2) 

rN2(jj,2) * terml + term2 

1822 continue 

do 1823 jj*jmid-3, jmid 
c do 1823 jj-0,k2 

terml “ 0.0 

if (rN2(jj, 2). ne.0.0. and. 

(t(jj+3)-t(jj)).ne.0.) 

terml - (tt-t(jj)) / (t(j j+3)-t(j j) ) * rN2(jj,2) 
term2 - 0.0 

if (rN2(jj+l, 2). ne.0.0. and. 

(t(jj+4)-t(jj+l)).ne.0.) 

term2 - (t(jj+4)-tt) / (t(j j+4)-t(jj+l)) * rK2(jj+l,2) 

rN3(jj,2) * terml + term2 

1823 continue 


suml-Q . 0 
sum2*0 . 0 
sum3-0 . 0 
eum4“0.0 

do 1830 i-imid-3,imid 
do 1830 J»jmid-3, jmid 

bis - rN3(i, 1) 
bjt - rN3 ( J » 2 ) 

if (bis .ne . 0 . O.and .bjt .ne . 0 . 0) then 
suml - suml+v(i, j)*x(i, j)*bis+bjt 
svbb 2 - sua2+w ( i » j ) +y ( i # j ) *bl e +b j t 
sum3 - sum3+v(l,j)*z(i,j)+bis*bjt 
sum4 » sum4+w(i, j)*bisebjt 
end if 
1830 cont inue 

xs - suml/sum4 
ys - sum2/sum4 
xs - sum3/sum4 


return 
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•ad 


subroutine r*adig*s(f ilnm) 

p*r*m*t*r (mxdat a*100000 ,«utdim“l 100) 
r**l data(mxdata) ,t (-5:mxdim) ,*(-5 :mxdim) , 
w(0:mxdim,0 ;aucdim) , 

*(0:mxdim,0:mxdia) ,y(0:mxdim,0:mxdim) . 
s(0:mxdim,0:nudia) ,u(0: 1) ,v(0: 1) 

charact«r*80 filnm 

charact*r*70 chlin*, pars*. cht**p, blank 
charact*r*70 alio* 

common/ »urf data/ *,t,v,x,y,z 
c cannon/ par ami im/ u. v 
cozmon/mirbsiz*/ kl,k2 

blank(l:) - ’ 
blank (41 ; ) - 1 

op*n(unit* 2 ,fil*-f ilnm, status- 'unknown') 


♦♦* initializ* data array *** 

* 

e writ# (lout,*) 'Initializing data array ... ’ 

do 50 n-l,mxdata 
data(n) - -9999.0 
50 coat inu« 

* 

*** r*ad ICES film into data array *** 

* 

if ound ■ 0 

90 r*ad(2,99,*nd-100) alin* 

99 format (a70) 

if (alin«( :4) .*q. '128, ’ .and. ifound. *q.0) th*n 
print*, 'ICES typo 126 found... 1 
ifound " 1 
ityp* *• 126 
alin* * alin* (5:) 
idata ■ 1 

•Is* if (if ound.*q.0) th*n 
go to 90 
*nd if 

* 

*** par** out data from alin* *** 

♦ 

95 cht*mp * pars* (alin*) 

if (cht*sp.n* .blank) th*n 

r*ad(cht*mp, * ,*rr*100) data(idata) 
idata - idata * 1 

•Is* 

go to 90 
•ndif 
go to 95 

c r*ad(2,*,*rr-100) ityp* , (data(i) ,i-l,mzdata) 

100 clos*(2) 

writ* (lout ,*) 'ICES typ* *',ityp* 

c writ* (lout ,*) ' — — — 

ndata “ 0 

do 200 n*l,mxdata 

if (data(n) .n*. -9999.0) ndata - n 
200 cont inu* 

writ* (lout , *) ’ndata * ’ .ndata .data (ndata) 

If (ityp*. *q. 126) th*n 


*** ICES *ntity 128: B-Splin* Surfac* *** 


writ* (lout ,*) 'ICES *ntity 128: B-Splln* Surfac*' 

•*• assign variables from data array *** 

* 

iptr • 1 

kl - intCdata(iptr)) 
k2 * int(4ata(iptr+D) 

ml *> int(data(iptr+2)) 
m2 ■ int (data(iptr+3) ) 
ipropl - int (data(iptr+4) ) 
iprop2 ■ int (data (iptr+6) ) 
iprop3 - int ( dat a ( iptr+6 ) ) 
iprop4 * int (data (iptr+7)) 
iprop5 “ int (data (iptr+8)) 

nl • kl-ml+1 
n2 - k2-m2+l 
ia - nl+2*ml 
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ib - n2+2*m2 

ic - (kl+1) *(k2+l) 

writ* (lout ,*) *nl J n2 J a,b,e - ' ,d 1 ,n2,ia,ib, ic 

**+ knot points *** 

writ* (lout ,*) 'knot points’ 
iptr » 10 
do 1500 i*0, ia 

s(i) ■ data(iptr+i) 
c vrite(4,*) ’•«’,*(!) 

1500 continu* 

iptr ■ 11+ia 
do 1502 i-0,ib 

t(i) ■ data(iptr+i) 
c writ*(4,*) *t»',t(i) 

1602 continu* 

*♦* weights *** 

writ*(lout,*) ’weights’ 
iptr - 12+ia+ib 
do 1600 j*0,k2 
do 1600 i«0,kl 
w(i,j) - data(iptr) 

c writ* (4,*) 'i,j,w-M.j,w(i.J) 

iptr - iptr + 1 
1600 continu* 

*** control points *** 

writ* (lout , *) ’control pts ’ 
iptr - 12+ia+ib+ic 
do 1700 j«0,k2 
do 1700 i-0,kl 

*(i,j) - data(iptr) 
y(i,J) ■ data(iptr+l) 

*(i,J) ■ data(iptr+2) 

c writ*(4,*) i, j ,x(i, j) ,y(i, j) ,*(i,j) 

iptr * iptr + 3 
1700 continu* 

iptr - 12+ia+ib+4*ic 
u(0) ■ data(iptr) 

u(l) ■ data(iptr+l) 
v(0) ■ data(iptr+2) 

v(l) ■ data(iptr+3) 

write (lout ,*) 'End of data at: '.iptr* 3 

writ* (lout , *) ’uO.ul.vO.vi - ’ ,u(0) ,u(l) ,v(0) ,v(l) 

•Is* 

writ* (lout , *) ’ICES *ntity ’.itype,’ not supported.’ 
writ* (lout ,*) ’(fust b* surface (128)’ 
writ* (lout ,*) ’Program terminating...’ 
stop 
•nd if 

return 

•nd 


function xofs(ss) 

*** axial location for ss valu* at constant tt value ♦*• 
common/ ttval/ tt 

call got surf xyz(ss,tt ,xs ,ys , ss) 
xofs * xs 

return 

•nd 


function xofsm(ss) 
♦** negative of xofs *** 
common/ ttval/ tt 

xofsa - -xofs (ss ) 

return 

•nd 


subroutine letef ind(njpts) 
parameter (mxdim* 1 100) 

real t(-5:mxdia) ,*(-5:axdi») ,w(0:mxdim,0:mxdim) , 

• x (0: mxdim, 0: mxdim) ,y(0:mxdim,0;mxdim) , 

x(0:mxdim ( 0:mxdim) ,u(0:l) ,v(0:l), 

. xt* (mxdim) , yte (mxdim) ,zt* (mxdim) ,rte (mxdim) , 

xl* (mxdim) ,yle (mxdim) ,xl* (mxdim) ,rl* (mxdim) , 
ale (mxdim), ste (mxdim) ,tl* (mxdim) ,tte (mxdim) 

r*al csl* (mxdim) ,cst* (mxdim) , cxl* (mxdim) , ext e (mxdim) , 
. cyle (mxdim) , cyt* (mxdim) ,czle (mxdim) , ext* (mxdim) 
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raal at(15) ,xt(15) 


common/ surf data/ a.t.v.x.y.z 
common/ par ami ia/ u.v 
comaon/nurb* iza/ kl • k2 

common/ xyzadga a/ xla,xt*,yl*,yta,xla,zta,rla,rta. 

cxla,cxta ,cyla ,cyta,czla,czta 

co«»on/ttval/ tt 
axtarnal xofa.xofam 

op#n(unit«13,f ila- * la .pta ' .atntua- 'unknown * ) 
opan(unit-14 ,f ila* ' t« .pts 1 .xtatua* 'unknown* ) 
writa(13,*) njpta 
vrita(14 ,♦) njpta 

* 

aaa Find adgaa 
a 

tol - 0.001 

dalta2 - ((v(l)-O.OOOD-v(O)) / raal(njpta-l) 
do 1760 nj-1, njpta 

c writ* (lout ,*) 'Finding min/max for: ^nj,* of njpta 

tt - v(0) ♦ dal ta2araal (nj-1) 
if (tt.aq.O.) tt - l.a-4 

*** find trailing (max I) valua for tt 

aa brackat max X point *a 

do 100 i-1,13 

fact ■ raal(i-2) /10 . 
at(i)*u(0)+fact*(u(l)-u(0)) 
xt (i)*xofa(at(i) ) 

100 continu* 

do 110 i*3, 13 

if (xt(i-l) .gt.xt(i-2) -and. 
xt(i-l) .gt . xt (i) ) than 
al - at (i-2) 
a2 - at(i-l) 
a3 - *t(i) 
and if 
1 10 coot inua 

xmax-brant(al , a2 , a3 . xof am .tol . aaax) 

call gataurfxyz(amax,tt,xa »ya,za) 

xta(nj) ■ xa 
yta(nj) - y* 

zta(nj) • xa 

rta(nj) - aqrt(ya**2+za**2) 
ata(nj) ■ »aax 
tta(nj) ■ tt 

aaa find landing (min X) valua for tt *** 
do 120 i*3,13 

if (xt(i-l) .lt.xt(l-2) .and. 
xt(i-l) -lt.xt(i)) than 
al - at(l-2) 
a2 - at(i-l) 
a3 ■ at(i) 
and if 
120 contlnua 

xmin-brant(al,a2,a3,xofa,tol,amin) 
call gataurfxyz(*min,tt,xa,ya,za) 

xla(nj) ■ xa 
yla(nj) - y» 

zla(nj) - xa 

rla(nj) - aqrt (yaaa2+za**2) 
ala(nj) « aain 
tla(nj) * tt 

vrita(13,997) xla(nj) ,rla(nj ) ,0.0 
vrita(14,997) xta(nj) ,rta(nj ) ,0.0 
097 format (lx, 4 (f 12 .6, lx)) 

aa* naxt tt valua *** 

1760 continua 

cloaa(13) 

cloaa(14) 

rat urn 
and 
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spli.na routinas 

* 

subrout ina spcoaf (spt , spf , coaf t , nsp) 
raal coaf f (*), spt (*), spf (*) , 

b(200) ,h(200) ,uu(200) ,w(200) 

* spt (200) - t 

* spf (200) - y 

* coaff (200) - z 

do i-1, nsp-1 

h(i)-spt (i+l)-spt (i) 
b(i)-6.0*(spf (i+l)-spf (i))/h(i) 
and do 

uu(l)»2.0*(h(l)*h(2)) 

w(l)-b(2)-b(l) 

do i«2, nap-2 

uu(i)-2.0*(h(i)+h(i-l))-h(i-l)**2/uu(i-l) 
w(i)-b(i)-b(i-l)-h(i-l)*w(i-l)/uu(i-l) 
and do 

coaff (nsp-1 )-0.0 
do i-nsp-2,2,-1 

coaff (i)-(w(i)-h(i)*coaff (i+l))/uu(i) 
and do 

coaff (l)-O.O 

rat urn 
and 


function spaval(xx ,spt , spf , coaf f .nsp) 
raal coaff (*), spt (*), spf (*) 

♦ spt (200) » t 

♦ spf (200) - y 

♦ coaff (200) - z 

if (xx . It . spt (1) ) than 
C wita(lout ,*) xx , spt (1 ) ,nsp 

c vrita(lout ,*) 'Warning... outsida (for) orig data — axtrapolating' 

h-spt(2)-»pt(l) 
a-(coaff(2)-coaff(l))/6.0/h 
b-coaff (l)/2.0 

c— -h/6 . 0* coaff (2) -h/3 . 0* coaff (1 ) ♦ 

(spf (2)-spf (1) )/h 
dx-xx-*pt (1) 

spaval-spf (1) +dx*(c+dx*(b+dx*a) ) 
rat urn 
and if 

do i«l,nsp-l 

if (xx.ga.spt(i).and.xx.la.spt(i+l)) than 
h-spt(i+l)-spt (i) 
a- (coaff (i+1) -coaff (i))/6.0/h 
b-coaff (i)/2.0 

c— h/6.0*coaff (i+l)-h/3.0*coaff(i) + 

(spf (i-H)-spf (i))/h 
dx-xx-spt(i) 

spaval-spf ( i) +dx* (c+dx* (b+dx*a) ) 
rat urn 
and if 
and do 

if (xx.gt . spt (nsp)) than 
c vrlta(lout,*) xx, spt (nsp), nsp 

c vrita(lout ( *) 'Warning... outsida (aft) orig data — axtrapolating’ 

h-spt(nsp)-spt (nsp-1) 
a-(coaff (nsp) -coaf f (nsp-1)) /6 . 0/h 
b-coaff (nsp-l)/2.0 
c— h/6 . 0*coaf f (nsp) -h/3 .0* 

coaff (nsp-1) +(spf (nsp) -spf (nsp-1 ))/h 
dx-xx-spt (nsp- 1 ) 

spaval-spf (nsp-1 ) +dx* (c+dx* (b+dx*a)) 
ratum 
and If 

*** arror trap **• 

vrita(lout ,*) 'Error occurrad in SPEVAL rotuina!’ 
vrita(lout,*) 'Talus of splina indax - ’ ,xx 
stop 

and 
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FUNCTION BRENT( AX ,BX,CX, F ,TOL, XMIN) 

PARAMETER (ITMAX-200.CCOLD- . 3B19660 .ZEPS-l .OE-IO) 
raal brant , ax.bx.cx ,tol,xmin,f 
axtarnal { 

A-HIN(AX.CX) 

B-KAX(AX.CX) 

V-BX 

V-V 

X-V 

E-O. 

FX*F(X) 

FV-FX 

FW-FX 

DO 11 ITER-1, ITHAX 
IM-0.6*(A+B) 

TOLl-TOL*ABS (I) +ZEPS 
T0L2-2.+T0L1 

IF(ABS(X-XM).LE. (T0L2- .5*(B-A) )) COTO 3 
IF(ABS(E) .CT.TOll ) THEN 
R-(X-V)*(FX-FV) 

Q- (X-V) ♦ (FX-FV) 

P-(X-V)*Q-(X-W)*R 

Q»2.*(Q-R) 

IF(Q.CT.O.) P— P 
Q-ABS(Q) 

ETEMP-E 

E-D 

IF(ABS(P> .CE.A3SC .6*Q*ETEKP) .Ofi.P.LE.Q*(A-X) .OR. 

♦ P.CE.Q+(B-X)) COTO 1 

D-P/Q 
U-X+D 

IF(U-A.LT.T0L2 .OR. B-U.LT.T0L2) D-SICN(TOLl . XM-X) 
COTO 2 

END IF 

1 IF(X.CE.XM) THEN 

E-A-X 

ELSE 

E-B-X 

ENDIF 

D-CC0LD*E 

2 IF(ABS(D) .CE.TOL1) THEN 

U-X+D 

ELSE 

U-X+5ICN(TQL1 ,D) 

ENDIF 

FU-FCU) 

IF(FU.LE.FX) THEN 
IF(U.CE.X) THEN 
A-X 
ELSE 
B-X 
ENDIF 
V-W 
FV-FW 
W-I 
FW-FX 
X-U 
FX-FU 
ELSE 

IF(U.LT.X) THEN 

A-U 

ELSE 

B-U 

ENDIF 

IF(FU.LE. FV .OR. W.EQ.X) THEN 

V-W 

FV-FW 

w-u 

FW-FU 

ELSE IF(FU.LZ.FV .OR. V.EQ.X .OR. V.EQ.W) THEN 

v-u 

FV-FU 

ENDIF 

ENDIF 

11 CONTINUE 

PAUSE ’Brant axcaad maximum iteration*.' 

3 XMIN-X 
BRENT-FX 
RETURN 
END 


SUBROUTINE VINOKUR (S.LMAX.SMIN.SMAX.DSAE.DSBE.KASE.IERR) 

C 

C VI NX ITER - CONTROL ROUTINE TO SATISFY STRETCHING CONSTRAINTS IN 
C VINOKUR’S FUNCTION EXACTLY 

C 

REAL S(«) 

TOLHIN - 1.0E-6 
T0LMIN2 - 1.0E-6 
I ERR - 0 
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c 

C VINOKUR’S FUNCTION CREATES A DISTRIBUTION OF GRID POINTS WHICH 
C SATISFY A SPECIFIED DERIVATIVE CONDITION, BUT WE REQUIRE A DELTA-S 
C CONSTRAINT INSTEAD. THESE TWO VALUES ARE EQUIVALENT ONLY TO FIRST 
C ORDER, AND HENCE, WE RESORT TO AN ITERATIVE PROCEDURE TO OBTAIN 
C HORE ACCURATE DELTA-S’ S. UP TO TEN ITERATIVE SWEEPS ARE 
C MADE. THE FIRST GUESS SETS DS/DXI - DELTA-S. THE NEXT CUESS 
C RECALCULATES DS/DXI USINC THE LEADING TERM IN THE TRUNCATION ERROR 
C (D2S/D(XI)**2) . THE NEXT EIGHT ITERATIONS USE A 2-D SECANT 
C ALGORITHM TO HOME IN ON THE DS/DXI ’S AT BOTH ENDS WHICH WILL CIVE 
C THE CORRECT DELTA-S. 

C 

C IN THE CASES WHERE A SINGLE-SIDED STRETCHING FUNCTION IS REQUIRED, 

C (KASE - 1 OR 2) A SECANT ALGORITHM IN 1-D IS APPLIED INSTEAD 
C 

C KASE - 0 

C STRETCHING ON BOTH ENDS, SO USE A 2-D SECANT METHOD TO ARRIVE AT THE 

C VALUES OF DSA AND DSB WHICH WILL SATISFY DS1E AND DS2E WITHIN 

C ROUNDOFF. 

C 

IF (KASE .EQ. 0 )THEN 
C 

c INITIAL CUESS - AN - DSAE, BN - DSBE 

AN2 - DSAE 
BN2 - DSBE 
C 

CALL VINK (S,LMAX,SMIN.SMAX,AN2,BN2, ESA, ESB, KASE ) 

C 

FN2 - ESA/DSAE - 1 
CN2 - ESB/DSBE - 1 
C 

C SECOND CUESS - CALCULATE DS1 AND DS2 FROM A TRUNCATED TAYLOR SERIES 

DSSA - 2.*S< l)-6. *S( 2)+4.*S( 3) -S( 4) 

DSSB - 2 . ♦S(LMAX) -5 . *S (LMAX-1) +4 . *S (LMAX-2) -S(LHAX-3) 

AN1 - DSAE-0.5*DSSA 
BN1 « DSBE+O . 6*DSSB 


CALL VINK (S, LMAX, SMIN, SMAX, AN1.BN1, ESA. ESB, KASE ) 

C 

FN1 - ESA/DSAE - 1 

GN1 - ESB/DSBE - 1 

AN - AN1 

BN - BN1 

C 

C 3RD THRU 10TH GUESSES , USE 2-D SECANT METHOD 

DO N - 3,20 
C 

c CALCULATE OFFSET DERIVATIVES 

CALL VINK (S.LMAX, SKIN, SMAX, AN2, BN1 ,ESA21 ,ESB21 ,KASE ) 

CALL VINX(S, LMAX, SHIN, SHAX, AN1 ,BN2,ESA12, ESB12, KASE } 

FA - ( ESA - ESA21 )/DSAE 

FB - ( ESA - ESA12 )/DSAE 

GA - ( ESB - ESB21 )/DSBE 

CB - ( ESB - ESB 12 )/DSBE 

DEN - FA*CB - F3*CA 

DELA - -(AN1 - AN2) 

DELB - “(BN1 - BN2) 

C 

C STICK WITH LAST CUESS IF APPROACHINC ROUNDOFF 

IF (ABS(DEN) .LT.T0LMIN2) THEN 

CALL VINK (S.LMAX, SHIN, SMAX, AN, BN, ESA. ESB. KASE ) 
RETURN 
END IF 
C 

c CALCULATE NEXT DISTRIBUTION 

AN • AN1 + DELA*( CB*FN1 - FB*GN1 )/DEN 
BN - BN1 + DELB* (- CA*FN1 + FA*CN1 )/DEN 


CALL VINK(S. LMAX, SMIN, SMAX, AN.BN.ESA, ESB, KASE ) 

FN - ESA/DSAE - 1 
CN - ESB/DSBE - 1 

UPDATE N, N-l, N-2 AND CONTINUE 

AN2 - AN1 

BN2 - BN1 

AN1 - AN 

BN1 - BN 

FN1 - FN 

CN1 - CN 


C NEXT CUESS 

ENDDO 

C 

C KASE - 1 

C STRETCHINC AT THE LAST ENDPOINT ONLY, SO USE A i-D SECANT METHOD 
C TO ARRIVE AT THE VALUES OF DSB WHICH WILL SATISFY DSBE WITHIN 
C ROUNDOFF. 

C 

ELSEIF (KASE.EQ. 1) THEN 
C 

c INITIAL CUESS - BN - DSBE 
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c 

c. 


c 

c. 


c 

c. 


c 

c. 


c 

c. 


c 

c. 

c 

c. 

c 

c 

c 

c 

c 

c. 


c 

c, 


c 

c 


c 

c. 


c 

c. 


c 

c. 


c 

c. 

c 

c. 

c 

c. 


c 

c. 


c 

c- 

c 


BN2 - DSBE 

CALL VINK(S , LMAX, SMI N, SMAX, DS AE ,BN2, ESA, ESB, XASE ) 

FN2 - ESB/DSBE - 1 

SECOND CUESS - CALCULATE DS1 AND DS2 FROM A TRUNCATED TATLOR SERIES 
DSSB - 2 . *S(LMAX) -6 . •SCLMAX-l ) +4 . *S (LMAX-2) -S(LMAX-3) 

3N1 - DSBE-0.5*DSSB 

CALL VINX(S, LMAX, SHIN. SMAX, DSAE, BN1, ESA, ESB .KASE ) 

FN1 ■ ESB/DSBE - 1 
BN - BN1 

. . 3RD THRU 10TH CUESSES , USE 1-D SECANT METHOD 
DO N - 3,20 

STICK WITH LAST CUESS IF APPROACH INC ROUNDOFF 
DEN - FN1-FN2 

IF (ABS(DEN) .LT.T0LMIN2) THEN 

CALL VINK(S, LMAX, SMIN, SMAX, DSAE, BN, ESA, ESB. KASE ) 

RETURN 
END IF 

CALCULATE NEXT DISTRIBUTION 

BN - BN1 - FNi* (BN1-BN2) /DEN 

CALL VINK(S ,LMAX,SMIN , SMAX ,DSAE ,BN ,E5A ,ESB , KASE ) 

FN - ESB/DSBE - 1 

UPDATE N, N-l, N-2 AND CONTINUE 

BN2 - BN1 
BN1 - BN 
FN2 - FN1 
FN1 - FN 

.. NEXT CUESS 
END DO 

KASE ~ 2 

STRETCHING AT THE FIRST ENDPOINT ONLY, SO USE A 1-D SECANT METHOD 
TO ARRIVE AT THE VALUES OF DSA WHICH WILL SATISFY DS1E WITHIN 
ROUNDOFF. 

ELSEIF (KASE.EQ.2) THEN 

. . . INITIAL CUESS - AN - DSAE 
AN2 - DSAE 

CALL VINK(S,LMAX,SMIN, SMAX, AN2, DSBE. ESA. ESB, KASE ) 

FN2 - ESA/DSAE - 1 

SECOND CUESS - CALCULATE DS1 AND DS2 FROM A TRUNCATED TAYLOR SERIES 
DSSA - 2.*S( 1)-5.»S( 2)+4.*S( 3) -S( 4) 

AN1 - DSAE-0.6*DSSA 

CALL VINK (S.LMAX, SHIN, SMAX, AN 1, DSBE, ESA. ESB, KASE ) 

FN1 - ESA/DSAE - 1 

AN - AN1 

. . . 3RD THRU 10TH CUESSES , USE 1-D SECANT METHOD 
DO N - 3,20 

STICK WITH LAST CUESS IF APPROACHING ROUNDOFF 

DEN - FN1-FN2 

IF (ABS(DEN) .LT.T0LMIN2) THEN 
CALL VINKtS.LMAX, SHIN, SMAX, AN, DSBE, ESA, ESB, KASE ) 

RETURN 

ENDIF 

CALCULATE NEXT DISTRIBUTION 

AN - AN1 - FN1*(AN1-AN2)/DEN 

CALL VINK (S.LMAX, SMIN, SMAX, AN, DSBE, ESA. ESB, KASE ) 

FN - ESA/DSAE - 1 

UPDATE N, N-l, N-2 AND CONTINUE 

AN2 - AN1 

AN1 - AN 
m - FN1 
FN1 - FN 

NEXT CUESS 

ENDDO 

. END CASE TEST 
END IF 

. ERROR TEST 
IERR - 1 

WRITE(* ,*) ’ERROR IN VINOKUR’ 

WRITE(*,*> KASE.ABS(DEN) 

, RETURN TO CALLING ROUTINE 
RETURN 
END 
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C VINX - STRETCHES POINTS ON A LINE SO THAT SPECIFIED DERIVATIVES 
C AT THE EDCES ARE SATISFIED (TAKEN FROM NASA CR 3313 BY 
C VINOKUR (1980)) 

C 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

c 

c 

c 


SUBROUTINE VINX(S,LMAX.SMIN,SHAX,DS1 .DS2.ES1 .ES2.KASE ) 

VINOKUR’S ALGORITHM IS DESICNED TO DISTRIBUTE A SPECIFIED 
NUMBER OF POINTS ALONC A CURVE, CIVEN THE NUMBER OF POINTS, 

THE LENGTH OF THE CURVE, AND THE DERIVATIVE CONDITIONS AT 
BOTH ENDS OF THE CURVE. IN CFD APPLICATIONS, THE USER WOULD 
USUALLY RATHER SPECIFY THE DELTA-S’S AT THE ENDS OF THE CURVE, 
WHICH ARE EQUIVALENT TO THE DERIVATIVES ONLY TO FIRST ORDER. 
THEREFORE. THE USER MAY WISH TO APPLY THIS ALGORITHM ITERATIVELY 
TO OBTAIN AN EXACT DELTA-S SPECIFICATION. SUBROUTINE 
VINXITER WILL ITERATE ON THIS SCHEME UNTIL THE PROPER DELTA-S 
CONSTRAINTS ARE SATISFIED. 


INPUT: 

LMAX - NUMBER OF POINTS ON THE CURVE 
SMIN, SMAX - BECINNINC AND END VALUES OF S 
DS1, DS2 - THE DERIVATIVE END CONDITIONS INPUT INTO 
VINOKUR’S FUNCTION 

KASE - 0 - SATISFY DELTA-S ON BOTH ENDS 

- 1 - SATISFY DELTA-S ONLY AT XI -UMAX 

- 2 - SATISFY DELTA-S ONLY AT XI-XIMIN 

OUTPUT: 

S( XI ) - RESULTING S DISTRIBUTION FROM VINOKUR’S FUNCTION 

Kl - ( S (XIMIN+1)-S (XIMIN) ) <- CALCULATED DELTA-S 

ES2 - ( S(XIMAX) -S(XIMAI-l) ) <- 


ADDITIONALLY, THIS VERSION USES THE APPROXIMATE INVERSE SOLUTION 
FOR Y-SIN(X)/X AND Y*SINH(X)/X RATHER THAN A NEWTON ITERATION. THE 
APPROXIMATE SOLUTION WAS ALSO TAKEN FROM NASA CR 3313. 


COMMON /PIVAL/ PI 

DIMENSION S(*), D1(4,2),D2(4,2) 


C 

C CALCULATE CONSTANTS 

SDEL - SMAX-SMIN 
SO-SDEL/FLOAT (LMAX- 1 ) /DS 1 
S1»SDEL/FLQAT(LHAX-1)/DS2 
B-SQRT(S0*S1) 
A«SqRT(SO/Sl) 
c 

C USE VARIUOS KASE TYPE 

IF (KASE.EQ.l) THEN 
B-Sl 

ELSEIF (KASE.EQ .2} THEN 
B-SO 
END IF 
C 


C CALCULATE X BASED ON VALUE OF B 

IF ( B.LT.1.0 ) THEN 
C 

c X IS REAL 


IF ( B . LT . 0 . 26938972 ) THEN 

X - PI*(1. -B 4 B**2 - (1.4PI**2/6.)*B**3 

> + 6.794732*B**4 -13.205501*B**5 4 11 .726095*B**6) 

ELSE 

C -1,-B 

X - SQRT(6. *C)*(1. 

> 40.16*C 4 0 . 057321429*0**2 +0 . 048774238*C**3 

> -0.053337753*0**4 4 0 . 075845 134*C**5) 

ENDIF 

C 


C X IS ZERO 

ELSEIF ( B.EQ.1.0 ) THEN 
X-0. 

C 

c X IS IMAGINARY 


ELSE 

IF (B.LT. 2. 7829681) THEN 
C - B-l. 

X - SQRT(6. *C)*(1 . 

* -0.16*C 4 0.0573214290*0**2 - 0.0249072950*0**3 

* + 0. 0077424461 *C**4 - 0.0010794123*0* *5) 
ELSE 

V - ALOC(B) 

W - l./B - 0.028527431 
X ■ V 4 (1.4l./V)*AL0C(2.*V) -0.02041793 

* 4 0. 24902722*W 4 1 . 9496443*W**2 

* - 2. 6294647* W* *3 4 8 .56795911*W**4 
ENDIF 

ENDIF 

C 

c DISTRIBUTE POINTS ALONC EDCE 

IF ( KASE.EQ.l .OR. KASE.EQ. 2 ) THEN 
S(1 ) - 0.0 

S(LMAX) - SDEL 
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DO I - 2.LHAX-1 
J - LHA1+1-I 
XI - FLOAT(I-l) ( (LHAX-l) 

IF (B.CT. l.OOOl) THEN 

VI - 1. + TANK(X/2.*(XI-l.))/TANHCl/2.) 

ELSEIF (B.LT. 0.9999) THEN 

U1 - 1. + TAN (X/2.*(XI-1.))/TAN (X/2.) 

ELSE 

U1 - XI*(1.-.6*(B-1.)*(1- -XI)* (2. -XI) ) 

END IF 

U2 - SINH(XI*X)/SINH(I) 

IF (KASE.EQ.l) THEN 
FACT - ABS(DSl) 

S(J) - ((1.-FACTW1.-U1) + FACT*(1.-U2) ) ♦SDEL 
ELSEIF (KASE.EQ.2) THEN 
FACT - ABS(DS2) 

S(I) - ( { 1 . -FACT) • U1 + FACT* U2 ) *SDEL 

END IF 
ENDDO 
C 

c KASE - 0 

ELSE 

DO I-l.LMAX 

XI-FLOAT(I-l) /FLOAT (LHAX-1) 

CNUM-I*(XI-0.6) 

CDEN-X/2. 

IF (B.LT. 0.9999) THEN 

CC - TAN(CNUH)/TAN (CDEN) 

U - 0.6*(1.+CC) 

ELSEIF (B. CE. 0.9999. AND. B.LE. 1.0001) THEN 
U»XI*(1.+2.*(B-1 . )*(XI-0.6)*(1.-XI)) 

ELSEIF (B.CT. 1.0001) THEN 
CC - TANH(CNUM)/TANH(CDEN) 

V - 0 . 5*(1 - +CC) 

END IF 

S(I) - V*SDEL/(A+(1.-A)*U) 

ENDDO 
END IF 
C 

DO L - l.LHJUC 

S(L) - S(L) ♦ SMIN 
ENDDO 

ESI - S( 2)-S( 1) 

ES2 - S(LMAX)-S (LKAX-1) 

C 

c RETURN TO CALL INC ROUTINE 

RETURN 

END 

charactor*70 function PARSE(a) 



C 

C parao part* out a aubatring 
C 

C 901017 cja 
C 

C 

character a*70, atiztp*70> c*l 
intogor alon, l, j, *1 
logical dolim 
c 

c > Strip loading vhitoapaco and dolimitora 

100 if ( ( .not . dolin(a(l:l))) .or. (alon(a) .oq.O) ) goto 101 

•tap ■ *(2;alon(s)) 

a * atxp 

goto 100 

101 continuo 
c 

C > Soarch for tho noxt dolizitor 

c 

al - alon(a) 

i « cl +1 
do 110 j - al ,2,-1 
c • a (J : j ) 

if (dolia(c)) i » j -1 
110 continuo 
c 

c > Rotum tho aubatring and roduco a 

c 

if (i.gt.O) than 
parao * a(l:i) 
atzp ■ a ( i+1 : lon(a) ) 

a * atop 

olaa 

parao ■ ’ ' 

•( 1 : 1 ) - ' ’ 
ondif 
rotum 
and 
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logical function DELIM(c) 


C 

C 

C Raturn trua if tha variabla c is a dalimitar 
C 

C cjm 901214 
C 

C 

charactar c*l 
intagar tstc(6), i 

data tstc/32, 9, 44, 47, 66, 69/ 

c 

dalim - .falsa, 
do 100 i - 1,6 

if ( ichar(c) .aq. tstc(i) ) dalim ■ .trua. 

100 continua 
rat urn 

and 

intagar function SLEN(s) 

C 

C 

C raturn tha actual string langth (i.a., strip trailing spacas) 

C 

C Jch 901017 
C 

C 

charactar s*70 
intagar i, iara(70) 
logical dalim 
C 

C strip trailing whita spaca and daliaitars to find actual langth 
i ■ lan(s) +1 
100 continua 
i - i-1 

if ( (i.gt.0) .and. dalim(s(i :i>) ) goto 100 

slan " i 

raturn 

and 
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